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 1,045 records · Page 58

Numerical study of variable density turbulence interaction with a normal shock wave

Accurate numerical simulations of shock–turbulence interaction (STI) are conducted with a hybrid monotonicity-preserving–compact-finite-difference scheme for a detailed study of STI in variable density flows. Theoretical and numerical assessments of data confirm that all turbulence scales as well as the STI are well captured by the computational method. Linear interaction approximation (LIA) convergence tests conducted with the shock-capturing simulations exhibit a similar trend of converging to LIA predictions to shock-resolving direct numerical simulations (DNS). The effects of density variations on STI are studied by comparing the results corresponding to an upstream multi-fluid mixture with the single-fluid case. Here the results show that for the current parameter ranges, the turbulence amplification by the normal shock wave is much higher and the reduction in turbulence length scales is more significant when strong density variations exist. Turbulent mixing enhancement by the shock is also increased and stronger mixing asymmetry in the postshock region is observed when there is significant density variation. The turbulence structure is strongly modified by the shock wave, with a differential distribution of turbulent statistics in regions having different densities. The dominant mechanisms behind the variable density STI are identified by analysing the transport equations for the Reynolds stresses, vorticity, normalized mass flux and density specific volume covariance.

42 ENGINEERING↗

Density Functionals for Hydrogen Storage: Defining the H2Bind275 Test Set with Ab Initio Benchmarks and Assessment of 55 Functionals

Efficient and high-capacity storage materials are indispensable for a hydrogen-based economy. In silico tools can accelerate the process of discovery of new adsorbent materials with optimal hydrogen adsorption enthalpies. Density functional theory is well-poised to become a very useful tool for enabling high-throughput screening of potential materials. In this work, we have identified density functional approximations that provide good performance for hydrogen binding applications following a two-pronged approach. First, we have compiled a data set (H2Bind275) that comprehensively represents the hydrogen binding problem capturing the chemical and mechanistic diversity in the binding sites encountered in hydrogen storage materials. We have also computed reference interaction energies for this data set using coupled-cluster theory. Second, we have assessed the performance of 55 density functional approximations for predicting H 2 interaction energies and have identified two hybrid density functionals (ωB97X-V and ωB97M-V), two double hybrid density functionals (DSD-PBEPBE-D3(BJ) and PBE0-DH), and one semilocal density functional (B97M-V) as the best performing ones. We have recommended the addition of empirical dispersion corrections to systematically underbinding density functionals such as revPBE, BLYP, and B3LYP for improvements in performance at negligible additional cost. We have also recommended the usage of the def2-TZVPP basis set as it represents a good compromise between accuracy and cost, limiting the finite basis set errors to less than 1 kJ/mol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

RDX Compression, α→ γ Phase Transition, and Shock Hugoniot Calculations from Density-Functional-Theory-Based Molecular Dynamics Simulations

Prediction of the density and lattice compression properties of the α and γ phases of the hexahydro-1,3,5-trinitro-1,3,5-s-triazine (RDX) crystal and of the low-pressure α → γ phase transition upon pressure increase are general tests used to assess the accuracy of density-functional-theory- (DFT-) based computational methods and to identify the essential parameters that govern the behavior of this high-energy-density material under extreme conditions. The majority of previous DFT studies have analyzed such issues under static optimization conditions by neglecting the corresponding temperature effects. In this paper, we extend previous investigations and analyze the performance of dispersion-corrected density functional theory to predict the compression of RDX in the pressure range of 0–9 GPa and the corresponding α → γ phase transition under realistic temperature and pressure conditions. We demonstrate that, by using static dispersion-corrected density functional theory calculations, direct interconversion between the α and γ phases upon compression is not observed. This limitation can be addressed by using isobaric–isothermal molecular dynamic simulations in conjunction with DFT-D2-calculated potentials, an approach that is shown to provide an accurate description of both the crystallographic RDX lattice parameters and the dynamical effects associated with the α→ γ phase transformation. An even more comprehensive and demanding analysis was done by predicting the corresponding shock Hugoniot curve of RDX in the pressure range of 0–9 GPa. It was found that the theoretical results reproduce reasonably well the available experimental Hugoniot shock data for both the α and γ phases. Finally, the results obtained demonstrate that a satisfactory prediction of the shock properties in high-energy-density materials undergoing crystallographic and configurational transformations is possible through the combined use of molecular dynamics simulations in the isobaric–isothermal ensemble with dispersion-corrected density functional theory methods.

74 ATOMIC AND MOLECULAR PHYSICS↗

Multireference Density Functional Theory with Generalized Auxiliary Systems for Ground and Excited States

To describe static correlation, we create a new approach to density functional theory (DFT), which uses a generalized auxiliary system that is of a different symmetry, such as particle number or spin, from that of the physical system. The total energy of the physical system consists of two parts: the energy of the auxiliary system, which is determined with a chosen density functional approximation (DFA), and the excitation energy from an approximate linear response theory that restores the symmetry to that of the physical system, thus rigorously leading to a multideterminant description of the physical system. The electron density of the physical system is different from that of the auxiliary system and is uniquely determined from the functional derivative of the total energy with respect to the external potential. Our energy functional is thus an implicit functional of the physical system density, but an explicit functional of the auxiliary system density. We demonstrate that the total energy minimum and stationary states, describing the ground and excited states of the physical system, can be obtained by a self-consistent optimization with respect to the explicit variable, the generalized Kohn–Sham noninteracting density matrix. We have developed the generalized optimized effective potential method for the self-consistent optimization. Among options of the auxiliary system and the associated linear response theory, reformulated versions of the particle–particle random phase approximation (pp-RPA) and the spin-flip time-dependent density functional theory (SF-TDDFT) are selected for illustration of principle. Numerical findings show that our multireference DFT successfully describes static correlation in bond dissociation and double bond rotation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Substrate Engineered Interconnected Graphene Electrodes with Ultrahigh Energy and Power Densities for Energy Storage Applications

Supercapacitors combine the advantages of electrochemical storage technologies like high energy density batteries and high power density capacitors. At 5-10 W h Kg -1 , the energy densities of current supercapacitors are still significantly lower than energy densities of lead acid (20-35 W h Kg -1 ), Ni-metal hydride (40-100 W h Kg -1 ) and Li-ion (120-170 W h Kg -1 ) batteries. Recently, graphene based supercapacitors have shown energy density of 40-80 W h kg -1 . However, their performance is mainly limited due to the reversible agglomeration and restacking of individual graphene layers caused by π-π interactions. The restacking of graphene layers leads to significant decrease of ion-accessible surface area and the low capacitance of graphene based supercapacitors. Here, we introduce a microstructure substrate based method to produce fully delaminated and stable interconnected graphene structure using flash reduction of GO in a few seconds. With this structure, here we achieve the highest amount of volumetric capacitance obtained so far by any type of pure carbon based material. The affordable and scalable production method is capable of producing electrodes with energy density of 0.37 Wh cm -3 and power density of 416.6 W cm -3 . This electrode maintained more than 91 % of its initial capacitance after 5000 cycles. Moreover, combining with ionic liquid, this solvent free graphene electrode material is highly promising for on-chip electronics, micro supercapacitors as well as high power applications.

25 ENERGY STORAGE↗

Forest Canopy Density Effects on Snowpack Across the Climate Gradients of the Western United States Mountain Ranges

Limitations inherent to field studies such as limited spatial and temporal coverage has made it very challenging, if not impossible, to explicitly disentangle the effect of canopy density from climate on snowpack dynamics. Hence, our understanding of how variations in canopy density can affect snow processes under diverse climate conditions has been limited. To address this knowledge gap, this study uses a physics-based modeling approach to represent most key processes that drive snow accumulation and melt in the open and forest with varied forest density, across the climate gradients of the Western U.S. mountain ranges, as represented by 228 Snow Telemetry (SNOTEL) locations. Simulations suggest that, under most winter climate conditions, canopy density exerts a strong influence on under-canopy snowpack dynamics and their controlling processes. Changing canopy density can shift the direction of the modeled differential peak snow accumulation and snow disappearance date (SDD) between the open and forest. For example, while simulated snowpack lasts longer in the open than the high-density forest at most wet/warm and dry/cold locations by up to 55 days, about 90% of all locations show a longer snow duration in the low-density forest by up to 27 days as compared to the open. Interannual variability in climate can also change the magnitude and direction of canopy impact on snow processes. By comparing model simulations to published field observations, and we also identify conditions that require additional field observations to improve process-scale understanding of forest-snow interaction.

54 ENVIRONMENTAL SCIENCES↗

Inhomogeneous high temperature melting and decoupling of charge density waves in spin-triplet superconductor UTe 2

Charge, spin and Cooper-pair density waves have now been widely detected in exotic superconductors. Understanding how these density waves emerge — and become suppressed by external parameters — is a key research direction in condensed matter physics. Here we study the temperature and magnetic-field evolution of charge density waves in the rare spin-triplet superconductor candidate UTe 2 using scanning tunneling microscopy/spectroscopy. We reveal that charge modulations composed of three different wave vectors gradually weaken in a spatially inhomogeneous manner, while persisting to surprisingly high temperatures of 10–12 K. We also reveal an unexpected decoupling of the three-component charge density wave state. Our observations match closely to the temperature scale potentially related to short-range magnetic correlations, providing a possible connection between density waves observed by surface probes and intrinsic bulk features. Importantly, charge density wave modulations become suppressed with magnetic field both below and above superconducting T c in a comparable manner. Our work points towards an intimate connection between hidden magnetic correlations and the origin of the unusual charge density waves in UTe 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Charge density wave in a band insulator

Charge-density-wave implies a periodic modulation of the charge density. Typically observed in metallic systems, charge-density-wave arise from Fermi surface instabilities, resulting in the total or partial gapping of the Fermi surface. Here, we present experimental evidence for a charge-density wave state emerging in a band insulator which has no Fermi surface. The bulk and surface of our material platform, 𝛂-Bi 4 Br 4 , is gapped over the entire Brillouin zone. Through topographic and spectroscopic imaging at low temperatures, we unveil an unexpected unidirectional charge modulation in 𝛂-Bi 4 Br 4 , breaking the lattice translation symmetry. The charge-density-wave develops at temperatures below 40 K and adds an energy gap atop the existing insulating gap of 𝛂-Bi 4 Br 4 . Furthermore, our transport measurements reveal nonlinear electrical conduction, a phenomenon conventionally associated with the sliding or phason mode of incommensurate charge-density-waves. These unusual observations represent a different type of charge-density-wave

Phase transitions and critical phenomena↗

A representation-independent electronic charge density database for crystalline materials

Abstract In addition to being the core quantity in density-functional theory, the charge density can be used in many tertiary analyses in materials sciences from bonding to assigning charge to specific atoms. The charge density is data-rich since it contains information about all the electrons in the system. With the increasing prevalence of machine-learning tools in materials sciences, a data-rich object like the charge density can be utilized in a wide range of applications. The database presented here provides a modern and user-friendly interface for a large and continuously updated collection of charge densities as part of the Materials Project. In addition to the charge density data, we provide the theory and code for changing the representation of the charge density which should enable more advanced machine-learning studies for the broader community.

36 MATERIALS SCIENCE↗

Computational studies on scattering of radio frequency waves by density filaments in fusion plasmas

In modern magnetic fusion devices, such as tokamaks and stellarators, radio frequency (RF) waves are commonly used for plasma heating and current profile control, as well as for certain diagnostics. The frequencies of the RF waves range from ion cyclotron frequency to the electron cyclotron frequency. The RF waves are launched from structures, like waveguides and current straps, placed near the wall in a very low density, tenuous plasma region of a fusion device. The RF electromagnetic fields have to propagate through this scrape-off layer before coupling power to the core of the plasma. The scrape-off layer is characterized by turbulent plasmas fluctuations and by blobs and filaments. The variations in the edge density due to these fluctuations and filaments can affect the propagation characteristics of the RF waves—changes in density leading to regions with differing plasma permittivity. Analytical full-wave theories have shown that scattering by blobs and filaments can alter the RF power flow into the core of the plasma in a variety of ways, such as through reflection, refraction, diffraction, and shadowing [see, for example, Ram and Hizanidis, Phys. Plasmas 23, 022504 (2016), and references therein]. There are changes in the wave vectors and the distribution of power–scattering leading to coupling of the incident RF wave to other plasma waves, side-scattering, surface waves, and fragmentation of the Poynting flux in the direction towards the core. However, these theoretical models are somewhat idealized. In particular, it is assumed that there is step-function discontinuity in the density between the plasma inside the filament and the background plasma. In this paper, results from numerical simulations of RF scattering by filaments using a commercial full-wave code are described. The filaments are taken to be cylindrical with the axis of the cylinder aligned along the direction of the ambient magnetic field. The plasma inside and outside the filament is assumed to be cold. There are three primary objectives of these studies. The first objective is to validate the numerical simulations by comparing with the analytical results for the same plasma description—a step-function discontinuity in density. A detailed comparison of the Poynting flux shows that numerical simulations lead to the same results as those from the theoretical model. The second objective is to extend the simulations to take into account a smooth transition in density from the background plasma to the interior of the filament. The ensuing comparison shows that the deviations from the results of the theoretical model are quite small. The third objective is to consider the scattering process for situations well beyond a reasonable theoretical analysis. This includes scattering off multiple filaments with different densities and sizes. In conclusion, simulations for these complex arrangements of filaments show that, in spite of the obvious limitations, the essential physics of RF scattering is captured by the analytical theory for a single filament

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Mitigating laser-imprint effects in direct-drive inertial confinement fusion implosions with an above-critical-density foam layer

Low-density foams of low-/mid-Z materials have been previously proposed to mitigate laser imprint for direct-drive inertial confinement fusion (ICF). For foam densities above the critical density of the drive laser, the mechanism of laser-imprint mitigation relies on the reduced growth rate of Rayleigh–Taylor instability because of the increased ablation velocity and density scale length at the ablation surface. Experimental demonstration of this concept has been limited so far to planar-target geometry. The impact of foams on spherical implosions has not yet been explored in experiments. To examine the viability of using an above-critical-density foam layer to mitigate laser-imprint effects in direct-drive ICF implosions on OMEGA, we have performed a series of 2-D DRACO simulations with state-of-the-art physics models, including nonlocal thermal transport, cross-beam energy transfer, and first-principles equation-of-state tables. The simulation results indicate that a 40-μm-thick CH or SiO 2 foam layer with a density of p = 40 mg/cm 3 added to a D2-filled polystyrene (CH) capsule can significantly improve the moderate-adiabat (α ≈ 3) implosion performance. In comparison with the standard CH target implosion, an increase of neutron yield by a factor of 4 to 8 and the recovery of 1-D compression pR are predicted by DRACO simulations for a foam-target surface roughness of σ rms ≤ 0.5 μm. These encouraging results could readily facilitate experimental demonstrations of laser-imprint mitigation with an above-critical-density foam layer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Overlapped embedded fragment stochastic density functional theory for covalently-bonded materials

The stochastic density functional theory (DFT) is a valuable linear-scaling approach to Kohn-Sham DFT that does not rely on the sparsity of the density matrix. Linear (and often sub-linear) scaling is achieved by introducing a controlled statistical error in the density, energy, and forces. The statistical error (noise) is proportional to the inverse square root of the number of stochastic orbitals and thus decreases slowly; however, by dividing the system into fragments that are embedded stochastically, the statistical error can be reduced significantly. This has been shown to provide remarkable results for non-covalently-bonded systems; however, the application to covalently bonded systems had limited success, particularly for delocalized electrons. Here, we show that the statistical error in the density correlates with both the density and the density matrix of the system and propose a new fragmentation scheme that elegantly interpolates between overlapped fragments. In conclusion, we assess the performance of the approach for bulk silicon of varying supercell sizes (up to N e = 16 384 electrons) and show that overlapped fragments reduce significantly the statistical noise even for systems with a delocalized density matrix.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Carbon ablator areal density at fusion burn: Observations and trends at the National Ignition Facility

For inertial confinement fusion experiments, the pusher is composed of a high-density deuterium tritium cyrogenic fuel layer and an ablator, often made of carbon. In an ideal, no-mix implosion, increasing the areal density of the pusher transfers more pressure to the hot spot and increases the hot spot confinement time. There has been a lack of knowledge about the final compressed state of the ablator for implosions at the National Ignition Facility. 14 MeV fusion neutrons inelastically scattering on the remaining carbon ablator excites a nuclear metastable state that emits a prompt 4.4 MeV gamma ray. The gamma reaction history diagnostic data, when reduced by a new data analysis technique, can isolate and measure the carbon gamma rays, which are proportional to the areal density of the ablator during fusion burn. The trends over many National Ignition Facility campaigns show that the ablator areal density is weakly sensitive to the maximum shell velocity, the cold fuel radius, the ablator mass remaining, or the laser picket intensity. Controlled parameter scans reveal that, for specific campaigns, ablator compression has a strong dependence on laser coast time, high Z dopants, and the laser drive foot duration. Using a model of the compressed ablator density profile reveals that the greatest variation of the ablator areal density comes from its thickness, with highly compressed, thin layers having high areal density values. The compression and thickness of the ablator are other metrics that designers should understand to differentiate the types of capsule degradation and maximize the inertial confinement fusion performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Forward and backward stimulated Raman scattering in multi-speckled beams: Density dependence and effects on cross-beam energy transfer

Particle-in-cell simulations of cross-beam energy transfer (CBET) between multi-speckled beams reveal a potentially prominent energy loss from forward stimulated Raman scattering (FSRS) in the amplified seed beam [Yin et al., Phys. Plasmas 26, 082708 (2019)]. In this work, an intensity scan is used to examine this energy loss in the CBET-amplified seed beam at intensities of a few times 1014 W/cm2. Electron tracers demonstrate both the trapping in the electron plasma waves that is sustained over multiple speckles and the subsequent acceleration to velocities well beyond the phase velocity. Motivated by expanding our understanding of this energy loss mechanism, we isolate the FSRS density dependence by using a 2×1014 W/cm2 pulse to scan densities between 0.06 ncr and 0.12 ncr; a spectral analysis characterizes the energy losses from scattering with the changing densities. At saturation, both greater beam attenuation is measured at the laser frequency and larger scattering FSRS angles are observed as density increases, along with a rapid transition to primarily backscatter SRS between 0.10 ncr and 0.12 ncr, earlier than predicted by linear theory. Density gradients are shown to effectively mitigate oblique FSRS and beam energy losses in the two lower density cases, whereas 0° forward scattering is primarily suppressed in the higher density cases. Finally returning to the two-beam CBET system of interest, the interplay and mitigation of the secondary instabilities at larger domain and beam spatial scales are discussed, here also examining backward stimulated Brillouin scattering.

Stark, D. J. (ORCID:0000000203642733)↗

Quantitative spatiotemporal density evolution of aluminum heated purely by monochromatic electrons

A spatially resolved air-wedge shearing interferometer and shadowgraph diagnostic provides measurements of electron density with a resolution of ∼40 μm. A ∼100-ns-long, monoenergetic electron bunch at 19.8 MeV and a current of 1.4 kA (8.5×1014 e−) heats 100-μm-thick aluminum (Al) foils in a 1-mm-spot to Te∼1 eV. A 5-ns-long, ∼60 mJ, frequency doubled Nd:YAG laser probes the dense Al plasma. Electron densities up to 1020cm−3 are resolved; the maximum resolvable density is limited by opacity, transmission, and spatial fringe density achievable with the detector. This diagnostic provides measurements of the total phase shift, transmission, and electron density. Several measurements at different time slices provide the ability to determine the velocity of the leading edge of the shadowgraph and compare it to the motion of different density shells. These measurements are also compared to radiation hydrodynamics simulations. A rough quantitative agreement is shown between the hydro simulations and the measurements; there are differences in the exact density distributions.

Coleman, J. E. (ORCID:0000000308873757)↗

Impact of pedestal density gradient and collisionality on ELM dynamics

BOUT++ turbulence simulations are conducted to capture the underlying physics of small ELM characteristics achieved by increasing separatrix density via controlling strike points from vertical to horizontal divertor plates for three EAST discharges. BOUT++ linear simulations show that the most unstable modes change from high-n ideal ballooning modes to intermediate-n peeling–ballooning modes and eventually to peeling–ballooning stable plasmas in the pedestal. Nonlinear simulations show that the fluctuation is saturated at a high level for the lowest separatrix density. The ELM size decreases with increasing separatrix density, until the fraction of this energy lost during the ELM crash becomes less than 1% of the pedestal stored energy, leading to small ELMs. Simulations indicate that small ELMs can be triggered either by the marginally peeling–ballooning instability near the peak pressure gradient position inside the pedestal or by a local instability in the pedestal foot with a larger separatrix density gradient. The pedestal collisionality scan for type-I ELMs with steep pedestal density gradient shows that both linear growth rate and ELM size decrease with increasing collisionality. On the contrary, the pedestal collisionality and pedestal density width scan with a weak pedestal density gradient indicate small ELMs can either be triggered by a high-n ballooning mode or by a low-n peeling mode in a low collisionality region 0.04–0.1. We report the simulations indicate the weaker the linear unstable modes near marginal stability with small linear growth rate, the lower nonlinearly saturated fluctuation intensity and the smaller turbulence spreading from the linear unstable zone to stable zone in the nonlinear saturation phase, leading to small ELMs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Al–W gradient density materials—Processing and dynamic ramp compression

Materials with high-density gradients are desired for controlling loading paths in dynamic compression, important for studying material properties in extreme conditions and inertial confinement fusion. The large density difference between Al and W makes them ideal choices for producing gradient density materials, but their extremely different melting temperatures make them challenging to fabricate simultaneously. We report a method for producing Al–W porosity-free materials with a fourfold increase in density (2.7–11 g/cm 3 ) across the composition range, from Al-rich to W-rich, without intermetallic phase formation. This was achieved by understanding the aluminum-dominated densification behavior and examining the influence of pressure and temperature on the densification of Al–W composites. Dynamic compression experiments conducted with the Al–W gradient density material produced shock ramp compressions as expected based on the designed composition, and the performed hydrodynamics simulations showed excellent agreement with experimental results. The results demonstrate that current activated pressure-assisted densification allows for the easy and rapid fabrication of gradient density materials with significant density gradients and tailored compositions, facilitating precise control of the loading paths. These materials have the potential to create customized pressure drives for advancing the fields of material science in extreme environments and dynamic compression.

Alloys↗

Medium-density amorphous ice unveils shear rate as a new dimension in water’s phase diagram

Recent experiments revealed a new amorphous ice phase, medium-density amorphous ice (MDA), formed by ball-milling ice I h at 77 K [Rosu-Finsen et al., Science 379, 474–478 (2023)]. MDA has density between that of low-density amorphous (LDA) and high-density amorphous (HDA) ices, adding to the complexity of water’s phase diagram, known for its glass polyamorphism and two-state thermodynamics. The nature of MDA and its relation to other amorphous ices and liquid water remain unsolved. Here, we use molecular simulations under controlled pressure and shear rate at 77 K to produce and investigate MDA. Here. we find that MDA formed at constant shear rate is a steady-state nonequilibrium shear-driven amorphous ice (SDA), that can be produced by shearing ice I h , LDA, or HDA. Our results suggest that MDA could be obtained by ball-milling water glasses without crystallization interference. Increasing the shear rate at ambient pressure produces SDAs with densities ranging from LDA to HDA, revealing shear rate as a new thermodynamic variable in the nonequilibrium phase diagram of water. Indeed, shearing provides access to amorphous states inaccessible by controlling pressure and temperature alone. SDAs produced with shearing rates as high as 10 6 s −1 sample the same region of the potential energy landscape than hyperquenched glasses with identical density, pressure, and temperature. Intriguingly, SDAs obtained by shearing at ~10 8 s −1 have density, enthalpy, and structure indistinguishable from those of water “instantaneously” quenched from room temperature to 77 K over 10 ps, making them good approximants for the “true glass” of ambient liquid water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗