Search NASA⌕ Search

SEARCH · Search NASA

Results for “density”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
NASA NTRS2012

Reliability of High I/O High Density CCGA Interconnect Electronic Packages under Extreme Thermal Environment

This paper provides the experimental test results of advanced CCGA packages tested in extreme temperature thermal environments. Standard optical inspection and x-ray non-destructive inspection tools were used to assess the reliability of high density CCGA packages for deep space extreme temperature missions. Ceramic column grid array (CCGA) packages have been increasing in use based on their advantages such as high interconnect density, very good thermal and electrical performances, compatibility with standard surface-mount packaging assembly processes, and so on. CCGA packages are used in space applications such as in logic and microprocessor functions, telecommunications, payload electronics, and flight avionics. As these packages tend to have less solder joint strain relief than leaded packages or more strain relief over lead-less chip carrier packages, the reliability of CCGA packages is very important for short-term and long-term deep space missions. We have employed high density CCGA 1152 and 1272 daisy chained electronic packages in this preliminary reliability study. Each package is divided into several daisy-chained sections. The physical dimensions of CCGA1152 package is 35 mm x 35 mm with a 34 x 34 array of columns with a 1 mm pitch. The dimension of the CCGA1272 package is 37.5 mm x 37.5 mm with a 36 x 36 array with a 1 mm pitch. The columns are made up of 80% Pb/20%Sn material. CCGA interconnect electronic package printed wiring polyimide boards have been assembled and inspected using non-destructive x-ray imaging techniques. The assembled CCGA boards were subjected to extreme temperature thermal atmospheric cycling to assess their reliability for future deep space missions. The resistance of daisy-chained interconnect sections were monitored continuously during thermal cycling. This paper provides the experimental test results of advanced CCGA packages tested in extreme temperature thermal environments. Standard optical inspection and x-ray non-destructive inspection tools were used to assess the reliability of high density CCGA packages for deep space extreme temperature missions. Keywords: Extreme temperatures, High density CCGA qualification, CCGA reliability, solder joint failures, optical inspection, and x-ray inspection.

high density CCGA qualifications↗
NASA NTRS2014

Core and Wing Densities of Asymmetric Coronal Spectral Profiles: Implications for the Mass Supply of the Solar Corona

Recent solar spectroscopic observations have shown that coronal spectral lines can exhibit asymmetric profiles, with enhanced emissions at their blue wings. These asymmetries correspond to rapidly upflowing plasmas at speeds exceeding approximately equal to 50 km per sec. Here, we perform a study of the density of the rapidly upflowing material and compare it with that of the line core that corresponds to the bulk of the plasma. For this task, we use spectroscopic observations of several active regions taken by the Extreme Ultraviolet Imaging Spectrometer of the Hinode mission. The density sensitive ratio of the Fe(sub XIV) lines at 264.78 and 274.20 Angstroms is used to determine wing and core densities.We compute the ratio of the blue wing density to the core density and find that most values are of order unity. This is consistent with the predictions for coronal nanoflares if most of the observed coronal mass is supplied by chromospheric evaporation driven by the nanoflares. However, much larger blue wing-to-core density ratios are predicted if most of the coronal mass is supplied by heated material ejected with type II spicules. Our measurements do not rule out a spicule origin for the blue wing emission, but they argue against spicules being a primary source of the hot plasma in the corona. We note that only about 40% of the pixels where line blends could be safely ignored have blue wing asymmetries in both Fe(sub XIV) lines. Anticipated sub-arcsecond spatial resolution spectroscopic observations in future missions could shed more light on the origin of blue, red, and mixed asymmetries.

coronal spectral profiles↗
NASA NTRSDate not supplied

Urban Air Mobility Airspace Dynamic Density

Airspace safety must be assured for Urban Air Mobility (UAM) to become reality. Emerging operations need to safely integrate with existing and future air traffic. UAM is anticipated to evolve in stages. Early stages, characterized by low-density operations, may be accommodated by traditional air traffic management approaches. However, to enable the scale of UAM operations necessary to reach ubiquitous integration into daily life, new paradigms that utilize collaborative and automated systems are being considered. Correspondingly, traditional approaches for determining a safe number of simultaneous operations in an airspace must be adapted to suit evolving traffic management paradigms. In this paper, we examine dynamic density as an approach to determine when an airspace is excessively populated, and safety may become compromised. We review previous work on dynamic density as it relates to air traffic controller workload and propose factors that may be more fitting for UAMs. Rather than considering controller workload, we propose factors that suggest increased likelihood of loss of separation and may lead to conflicts that require tactical avoidance maneuvers. In addition to monitoring dynamic density for tactical decision making, we endeavor to predict airspace dynamic density with adequate look-ahead to allow for strategic route selection.

UAM↗
NASA NTRSDate not supplied

Urban Air Mobility Airspace Dynamic Density

Effective flight planning requires information about a variety of potential threats, such as adverse weather or airspace restrictions, and about alternatives available if unforeseen events occur. Expected traffic along the route of flight is also essential to a safe outcome so that, for example, adequate fuel/energy supply can be loaded prior to flight. A dynamic density (DD) metric is introduced for the emerging urban air mobility (UAM) concept to predict airspace congestion that may lead to loss of separation between aircraft or less efficient operations. Using inspiration from dynamic density metric research for traditional air traffic management and a two-way highway analogy, we develop a dynamic density metric for a portion of airspace (aUAM corridor) that aggregates the impact from five factors: aircraft density, density of populous clusters, mean number of aircraft in populous clusters, mean distance between aircraft, and minimum distance between aircraft. This works describes our methodology, rationale, use cases, and visualization techniques to efficiently present the DD metric to an operator for informed decision making. We also present an approach for validating the metric. However, validation remains part of future work.

UAM↗
DOE OSTI2022

Experimental Steps toward a Density Law for Chlorine-Crediting Criticality Models of Aqueous Plutonium solutions

Criticality accident prevention is an essential safety consideration for all operations involving fissionable and fissile materials. Predictive criticality calculations using advanced neutron transport codes, such as Monte Carlo N-Particle Transport (MCNP), are invaluable tools for designing and implementing operational limits. While indispensable, these tools are limited by the quality and accuracy of the inputs that the user provides to define the modeled system. Parameters such as atomic composition, shape, and material density must be accurately defined to obtain a meaningful result. In the case of material property details for a plutonium model, accurate material characterization data is sometimes sparse. Thus limitations in calculation accuracy can be reduced by improving our knowledge of material properties in the targeted fissile systems. One of the biggest remaining challenges to accurately defining fissile systems is the description of aqueous fissile solutions. Even simple properties, such as density, are not well-known for fissile solutions relevant to nuclear energy and security. Described here are initial efforts undertaken to improve criticality calculation inputs for fissile plutonium chloride solutions in water. This effort is focused on experimentally determining accurate solution characteristics for ternary plutonium chloride/hydrochloric acid/water systems, by measurement of water activity and solution density. The effect of inputting experimental densities for these solutions into MCNP criticality calculations is compared to the traditional approach of modeling an idealized (and fictitious) plutonium-water mixture. Expansion of these efforts to a working density law for aqueous plutonium chloride solutions is also discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

CHARGED POINT DEFECTS AND DEFECT COMPLEXES IN PHOSPHORENE AND STRUCTURAL PREDICTION OF TAU-11, A DENSITY FUNCTIONAL THEORY APPROACH

With the advent of powerful computers, first principle calculations have become a robust and reliable tool to complement experiments in the study of materials. First principle calculations such as density functional theory can be used to guide experiments to save time and effort when it comes to new and advanced materials discovery. Study of charge defects in materials and crystal structure prediction are such areas in which density functional theory calculations facilitate the prediction of defect formation energy, charge transition levels and elemental site occupations. Experimentally, studying charge defects requires inferring results from various techniques. Likewise, determining site occupancies for similar elements using experimental techniques is challenging and sometimes impossible. In this work, we study the charged intrinsic and extrinsic defects and defect complexes in phosphorene and determine the site occupancies of Al, Fe and Si in the τ11 phase using density functional theory. In the defect work, we calculate the defect formation energy, charge transition level, binding energies and Stokes shift for vacancy, dopant substitution and dopant-vacancy defect complexes in phosphorene. We found that vacancy defect in phosphorene becomes negatively charged in n-type doping and may passivate the dopants and reduce carrier concentration and mobility. For non-metal dopants in phosphorene, we predict that O, S and Si prefer to form dopant-vacancy complexes removing the vacancy defect states from the band gap. Mn dopant-vacancy defect complex exhibits possibility of switching between two magnetic spin states. Lastly, using density functional theory, we complement the results from neutron powder diffraction to determine the site occupancies of Al and Si in τ11 phase which was used to determine the alloying element to develop a high temperature and low density alloy.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Evaluation of Thermal Scattering Law for U 3 Si 2 Using Temperature-Dependent Phonon Density of States

The thermal scattering law for uranium silicide, a candidate material for high-density nuclear fuel, was evaluated from the temperature-dependent phonon density of states obtained using the quasi-harmonic approximation. For uranium silicide model systems of various volumes near the equilibrium volume, Hellmann-Feynman forces were calculated by performing ab initio lattice dynamics simulations based on density functional theory. Subsequently, the quasi-harmonic approximation was applied by interpolating the results of the harmonic approximation for each volume system. As a sensitivity test, the phonon density of states was calculated for the structures with temperature-dependent volumes at 298 K and 1000 K, obtained by the quasi-harmonic approximation, and the corresponding thermal neutron scattering cross section was evaluated. The temperature-dependent volume and heat capacity of uranium silicide obtained by the quasi-harmonic approximation showed reasonable agreement with experimental data while lattice thermal expansion showed deviation from experimental observations. Thermal scattering cross sections evaluated from 0 K and temperature-dependent phonon density of states exhibited differences of up to 4.1% at 1000 K. This study examines the impact of temperature on the thermal scattering law, employing ab initio lattice dynamics simulations through the quasi-harmonic approximation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Simple hydrogenic estimates for the exchange and correlation energies of atoms and atomic ions, with implications for density functional theory

Exact density functionals for the exchange and correlation energies are approximated in practical calculations for the ground-state electronic structure of a many-electron system. An important exact constraint for the construction of approximations is to recover the correct non-relativistic large-Z expansions for the corresponding energies of neutral atoms with atomic number Z and electron number N = Z, which are correct to leading order (–0.221 Z 5/3 and –0.021 Z ln Z respectively) even in the lowest-rung or local density approximation. We find that hydrogenic densities lead to E x (N,Z) ≈ –0.354 N 2/3 Z (as known before only for Z >> N >> 1) and E c ≈ –0.02 N ln N . These asymptotic estimates are most correct for atomic ions with large N and Z >> N, but we find that they are qualitatively and semi- quantitatively correct even for small N and for N ≈ Z. The large-N asymptotic behavior of the energy is pre-figured in small-N atoms and atomic ions, supporting the argument that widely-predictive approximate density functionals should be designed to recover the correct asymptotics. Here, it is shown that the exact Kohn-Sham correlation energy, when calculated from the pure ground-state wavefunction, should have no contribution proportional to Z in the Z → ∞ limit for any fixed N.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Density Dependence in Three Snake River Sockeye Salmon Nursery Lakes in Central Idaho

Abstract Snake River Sockeye Salmon Oncorhynchus nerka, listed as an endangered species in 1991, currently inhabit three nursery lakes (Redfish, Pettit, and Alturas lakes) in the Sawtooth Valley, Idaho. Conspecific kokanee (lacustrine Sockeye Salmon) are also present in the lakes. Snake River Sockeye Salmon recovery efforts, initially focused on genetic conservation, are now attempting to rebuild naturally spawning populations using hatchery supplementation. However, in Sockeye Salmon nursery lakes, density dependence is frequently observed when elevated O. nerka abundance leads to declines in zooplankton biomass, body size, and shifts in community composition. In turn, these changes lead to reductions in juvenile O. nerka growth rates, survival, and adult returns. We examined a long-term data set of O. nerka population metrics and associated zooplankton community metrics. We found evidence of density dependence within and among nursery lakes. We detected differences in zooplankton biomass, lengths of preferred zooplankton prey (Daphnia spp. and cyclopoid copepods), parr growth rates, and age-1 smolt size among the three lakes. We found negative relationships between O. nerka density and zooplankton biomass and size. We identified positive relationships between zooplankton biomass and two response variables: smolt size at migration and growth rates of hatchery parr. The relationships were generally similar among lakes. Variable outcomes were a result of differences in O. nerka density (or zooplankton biomass), controlled primarily by the relative proportion of spawning and rearing habitat in each lake. Understanding unique lake habitats, ecological interactions, and the role of density dependence is germane to management of Snake River Sockeye Salmon populations.

Griswold, Robert G.↗
DOE OSTI2024

Spatial patterns in occupancy and density of larval lampreys in freshwater habitats restored to a Stage 0 condition

Abstract We examined occupancy and density of larval lampreys ( Entosphenus tridentatus and Lampetra spp.) in two rivers in Oregon (USA) restored to a Stage 0 condition 1–5 years prior, using a multiscale occupancy model and a zero‐inflated Poisson mixture model. We sampled lampreys using backpack electrofishing in randomly distributed, paired, 1‐m 2 quadrats and recorded environmental data. Probabilities of occupancy and density were higher when water velocity was low, the substrate was noncompacted, and sediment was dominated by fines (<4 mm). At mean water depth (0.34 m) and velocity (0.09 m/s), estimated densities in occupied quadrats were 4.8 lampreys/m 2 (95%: 3.4–6.9) when the substrate was compacted, and fines were not dominant, and 21.1 lampreys/m 2 (95%: 17.7–25.3) when the substrate was noncompacted and fines were dominant. Probabilities of detecting occupancy in a 1‐m 2 quadrat sampled by backpack electrofishing were 0.76 (95%: 0.64–0.87) when captured after visual observation and 0.80 (95%: 0.71–0.88) with blind sweeps (i.e., constantly moving the net regardless of observation). The probability of capturing a single lamprey in a quadrat sampled by blind sweeps was 0.32 (95%: 0.27–0.37). Sampling in paired 1‐m 2 quadrats facilitated concurrent examination of patterns in occupancy and density while accounting for capture probability, which could aid temporal monitoring of restored habitats. To the best of our knowledge, this is the first study to document occupancy and estimate densities of larval lampreys in habitats that underwent valley floor restoration to Stage 0. We observed both lamprey genera within 5 years of restoration. Aquatic restoration that increases low‐velocity, noncompacted, fine sediment habitats could benefit lampreys.

Harris, Julianne E.↗
DOE OSTI2023

The connection between nonzero density and spontaneous symmetry breaking for interacting scalars

We consider U(1)-symmetric scalar quantum field theories at zero temperature. At nonzero charge densities, the ground state of these systems is usually assumed to be a superfluid phase, in which the global symmetry is spontaneously broken along with Lorentz boosts and time translations. We show that, in d > 2 spacetime dimensions, this expectation is always realized at one loop for arbitrary non-derivative interactions, confirming that the physically distinct phenomena of nonzero charge density and spontaneous symmetry breaking occur simultaneously in these systems. We quantify this result by deriving universal scaling relations for the symmetry breaking scale as a function of the charge density, at low and high density. Moreover, we show that the critical value of μ above which a nonzero density develops coincides with the pole mass in the unbroken, Poincaré invariant vacuum of the theory. The same conclusions hold non-perturbatively for an O(N) theory with quartic interactions in d = 3 and 4, at leading order in the 1/N expansion. We derive these results by computing analytically the zero-temperature, finite-μ one-loop effective potential, paying special attention to subtle points related to the iε terms. We check our results against the one-loop low-energy effective action for the superfluid phonons in λΦ 4 theory in d = 4 previously derived by Joyce and ourselves, which we further generalize to arbitrary potential interactions and arbitrary dimensions. As a byproduct, we find analytically the one-loop scaling dimension of the lightest charge-n operator for the λΦ 6 conformal superfluid in d = 3, at leading order in 1/n, reproducing a numerical result of Badel et al. For a λΦ 4 superfluid in d = 4, we also reproduce the Lee-Huang-Yang relation and compute relativistic corrections to it. Finally, we discuss possible extensions of our results beyond perturbation theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Scattering Observables from One- and Two-body Densities: Formalism and Application to $\pmb \gamma $ ${}^3\hbox {He}$ Scattering

We introduce the transition-density formalism, an efficient and general method for calculating the interaction of external probes with light nuclei. One- and two-body transition densities that encode the nuclear structure of the target are evaluated once and stored. They are then convoluted with an interaction kernel to produce amplitudes, and hence observables. Here, by choosing different kernels, the same densities can be used for any reaction in which a probe interacts perturbatively with the target. The method therefore exploits the factorisation between nuclear structure and interaction kernel that occurs in such processes. We study in detail the convergence in the number of partial waves for matrix elements relevant in elastic Compton scattering on 3 He. The results are fully consistent with our previous calculations in Chiral Effective Field Theory. But the new approach is markedly more computationally efficient, which facilitates the inclusion of more partial-wave channels in the calculation. We also discuss the usefulness of the transition-density method for other nuclei and reactions. Calculations of elastic Compton scattering on heavier targets like 4 He are straightforward extensions of this study, since the same interaction kernels are used. And the generality of the formalism means that our 3 He densities can be used to evaluate any 3 He elastic-scattering observable with contributions from one- and two-body operators.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Investigation of ventilation-coupled high energy density sensible thermal energy storage

Low-cost energy storage will play an important role in supporting the decarbonization of the energy sector. Here, a novel approach to thermal energy storage for buildings is proposed, in which a tank of antifreeze liquid can be used to heat outdoor air in mechanical ventilation systems. If heated to a high enough temperature, the fluid can potentially undergo temperature swings of nearly 100°C during discharge in cold climates. Because the energy density of sensible storage systems scales proportionally to the system temperature change, this concept has the potential to offer higher energy densities than other liquid-based sensible storage devices used in building applications. A one-dimensional numerical model of the storage system was developed and experimentally validated using 30 wt% (wt%) potassium acetate as an aqueous antifreeze solution. Using a commercial hot water tank, energy densities of 47.0 kW-hour per cubic meter (kWh/m 3 ) were demonstrated in a laboratory setting. Material energy densities of approximately 74.7 (kWh/m 3 ) were measured. Design improvements may boost this energy storage density even further. Because of the system's relative simplicity, ventilation-coupled sensible storage has the potential to be an easily deployable, low-cost energy storage solution for building systems.

25 ENERGY STORAGE↗
DOE OSTI2025

Packing density optimization and recyclability of multi-modal SiC powder feedstocks in binder jet additive manufacturing

Achieving high powder packing density is critical in binder jet additive manufacturing (BJAM), as it directly influences the final part density, mechanical properties, and sintering behavior. Multi-modal powder blends, which combine particles of different sizes, have been explored as a strategy to optimize packing efficiency and minimize defects. In this study, bimodal and trimodal powder blends were obtained by mixing silicon carbide powder in three different sizes. These results show that increased powder density is achievable with bimodal powder blends but is reduced in trimodal blends, and it was found that a 13 % increase in the powder tap density was achieved using a bimodal blend of powder. The powder size distribution of the bimodal blend was measured at various stages during binder jet additive manufacturing and, no measurable powder separation occurred even after eight prints. Altogether, this study shows limited advantage to trimodal powder blends but good promise for bimodal blends for increasing printed density while maintaining reusability in the binder jet process.

Aguirre, Trevor G. [Oak Ridge National Laboratory ↗
DOE OSTI2021

An efficient computational framework for charge density estimation in twisted bilayer graphene

Electronic properties such as band structure and Fermi velocity in low-angle twisted bilayer graphene (TBG) are intrinsically dependent on the atomic structure. Rigid rotation between individual graphene layers provides an approximate description of the bilayer symmetry. Upon relaxation, in-plane displacement of the atoms in low angle TBG causes a change in the symmetry through the enlargement of the AB stacking regions and the reduction in size of AA and SP stacking regions. However, the effect of this in-plane relaxation on the charge density remains unexplored, because the necessary electronic structure calculations of such large supercells of low twist angle TBG are computationally infeasible. Therefore, we develop a computationally efficient framework that enables the exploration of the charge density symmetry of the low twist angle TBG. This framework is based on the Fourier representation of the charge density which presents high intensity Bragg peaks. Here we find that with the decrease of twist angle, low intensity satellite peaks also become apparent. Our framework incorporates these satellite peaks which reveals transformation of symmetry in the charge density distribution from high to low twist angle TBG. One striking outcome is the demonstration of the electron localization in the AA region of low twist angle TBG. Our framework helps to explain the effect of the atomistic relaxation on the charge density distribution and thus, it provides information about exotic electronic properties of low twist angle TBG at a low computational expense.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Herbicide, fertilization, and planting density effects on intensively managed loblolly pine early stand development

Production forestry in the southeast US has been partially transitioned to intensively managed short rotations (~10 years), in which multiple silvicultural interventions are performed during forest development. Understanding the responses to silvicultural practices and continued refinement of site-specific recommendations is critical to sustainably maximize forest production. We evaluated the effects of silvicultural practices (herbicide, fertilization, and planting density) on growth, stand homogeneity, and above- and belowground biomass accumulation and partitioning of loblolly pine (Pinus taeda) throughout early stand development (age 5 years) in the southeastern US. Five treatments with eight replications each were tested: no herbicide and no fertilization (C); herbicide only (H); herbicide and half-reduced fertilization rate (R); herbicide and full fertilization (F); and increased stand density (60+ %; 1346 vs. 2152 trees per hectare) with herbicide and full fertilization rate (D). Allometric equations generated from destructive harvests were applied to annual diameter measurements to estimate plot-level biomass and allocation. Herbicide was crucial to promote stand uniformity and increase yield (~600+ % stem biomass compared with C at age 5). Aboveground biomass was similar in R and F treatments, which was ~25% higher than in H at age 5. Increasing planting density along with multiple herbicide and fertilizer applications yielded higher biomass without compromising individual tree size (diameter and height). There was little effect of silviculture practices on allocation patterns. Our results parallel what was found for fertilization with herbicide from a number of loblolly stands under similar conditions and indicate a ~28% volume gain with fertilization during early stand development. Similarly, our results were consistent with other studies implementing similar differences in planting density and suggest a ~26% volume gain through early stand development with an initial 60% increase in planting density. Furthermore, our study helps to understand complex relationships between production and silvicultural practices during early stand development and demonstrates that silvicultural prescriptions can be optimized to increase sustainability of production forestry

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Probability Density Function Control for Stochastic Nonlinear Systems using Monte Carlo Simulation

This paper presents a fast implementation framework of output probability density function (PDF) control for a class of stochastic nonlinear systems which are subjected to non-Gaussian noises. The statistical properties of the system outputs can be adjusted by shaping the dynamic output probability density function to track the reference stochastic distribution. However, the dynamic probability density function evolution is very difficult to obtain analytically even if the system model and the stochastic distributions of the noises are known. Motivated by Monte Carlo simulation, the dynamic probability density function can be estimated by sampling data which forms the contribution of this paper. In particular, the sampling points are generated following the stochastic distribution of the noise for each instant. These points go through the system would generate the histogram for system outputs, then the dynamic model can be established based on the dynamic histogram which reflects the randomness and the nonlinear dynamics of the investigated system. Based on the established model, the output probability density function tracking can be achieved and the simulation results and discussions show the effectiveness and benefits of the presented framework.

Zhang, Qichun↗
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↗