Search NASA⌕ Search

SEARCH · Search NASA

Results for “Mass density ratio”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A sharp interface Lagrangian-Eulerian method for rigid-body fluid-structure interaction

Herein this paper introduces a sharp interface method to simulate fluid-structure interaction (FSI) involving rigid bodies immersed in viscous incompressible fluids. The capabilities of this methodology are benchmarked using a range of test cases and demonstrated using large-scale models of biomedical FSI. The numerical approach developed herein, which we refer to as an immersed Lagrangian-Eulerian (ILE) method, integrates aspects of partitioned and immersed FSI formulations by solving separate momentum equations for the fluid and solid subdomains, as in a partitioned formulation, while also using non-conforming discretizations of the dynamic fluid and structure regions, as in an immersed formulation. A simple Dirichlet-Neumann coupling scheme is used, in which the motion of the immersed solid is driven by fluid traction forces evaluated along the fluid-structure interface, and the motion of the fluid along that interface is constrained to match the solid velocity and thereby satisfy the no-slip condition. To develop a practical numerical method, we adopt a penalty approach that approximately imposes the no-slip condition along the fluid-structure interface. In the coupling strategy, a separate discretization of the fluid-structure interface is tethered to the volumetric solid mesh via stiff spring-like penalty forces. Our fluid-structure coupling scheme relies on an immersed interface method (IIM) for discrete geometries, which enables the accurate determination of both velocities and stresses along complex internal interfaces. Numerical methods for FSI can suffer from instabilities related to the added mass effect, but computational tests indicate that the methodology introduced here remains stable for selected test cases across a broad range of solid-fluid mass density ratios, including extremely small, nearly equal, equal, and large density ratios. Biomedical FSI demonstration cases include results obtained using this method to simulate the dynamics of a bileaflet mechanical heart valve in a pulse duplicator, and to model transport of blood clots in a patient-averaged anatomical model of the inferior vena cava.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A sharp interface Lagrangian-Eulerian method for flexible-body fluid-structure interaction

This paper introduces a sharp-interface approach to simulating fluid-structure interaction (FSI) involving flexible bodies described by general nonlinear material models and across a broad range of mass density ratios. This new flexible-body immersed Lagrangian-Eulerian (ILE) scheme extends our prior work on integrating partitioned and immersed approaches to rigid-body FSI. Our numerical approach incorporates the geometrical and domain solution flexibility of the immersed boundary (IB) method with an accuracy comparable to body-fitted approaches that sharply resolve flows and stresses up to the fluid-structure interface. Unlike many IB methods, our ILE formulation uses distinct momentum equations for the fluid and solid subregions with a Dirichlet-Neumann coupling strategy that connects fluid and solid subproblems through simple interface conditions. As in earlier work, we use approximate Lagrange multiplier forces to treat the kinematic interface conditions along the fluid-structure interface. This penalty approach simplifies the linear solvers needed by our formulation by introducing two representations of the fluid-structure interface, one that moves with the fluid and another that moves with the structure, that are connected by stiff springs. This approach also enables the use of multi-rate time stepping, which allows us to use different time step sizes for the fluid and structure subproblems. Our fluid solver relies on an immersed interface method (IIM) for discrete surfaces to impose stress jump conditions along complex interfaces while enabling the use of fast structured-grid solvers for the incompressible Navier-Stokes equations. The dynamics of the volumetric structural mesh are determined using a standard finite element approach to large-deformation nonlinear elasticity via a nearly incompressible solid mechanics formulation. This formulation also readily accommodates compressible structures with a constant total volume, and it can handle fully compressible solid structures for cases in which at least part of the solid boundary does not contact the incompressible fluid. Selected grid convergence studies demonstrate second-order convergence in volume conservation and in the pointwise discrepancies between corresponding positions of the two interface representations as well as between first and second-order convergence in the structural displacements. The time stepping scheme is also demonstrated to yield second-order convergence. To assess and validate the robustness and accuracy of the new algorithm, comparisons are made with computational and experimental FSI benchmarks. Test cases include both smooth and sharp geometries in various flow conditions. Furthermore, we also demonstrate the capabilities of this methodology by applying it to model the transport and capture of a geometrically realistic, deformable blood clot in an inferior vena cava filter.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

The Milky Way tomography with APOGEE: intrinsic density distribution and structure of mono-abundance populations

ABSTRACT The spatial distribution of mono-abundance populations (MAPs, selected in [Fe/H] and [Mg/Fe]) reflect the chemical and structural evolution in a galaxy and impose strong constraints on galaxy formation models. In this paper, we use APOGEE data to derive the intrinsic density distribution of MAPs in the Milky Way, after carefully considering the survey selection function. We find that a single exponential profile is not a sufficient description of the Milky Way’s disc. Both the individual MAPs and the integrated disc exhibit a broken radial density distribution; densities are relatively constant with radius in the inner Galaxy and rapidly decrease beyond the break radius. We fit the intrinsic density distribution as a function of radius and vertical height with a 2D density model that considers both a broken radial profile and radial variation of scale height (i.e. flaring). There is a large variety of structural parameters between different MAPs, indicative of strong structure evolution of the Milky Way. One surprising result is that high-α MAPs show the strongest flaring. The young, solar-abundance MAPs present the shortest scale height and least flaring, suggesting recent and ongoing star formation confined to the disc plane. Finally we derive the intrinsic density distribution and corresponding structural parameters of the chemically defined thin and thick discs. The chemical thick and thin discs have local surface mass densities of 5.62 ± 0.08 and 15.69 ± 0.32 M⊙pc−2, respectively, suggesting a massive thick disc with a local surface mass density ratio between thick to thin disc of 36 per cent.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

MaNGA DynPop – I. Quality-assessed stellar dynamical modelling from integral-field spectroscopy of 10K nearby galaxies: a catalogue of masses, mass-to-light ratios, density profiles, and dark matter

ABSTRACT This is the first paper in our series on the combined analysis of the Dynamics and stellar Population (DynPop) for the MaNGA survey in the final SDSS Data Release 17 (DR17). Here, we present a catalogue of dynamically determined quantities for over 10 000 nearby galaxies based on integral-field stellar kinematics from the MaNGA survey. The dynamical properties are extracted using the axisymmetric Jeans Anisotropic Modelling (JAM) method, which was previously shown to be the most accurate for this kind of study. We assess systematic uncertainties using eight dynamical models with different assumptions. We use two orientations of the velocity ellipsoid: either cylindrically aligned JAMcyl or spherically aligned JAMsph. We also make four assumptions for the models’ dark versus luminous matter distributions: (1) mass-follows-light, (2) free NFW dark halo, (3) cosmologically constrained NFW halo, (4) generalized NFW dark halo, i.e. with free inner slope. In this catalogue, we provide the quantities related to the mass distributions (e.g. the density slopes and enclosed mass within a sphere of a given radius for total mass, stellar mass, and dark matter mass components). We also provide the complete models which can be used to compute the full luminous and mass distribution of each galaxy. Additionally, we visually assess the qualities of the models to help with model selections. We estimate the observed scatter in the measured quantities which decreases as expected with improvements in quality. For the best data quality, we find a remarkable consistency of measured quantities between different models, highlighting the robustness of the results.

Astronomy & Astrophysics↗
DOE OSTI2022

A consistent and conservative Phase-Field method for multiphase incompressible flows

In the present study, a consistent and conservative Phase-Field method, including both the model and scheme, is developed for multiphase flows with an arbitrary number of immiscible and incompressible fluid phases. The consistency of mass conservation and the consistency of mass and momentum transport are implemented to address the issue of physically coupling the Phase-Field equation, which locates different phases, to the hydrodynamics. These two consistency conditions, as illustrated, provide the “optimal” coupling because (i) the new momentum equation resulting from them is Galilean invariant and implies the kinetic energy conservation, regardless of the details of the Phase-Field equation, and (ii) failures of satisfying the second law of thermodynamics or the consistency of reduction of the multiphase flow model only result from the same failures of the Phase-Field equation but are not due to the new momentum equation. Physical interpretation of the consistency conditions and their formulations are first provided, and general formulations that are obtained from the consistency conditions and independent of the interpretation of the velocity are summarized. Then, the present consistent and conservative multiphase flow model is completed by selecting a reduction consistent Phase-Field equation. Several novel techniques are developed to inherit the physical properties of the multiphase flows after discretization, including the gradient-based phase selection procedure, the momentum conservative method for the surface force, and the general theorems to preserve the consistency conditions on the discrete level. Equipped with those novel techniques, a consistent and conservative scheme for the present multiphase flow model is developed and analyzed. The scheme satisfies the consistency conditions, conserves the mass and momentum, and assures the summation of the volume fractions to be unity, on the fully discrete level and for an arbitrary number of phases. All those properties are numerically validated. Finally, numerical applications demonstrate that the present model and scheme are robust and effective in studying complicated multiphase dynamics, especially for those with large-density ratios.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A Pair of 2D Quantum Liquids: Investigating the Phase Behavior of Indirect Excitons

Long-lived indirect excitons (IXs) exhibit a rich phase diagram, including a Bose–Einstein condensate (BEC), a Wigner crystal, and other exotic phases. Recent experiments have hinted at a “classical” liquid of IXs above the BEC transition. To uncover the nature of this phase, we use a broad range of theoretical tools and find no evidence of a driving force toward classical condensation. Instead, we attribute the condensed phase to a quantum electron–hole liquid (EHL), first proposed by Keldysh for direct excitons. Taking into account the association of free carriers into bound excitons, we study the phase equilibrium between a gas of excitons, a gas of free carriers, and an EHL for a wide range of electron–hole separations, temperatures, densities, and mass ratios. Our results agree reasonably well with recent measurements of GaAs/AlGaAs coupled quantum wells.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Dependence of Enstrophy Transport and Mixed Mass on Dimensionality and Initial Conditions in the Richtmyer–Meshkov Instability Induced Flows

Here we present a comparative study of the enstrophy budget and mixed mass between two- and three-dimensional flows induced by Richtmyer–Meshkov instability (RMI). Specifically, the individual contributions to the enstrophy budget due to the production from baroclinicity and from vortex stretching (which vanishes in two-dimensional (2D) flow) are delineated. This is enabled by a set of two- and three-dimensional computations at Atwood 0.5 having both narrow- and broad-band perturbations. A further three-dimensional (3D) computation is conducted at Atwood 0.9 using an identical narrowband perturbation to the Atwood 0.5 case to examine the sensitivity to density ratio. The mixed mass is also considered with the goal to obtain insight on how faithfully a simplified calculation performed in two dimensions can capture the mixed mass for an inertial confinement fusion (ICF) or other practical application. It is shown that the late time power law decay of variable density enstrophy is substantially different in two and three dimensions for the narrowband initial perturbation. The baroclinic production term is negligible in three dimensions (aside from the initial shock interaction), as vortex stretching is larger by two orders of magnitude. The lack of vortex stretching considerably reduces the decay rate in both narrowband and broadband perturbations in two dimensions. In terms of mixed mass, the lack of vortex stretching reduces the mixed mass in two dimensions compared to three in all cases. In the broadband cases, the spectral bandwidth in the 2D case is wider; hence, there is a longer time period of sustained linear growth which reduces the normalized mixed mass further.

42 ENGINEERING↗
DOE OSTI2022

Spectral evidence of solar neighborhood analogs in CALIFA galaxies

We introduce a novel nonparametric method to find solar neighborhood analogs (SNAs) in extragalactic integral field spectroscopic surveys. The main ansatz is that the physical properties of the solar neighborhood (SN) should be encoded in its optical stellar spectrum. We assume that our best estimate of such a spectrum is the one extracted from the analysis performed by the Code for Stellar properties Heuristic Assignment (CoSHA) from the MaStar stellar library. It follows that finding SNAs in other galaxies consist in matching, in a χ 2 sense, the SN reference spectrum across the optical extent of the observed galaxies. We applied this procedure to a selection of CALIFA galaxies, by requiring a close to face-on projection, relative isolation, and non-active galactic nucleus. We explore how the local and global properties of the SNAs (stellar age, metallicity, dust extinction, mass-to-light ratio, stellar surface mass density, star-formation density, and galactocentric distance) and their corresponding host galaxies (morphological type, total stellar mass, star-formation rate, and effective radius) compared with those of the SN and the Milky Way (MW). We find that SNAs are located preferentially in S(B)a–S(B)c galaxies, in a ring-like structure, which radii seem to scale with the galaxy size. Despite the known sources of systematics and errors, most properties present a considerable agreement with the literature on the SN. We conclude that the solar neighborhood is relatively common in our sample of SNAs. Our results warrant a systematic exploration of correlations among the physical properties of the SNAs and their host galaxies. We reckon that our method should inform current models of the galactic habitable zone in our MW and other galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Nested solitons in two-field fuzzy dark matter

ABSTRACT Dark matter as scalar particles consisting of multiple species is well motivated in string theory where axion fields are ubiquitous. A two-field fuzzy dark matter (FDM) model features two species of ultralight axion particles with different masses, m1 ≠ m2, which is extended from the standard one-field model with $m_a \sim 10^{-22} \, {\rm eV}$. Here we perform numerical simulations to explore the properties of two-field FDM haloes. We find that the central soliton has a nested structure when m2 ≫ m1, which is distinguishable from the generic flat-core soliton in one-field haloes. However, the formation of this nested soliton is subject to many factors, including the density fraction and mass ratio of the two fields. Finally, we study non-linear structure formation in two-field cosmological simulations with self-consistent initial conditions and find that the small-scale structure in two-field cosmology is also distinct from the one-field model in terms of DM halo counts and soliton formation time.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The SAMI Galaxy Survey: stellar population and structural trends across the Fundamental Plane

ABSTRACT We study the Fundamental Plane (FP) for a volume- and luminosity-limited sample of 560 early-type galaxies from the SAMI survey. Using r-band sizes and luminosities from new multi-Gaussian expansion photometric measurements, and treating luminosity as the dependent variable, the FP has coefficients a = 1.294 ± 0.039, b = 0.912 ± 0.025, and zero-point c = 7.067 ± 0.078. We leverage the high signal-to-noise ratio of SAMI integral field spectroscopy, to determine how structural and stellar population observables affect the scatter about the FP. The FP residuals correlate most strongly (8σ significance) with luminosity-weighted simple stellar population (SSP) age. In contrast, the structural observables surface mass density, rotation-to-dispersion ratio, Sérsic index, and projected shape all show little or no significant correlation. We connect the FP residuals to the empirical relation between age (or stellar mass-to-light ratio Υ⋆ ) and surface mass density, the best predictor of SSP age amongst parameters based on FP observables. We show that the FP residuals (anti)correlate with the residuals of the relation between surface density and Υ⋆ . This correlation implies that part of the FP scatter is due to the broad age and Υ⋆ distribution at any given surface mass density. Using virial mass and Υ⋆, we construct a simulated FP and compare it to the observed FP. We find that, while the empirical relations between observed stellar population relations and FP observables are responsible for most (75 per cent) of the FP scatter, on their own they do not explain the observed tilt of the FP away from the virial plane.

D’Eugenio, Francesco↗
DOE OSTI2022

A fully Lagrangian, non-parametric bias model for dark matter halos

Here, we present a non-parametric Lagrangian biasing model and fit the ratio of the halo and mass densities at the field level using the mass-weighted halo field in the AbacusSummit simulations at z=0.5. Unlike the perturbative halo bias model that has been widely used in interpreting the observed large-scale structure traced by galaxies, we find a non-negative halo-to-mass ratio that increases monotonically with the linear overdensity δ 1 in the initial Lagrangian space. The bias expansion, however, does not guarantee non-negativity of the halo counts, and may lead to rising halo number counts at negative overdensities. The shape of the halo-to-mass ratio is unlikely to be described by a polynomial function of δ 1 and other quantities. Especially for massive halos with 6×10 12 h -1 M ⊙ , the halo-to-mass ratio starts soaring up at δ 1 >0, substantially different from the predictions of the bias expansion. We show that for the halo masses we consider (M>3×10 11 h -1 M ⊙ ) a non-parametric halo-to-mass ratio as a function of δ 1 and its local derivative ∇^2δ 1 can recover the halo power spectra to sub-percent level accuracy for wavenumbers k=0.01-0.1 h Mpc -1 given a proper smoothing scale to filter the initial density field, even though we do not fit the power spectrum directly. However, there is mild dependence of the recovery of the halo power spectrum on the smoothing scale and other input parameters. At k<0.01 h Mpc -1 and for massive halos with M>6x10 12 h -1 M ⊙ , our non-parametric model leads to a few percent overestimation of the halo power spectrum, indicating the need for larger or multiple smoothing scales. The halo-to-mass ratios obtained qualitatively agree with intuitions from extended Press-Schechter theory. We compare our framework to the bias expansion and discuss possible extensions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Satellite Reentry Predictions During Sudden Stratospheric Warmings

Prediction of uncontrolled satellite reentry time and trajectory are essential to avoid damages to human being and properties over land. Reentry time and trajectory are largely controlled by aerodynamic drag and mass neutral density below 200 km. In this paper, using a numerical model of the whole atmosphere, WACCM‐X, and a Precision Orbital Determination software, we demonstrate the major sudden stratospheric warming (SSW) in January 2009, during the deep solar minimum of 2009, may cause 0.2–0.3 day of variation in reentry time. The SSW can alter the mass density and composition mixing ratio in the lower thermosphere below 200 km through changed wave forcing and general circulation in the lower thermosphere. This is comparable to the effect of a 2‐SFU (Solar Flux Unit) solar radiance variation on the neutral density. 0.1 day of reentry time uncertainty corresponds to a global shift in impact location, that is, thousands of kilometers distance. For comparison, a major SSW during solar maximum contributes little to reentry time variation, which is predominantly controlled by the solar and geomagnetic forcing. We suggest a few proxies such as O/N 2 in the lower thermosphere or mean meridional wind to account for the lower atmosphere impacts in an operational model. This needs to be investigated in rigorous future work.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Development of scalable design optimization parameters for bi-component protective systems

Bi-component protection-systems need to withstand a variety of extreme loading conditions such as high-velocity projectile impact to shock and blast mitigation. The large number of parameters at each stage of bi-component protection system design process results in a near-infinite number of possible variables for optimization. The aim of this project is to identify the scaling laws linking material and geometric scales of the problem to the final impact performance, further including the effects of wave propagation/interaction and inter-component bonding. Specifically, the impact dynamics for different optimal projectile/target density ratios and relative component mass fractions are quantified for bonded ceramic/titanium alloy systems. The impact behavior such as wave mechanics and damage propagation within the target system is characterized and quantified via digital image correlation or via x-ray phase-contrast imaging.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The role of supernova convection for the lower mass gap and the isolated binary formation of gravitational wave sources

Understanding astrophysical phenomena involving compact objects requires an insight about the engine behind core-collapse supernovae (SNe) and the fate of the stellar collapse of massive stars. In particular, this insight is crucial in developing an understanding of the origin and formation channels of the growing populations of detected black hole-black hole, black holeneutron star and neutron star-neutron star mergers. The timescale of convection growth may have a large effect on the strength of SN explosion and therefore also on the mass distribution of stellar remnants. We adopt new formulas for the relation between the pre-SN star properties and their remnants (Fryer et al. 2022) and check how they impact population of double compact object (DCO) mergers formed via isolated binary evolution. The new formulas give one ability to test wide spectrum of assumptions on the convection growth time. In particular, different variants allow for a smooth transition between having a deep mass gap and a remnant mass distribution filled by massive neutron stars and low mass black holes. We present distribution of masses, mass ratios and the local merger rate densities of DCO for different variants of new formulas and test them together with different approaches to other highly uncertain processes. We find that mass distribution of DCO mergers is sensitive to adopted assumption on SN convection growth timescale up to $m_1$ + $m_2$ ≲ 35$M_⊙$. Between two extreme tested variants the probability of compact object formation within the mass gap may differ up to ~ 2 orders of magnitude.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Acceleration of uranium beam to record power of 10.4 kW and observation of new isotopes at Facility for Rare Isotope Beams

The Facility for Rare Isotope Beams (FRIB) is a major nuclear physics facility for research with fast, stopped, and reaccelerated beams that was successfully commissioned in May 2022. A key capability of FRIB is the production of an acceleration of the uranium beam, but this capability requires the facility to work at the design limits of the lowest charge-to-mass ratio and the highest power density on the beam intercepting devices. This paper presents techniques for overcoming the significant challenges in accelerating the uranium beam, culminating in the demonstration of 10.4 kW on target, and the discovery of three new isotopes. The high-power uranium beam enabled us to produce and identify G 88 a , A 93 s , and S 96 e , within the first 24 h of operation. The successful uranium operation at FRIB sets a new record for accelerated uranium beam power above 10 kW and opens a new avenue of research with rare isotopes. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

The role of supernova convection for the lower mass gap in the isolated binary formation of gravitational wave sources

ABSTRACT Understanding the astrophysical phenomena involving compact objects requires an insight about the engine behind the core-collapse supernovae (SNe) and the fate of the stellar collapse of massive stars. In particular, this insight is crucial in developing an understanding of the origin and formation channels of the growing populations of the detected black hole–black hole, black hole–neutron star, and neutron star–neutron star mergers. The time-scale of convection growth may have a large effect on the strength of SN explosion and therefore also on the mass distribution of stellar remnants. We adopt new formulas for the relation between the pre-SN star properties and their remnants and check how they impact the population of double compact object (DCO) mergers formed via the isolated binary evolution. The new formulas give one the ability to test a wide spectrum of assumptions on the convection growth time. In particular, the different variants allow for a smooth transition between having a deep mass gap and a remnant mass distribution filled by massive neutron stars and low-mass black holes. We present the distribution of masses, mass ratios, and the local merger rate densities of DCO for the different variants of new formulas and test them together with different approaches to other highly uncertain processes. We find that the mass distribution of DCO mergers is sensitive to the adopted assumption on the SN convection growth time-scale up to $m_1+m_2 \lesssim 35 \, \mathrm{M}_{\odot }$. Between the two extreme tested variants the probability of compact object formation within the mass gap may differ by up to approximately two orders of magnitude.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The differences between mass- and light-derived structural parameters over time for MaNGA elliptical galaxies

ABSTRACT We apply stellar population synthesis analysis to obtain spatially resolved archaeological inferences for a large sample of ‘red and dead’ Elliptical galaxies (Classical Ellipticals, CLEs) from the MaNGA/SDSS-IV DR15 survey. From their 2D stellar light and mass maps, we explore the differences between the radial mass and light distributions in the rest-frame bands g, r, and i as functions of look-back time, tlb, or redshift, z. We characterize these differences through the ratios between the following mass- and light-derived global properties: sizes, concentrations, and effective surface densities. We find that the mass-to-light ratios of these properties change with tlb, more the more massive the galaxies are. The CLE galaxy archaeological progenitors are, on average, less compact, concentrated, and dense in light than in mass as z decreases. However, at later times, when also the evolution of the progenitors becomes passive at all radii, there is an upturn in these trends and the differences between mass and light in compactness/concentration decrease towards z ∼ 0. The trends in the ratios of mass-to-light sizes agree qualitatively with results from direct observations in galaxy surveys at different redshifts. We discuss the caveats and interpretations of our results, and speculate that the strong structural evolution found in some previous studies for early-type galaxies could be explained partially by photometric changes rather than by intrinsic structural changes.

Ibarra-Medel, H. (ORCID:0000000297906313)↗
DOE OSTI2025

An extension of the localized artificial diffusivity method for immiscible and high density ratio flows

The localized artificial diffusivity (LAD) method is widely regarded as the preferred multi-material regularization scheme for the compact finite difference method, because it is conservative, easy to implement, and generally robust for a wide range of multi-material problems. However, traditional LAD methods face significant challenges when applied to flows with large density ratios and when maintaining thermodynamic equilibrium across material interfaces. These limitations arise from the formulation of the artificial diffusivity flux and the reliance on enthalpy diffusion for interface regularization. Additionally, traditional LAD methods struggle to ensure stability under large density ratio conditions, fail to maintain a finite interface thickness, and are therefore unsuitable for modeling immiscible interfaces. Here, in this work, we discuss the origins of these issues in traditional LAD methods and propose modifications which enable the simulation of large density ratio and immiscible flows. The proposed method targets the artificial diffusion fluxes at gradients and ringing in the volume fraction, rather than the mass fraction in traditional methods, to consistently regularize large density ratio interfaces. Furthermore, the proposed method introduces an artificial bulk density diffusion term to enforce equilibrium conditions across interfaces. To address the challenge of modeling immiscible flows, a conservative diffuse interface term is incorporated into the formulation to ensure a finite interface thickness. Specific consideration is taken in the design of the method to ensure that these crucial properties are maintained for N -material flows. The effectiveness of the proposed method is demonstrated through a series of canonical test cases, and its accuracy is validated by comparison with experimental data on micro-bubble collapse in water. These results highlight the method’s robustness and its ability to overcome the limitations of traditional LAD approaches.

Artificial diffusivity↗