SEARCH · Search NASA

Results for “THEORIES - FLOW MIXING LENGTH”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

NASA NTRS1986

Remote sensing of Earth Terrain

The Mueller matrix and polarization covariance matrix for polarimetric radar systems were studied. The clutter is modelled by a layer of random permittivity, described by a three dimensional correlation function, with variance, and horizontal and vertical correlation lengths. This model is applied, using the wave theory with Born approximations carried to the second order, to find the backscattering elements of the polarimetric matrices. It is found that 8 out of 16 elements of the Mueller matrix are identically zero, corresponding to a covariance matrix with four zero elements. Theoretical predictions are matched with experimental data for vegetation fields. A general mixing formula is derived for discrete scatterers immersed in a host medium. The inclusion particles are assumed to be ellipsoidal. The electric field inside the scatterers is determined by quasistatic analysis, assuming the diameter of the inclusion particles to be much smaller than the wavelength. The results are applicable to general multiphase mixtures.

Kong, J. A.↗

DOE OSTI2024

Phase-space entropy cascade and irreversibility of stochastic heating in nearly collisionless plasma turbulence

We consider a nearly collisionless plasma consisting of a species of “test particles” in one spatial and one velocity dimension, stirred by an externally imposed stochastic electric field—a kinetic analog of the Kraichnan model of passive advection. The mean effect on the particle distribution function is turbulent diffusion in velocity space—known as stochastic heating. Accompanying this heating is the generation of fine-scale structure in the distribution function, which we characterize with the collisionless (Casimir) invariant C 2 ∝ ∫ ∫ d x d v 〈 f 2 〉 —a quantity that here plays the role of (negative) entropy of the distribution function. We find that C 2 is transferred from large scales to small scales in both position and velocity space via a phase-space cascade enabled by both particle streaming and nonlinear interactions between particles and the stochastic electric field. We compute the steady-state fluxes and spectrum of C 2 in Fourier space, with k and s denoting spatial and velocity wave numbers, respectively. In our model, the nonlinearity in the evolution equation for the spectrum turns into a fractional Laplacian operator in k space, leading to anomalous diffusion. Whereas even the linear phase mixing alone would lead to a constant flux of C 2 to high s (towards the collisional dissipation range) at every k , the nonlinearity accelerates this cascade by intertwining velocity and position space so that the flux of C 2 is to both high k and high s simultaneously. Integrating over velocity (spatial) wave numbers, the k -space ( s -space) flux of C 2 is constant down to a dissipation length (velocity) scale that tends to zero as the collision frequency does, even though the rate of collisional dissipation remains finite. The resulting spectrum in the inertial range is a self-similar function in the ( k , s ) plane, with power-law asymptotics at large k and s . Our model is fully analytically solvable, but the asymptotic scalings of the spectrum can also be found via a simple phenomenological theory whose key assumption is that the cascade is governed by a “critical balance” in phase space between the linear and nonlinear timescales. We argue that stochastic heating is made irreversible by this entropy cascade and that, while collisional dissipation accessed via phase mixing occurs only at small spatial scales rather than at every scale as it would in a linear system, the cascade makes phase mixing even more effective overall in the nonlinear regime than in the linear one. Published by the American Physical Society 2024

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Five percent measurement of the gravitational constant in the Large Magellanic Cloud

We perform a novel test of General Relativity by measuring the gravitational constant in the Large Magellanic Cloud (LMC). The LMC contains six well-studied Cepheid variable stars in detached eclipsing binaries. Radial velocity and photometric observations enable a complete orbital solution, and precise measurements of the Cepheids’ periods permit detailed stellar modelling. Both are sensitive to the strength of gravity, the former via Kepler’s third law and the latter through the gravitational free-fall time. We jointly fit the observables for stellar parameters and the gravitational constant. Performing a full Markov Chain Monte Carlo analysis of the parameter space including all relevant nuisance parameters, we constrain the gravitational constant in the Large Magellanic Cloud relative to the Solar System to be G LMC /G SS = $0.93{^{+0.05}_{-0.04}}$. We discuss the implications of this 5% measurement of Newton’s constant in another galaxy for dark energy and modified gravity theories. This result excludes one Cepheid, CEP-1812, which is an outlier and needs further study: it is either a highly unusual system to which our model does not apply, or it prefers G LMC < G SS at 2.6σ. We also obtain new bounds on critical parameters that appear in semi-analytic descriptions of stellar processes. In particular, we measure the mixing length parameter to be $α = 0.90{^{+0.36}_{-0.26}}$ (when assumed to be constant across our sample), and obtain constraints on the parameters describing turbulent dissipation and convective flux.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

LLNL FESP Theory Highlights: March 2024

Predictions of heat load widths λ q based on particle orbits alone are very pessimistic. This paper shows that pedestal peeling-ballooning (P-B) MHD turbulence broadens the stable SOL by the transport, or spreading, of fluctuation energy from the pedestal. λ q is seen to increase with Γ ε , the fluctuation energy density flux. We elucidate the fundamental physics of the spreading process. Γ ε increases with pressure fluctuation correlation length. P-B turbulence is seen to be especially effective at spreading, on account of its large effective mixing length. Spreading is shown to be a multiscale process, which is enhanced by the synergy of large and small-scale modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Convective shells in the interior of Cepheid variable stars: Overshooting models based on hydrodynamic simulations

Context. Because Cepheid variable stars have long been used as a cosmic benchmark for scaling distances in our Galaxy and beyond, the accuracy of stellar evolution models for Cepheids have wide-reaching effects. However, our understanding of the dynamics in the interiors of these physically complex stars is limited. Aims. Our goal is to provide a detailed multi-dimensional picture of hydrodynamic convection and convective boundary mixing in the interior of Cepheids. Methods. Using the Modules for Experiments in Stellar Astrophysics (MESA), we studied the structure of intermediate-mass stars that cross the instability strip. Then, we performed two-dimensional hydrodynamic simulations of six stars with the fully compressible Multidimensional Stellar Implicit Code (MUSIC). Our simulations did not model the radial pulsations but focused on the interior structure of this family of stars. We developed and applied a new statistical analysis to examine convection and convective boundary mixing in the interior of these stellar simulations. Results. Based on a grid of MESA models, we demonstrated that a common structure for intermediate mass Cepheids includes an interior convective shell as well as a thin outer convective envelope. Using the extreme value theory approach to analyze our MUSIC simulation data, we found that overshooting above the convective shell fills the space between these convectively unstable layers. We developed a new statistical analysis that provides a clearer picture of how overshooting fills this layer; it also allowed us to formulate a detailed comparison between overshooting above and below the convective shell. Our analysis effectively decomposes the overshooting layer into two layers: a weak overshooting layer and a strong overshooting layer. Statistically, this is accomplished by decomposing the strongly non-Gaussian probability density function into a mixture of gamma distributions. Using our mixture model, we showed that the ratio of overshooting lengths above and below the convective shell depends directly on the radial extent of the convective shell as well as its depth in the star. We proposed a new form for the diffusion coefficient that addresses the need for overlapping overshooting layers between convective shells. We introduced the idea of a “super-mixing layer” where overshooting from both the convective shell and the convective envelope results in efficient mixing and could be viewed as merging the two adjacent convective zones.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1987

Interlaminar shear fracture toughness and fatigue thresholds for composite materials

Static and cyclic end notched flexure tests were conducted on a graphite epoxy, a glass epoxy, and graphite thermoplastic to determine their interlaminar shear fracture toughness and fatigue thresholds for delamination in terms of limiting values of the mode II strain energy release rate, G-II, for delamination growth. The influence of precracking and data reduction schemes are discussed. Finite element analysis indicated that the beam theory calculation for G-II with the transverse shear contribution included was reasonably accurate over the entire range of crack lengths. Cyclic loading significantly reduced the critical G-II for delamination. A threshold value of the maximum cyclic G-II below which no delamination occurred after one million cycles was identified for each material. Also, residual static toughness tests were conducted on glass epoxy specimens that had undergone one million cycles without delamination. A linear mixed-mode delamination criteria was used to characterize the static toughness of several composite materials; however, a total G threshold criterion appears to characterize the fatigue delamination durability of composite materials with a wide range of static toughness.

Obrien, T. Kevin↗

NASA NTRS1989

Interlaminar shear fracture toughness and fatigue thresholds for composite materials

O'Brien, T. Kevin↗

NASA NTRS1994

On sound generation by turbulent convection: A new look at old results

We have revisited the problem of acoustic wave generation by turbulent convection in stellar atmospheres. The theory of aerodynamically generated sound, originally developed by Lighthill and later modified by Stein to include the effects of stratification, has been used to estimate the acoustic wave energy flux generated in solar and stellar convection zones. We correct the earlier computations by incorporating an improved description of the spatial and temporal spectrum of the turbulent convection. We show the dependence of the resulting wave fluxes on the nature of the turbulence, and compute the wave energy spectra and wave energy fluxes generated in the Sun on the basis of a mixing-length model of the solar convection zone. In contrast to the previous results, we show that the acoustic energy generation does not depend very sensitively on the turbulent energy spectrum. However, typical total acoustic fluxes of order F(sub A) = 5 x 10(exp 7) ergs/sq cm/s with a peak of the acoustic frequency spectrum near omega = 100 mHz are found to be comparable to those previously calculated. The acoustic flux turns out to be strongly dependent on the solar model, scaling with the mixing-length parameter alpha as alpha(exp 3.8). The computed fluxes most likely constitute a lower limit on the acoustic energy produced in the solar convection zone if recent convection simulations suggesting the presence of shocks near the upper layers of the convection zone apply to the Sun.

Musielak, Z. E.↗

DOE OSTI2020

Dispersity-Driven Stabilization of Coexisting Morphologies in Asymmetric Diblock Copolymer Thin Films

Despite decades of research, characterization of the effects of polymer chain dispersity on the structural properties of block copolymer thin films remains challenging. Here we present an integrated experimental and modeling approach to characterize the morphology of thin films containing asymmetric diblock copolymers. Specifically, we used synergistic neutron reflectivity (NR) and self-consistent field theory (SCFT)-based modeling to realize unexpected morphology of thin films containing asymmetric copolymers. Using NR, a highly stable and reproducible mixed phase of coexisting cylinders and lamellar domains was discovered in asymmetric poly(deuterated-styrene-b-n butyl methacrylate) (dPS-PBMA) copolymer thin films containing 34% volume fraction of dPS. SCFT reveals how to obtain such a thermodynamically stable morphology in the presence of disperse majority block and asymmetric interactions of polymer species with surfaces. Stabilization of the coexisting domains is a consequence of the depth segregation based on chain-length distribution. The asymmetric chains microphase-separate into cylindrical domains close to the substrate and near-symmetric chains form lamellar domains at the air interface. In the absence of dispersity, the coexistence of cylindrical and lamellar domains is thermodynamically unstable because of the absence of depth segregation. Overall, such an effect of dispersity on diblock copolymer thin-film morphology reveals a unique and powerful strategy to create coexisting nanoscale domains and tailor properties of thin films.

36 MATERIALS SCIENCE↗

NASA NTRS1976

Heat transfer to a full-coverage film-cooled surface with 30 degree slant-hole injection

Heat transfer behavior was studied in a turbulent boundary layer with full coverage film cooling through an array of discrete holes and with injection 30 deg to the wall surface in the downstream direction. Stanton numbers were measured for a staggered hole pattern with pitch-to-diameter ratios of 5 and 10, an injection mass flux ratio range of 0.1 to 1.3, and a range of Reynolds number Re sub x of 150,000 to 5 million. Air was used as the working fluid, and the mainstream velocity varied from 9.8 to 34.2 m/sec (32 to 112 ft/sec). The data were taken for secondary injection temperature equal to the wall temperature and also equal to the mainstream temperature. The data may be used to obtain Stanton number as a continuous function of the injectant temperature by use of linear superposition theory. The heat transfer coefficient is defined on the basis of a mainstream-to-wall temperature difference. This definition permits direct comparison of performance between film cooling and transpiration cooling. A differential prediction method was developed to predict the film cooling data base. The method utilizes a two-dimensional boundary layer program with routines to model the injection process and turbulence augmentation. The program marches in the streamwise direction, and when a row of holes is encountered, it stops and injects fluid into the boundary layer. The turbulence level is modeled by algebraically augmenting the mixing length, with the augmentation keyed to a penetration distance for the injected fluid.

Crawford, M. E.↗

DOE OSTI2026

Custom-trained Machine-learning Interatomic Potentials: ZnCl2 Aqueous Solution

This dataset was generated using an iterative active-learning strategy implemented in the ArcaNN software package (https://github.com/arcann-chem/arcann_training) to train machine-learning interatomic potentials for aqueous ZnCl2 solutions. Each active-learning cycle consisted of three stages: training, exploration, and labeling. The initial training set combined configurations generated in this work from enhanced-sampling ab initio molecular dynamics simulations with configurations from a previously reported neural-network-potential study of aqueous ZnCl2. The enhanced-sampling ab initio molecular dynamics simulations involved Zn–Cl separation and the chloride coordination number around Zn²? as collective variables. These configurations served as the seed dataset. Subsequent active-learning cycles expanded the training set by identifying and labeling configurations that were poorly represented by the current models, thereby improving coverage of ion-association states and changes in local coordination and charge-state environments relevant to the solution free-energy landscape. For all selected configurations, single-point calculations of the total energies and atomic forces were performed within density functional theory using the CP2K Quickstep module. Reference calculations employed the revPBE-D3 and r2SCAN exchange-correlation functionals. Motivated by recent work on aqueous Zn²?, the main revPBE calculations omitted D3 dispersion contributions involving Zn²?, while retaining the D3 correction for water and chloride. For comparison, fully dispersion-corrected revPBE-D3 reference calculations were also performed, with D3 applied to all species, including Zn²?. Valence electrons were treated explicitly, while core electrons were represented using norm-conserving Goedecker–Teter–Hutter pseudopotentials. The wave functions were expanded using the mixed Gaussian-and-plane-wave scheme with TZV2P-MOLOPT basis sets for all elements and a 600 Ry auxiliary plane-wave cutoff for the electron density. Self-consistent-field convergence was accelerated using the orbital-transformation and Direct Inversion in the Iterative Subspace algorithms, with a convergence threshold of 10?6. All single-point calculations were performed in periodic orthorhombic cells. The CELL_REF keyword in CP2K was used to define a fixed reference cell with a box length of 25 Å. This treatment ensured a consistent reference for configurations extracted from NpT trajectories with fluctuating cell dimensions. The resulting DFT energies and atomic forces constitute the ground-truth labels used to train the MLIPs. The resulting MLIP was trained for aqueous ZnCl2 solutions spanning concentrations from 0 to 30 molal and a broad pH range, from strongly acidic to strongly basic conditions. Representative examples of configurations included in the MLIP training dataset are provided below. These include 1) Representative configurations from the dataset labeled at the revPBE-D3 level, with D3 dispersion interactions involving Zn2+ excluded (revPBE-wo-D3). 2) Representative configurations from the dataset labeled at the fully dispersion-corrected revPBE-D3 level, with D3 interactions applied to all species, including Zn2+ (revPBE-D3). 3) Representative configurations from the dataset labeled at the r2SCAN level of theory (r2SCAN).

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗

NASA NTRS1983

Fracture mechanics for delamination problems in composite materials

A fracture mechanics approach to the well-known delamination problem in composite materials is presented. Based on the theory of anisotropic laminate elasticity and interlaminar fracture mechanics concepts, the composite delamination problem is formulated and solved. The exact order of the delamination crack-tip stress singularity is determined. Asymptotic stress and displacement fields for an interlaminar crack are obtained. Fracture mechanics parameters such as mixed-mode stress intensity factors, KI, KII, KIII, and the energy release rate, G, for composite delamination problems are defined. To illustrate the fundamental nature of the delamination crack behavior, solutions for edge-delaminated graphite-epoxy composites under uniform axial extension are presented. Effects of fiber orientation, ply thickness, and delamination length on the interlaminar fracture are examined.

Wang, S. S.↗

DOE OSTI2025

The role of intermediate ΔΔ states in nucleon–nucleon scattering in the large-Nc and unitary limits, and ΔΔ and ΩΩ scattering

We explore potential explanations for why using large-Nc (Nc is the number of colors) scaling to determine the relative size of few-nucleon low-energy operators agrees with experiment even when dynamical Δ’s are not explicitly included. Given that the large-Nc analysis is predicated on the nucleons and Δ’s being degenerate, this is a curious result. We show that for purely S-wave interactions the relationships dictated by large-Nc scaling are unaffected whether the Δ is included or not. In the case of higher partial waves that do not mix with S-waves, the impact of the Δ is perturbative, which makes the agreement with naive (Δ-less) large-Nc ordering unsurprising. For higher partial waves that mix with S-waves, the nucleon and Δ would need to decouple to get agreement with naive large-Nc ordering. We find all NN, ΔN, and ΔΔ low energy coefficients for leading-order baryon–baryon scattering in Δ-full pionless effective field theory in terms of the two independent parameters dictated by the SU(2F) spin-flavor symmetry that arises in the Nc → ∞ limit. Because of recent lattice quantum chromodynamics results and experimental interest, we extend our analysis to the three-flavor case to study ΩΩ scattering. We show that in the unitary limit (where scattering lengths become infinite) one of the two SU(2F) parameters is driven to zero, resulting in enhanced symmetries, which agree with those found in spin-1/2 entanglement studies.

Richardson, Thomas R↗

DOE OSTI2024

VOTCA: multiscale frameworks for quantum and classical simulations in soft matter

Many physical phenomena in liquids and soft matter are multiscale by nature and can involve processes with quantum and classical degrees of freedom occurring over a vast range of length- and timescales. Examples range from structure formation processes of complex polymers or even polymer blends (Svaneborg & Everaers, 2023) on the classical side to charge and energy transport and conversion processes (Lee et al., 2019) involving explicit electronic and, therefore, quantum information. The Versatile Object-oriented Toolkit for Coarse-graining Applications (VOTCA) provides multiscale frameworks built on a comprehensive set of methods for the development of classical coarse-grained potentials (VOTCA-CSG) as well as state-of-the art excited state electronic structure methods based on density-functional and many-body Green’s function theories, coupled in mixed quantum-classical models and used in kinetic network models (VOTCA-XTP).

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Fluid Mixing during Phase Separation in Crystallizing White Dwarfs

Abstract Accurate models of cooling white dwarfs must treat the energy released as their cores crystallize. This phase transition slows the cooling by releasing latent heat and also gravitational energy, which results from phase separation: liquid C is released from the solid C/O core, driving an outward carbon flux. The Gaia color–magnitude diagram provides striking confirmation of this theory by revealing a mass-dependent overdensity of white dwarfs, indicating slowed cooling at the expected location. However, the observed overdensity is enhanced relative to the models. Additionally, it is associated with increased magnetism, suggesting a link between crystallization and magnetic field generation. Recent works aimed at explaining an enhanced cooling delay and magnetic field generation employ a uniform mixing prescription that assumes large-scale turbulent motions; we show here that these calculations are not self-consistent. We also show that thermohaline mixing is most likely efficient enough to provide the required chemical redistribution during C/O phase separation, and that the resulting velocities and mixing lengths are much smaller than previous estimates. These reduced fluid motions cannot generate measurable magnetic fields, suggesting any link with crystallization needs to invoke a separate mechanism. Finally, this mixing alters the chemical profiles, which in turn affects the frequencies of the pulsation modes.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Morphology and Transport of Multivalent Cation-Exchanged Ionomer Membranes Using Perfluorosulfonic Acid–Ce Z+ as a Model System

Perfluorosulfonic acids (PFSAs) are commonly used as solid polymer electrolyte membranes (PEMs) in electrochemical energy devices, where they are vulnerable to attack by radical species during operation. A popular strategy to combat this problem is to introduce radical scavengers like cerium (Ce) ions that neutralize these species before they attack the PFSA. Such cation doping creates a multi-ion system, in which understanding the mechanisms of cation solvation and transport becomes important for the effective design and utilization of PFSA–cation systems. Ce ions also provide a representative model system for multication-exchanged ionomers in electrochemical systems. In this study, hydration and conductivity measurements, along with X-ray fluorescence and scattering, are employed to elucidate how Ce ion exchange alters PFSA’s ionic solvation, as well as nano- and mesoscale morphologies, which ultimately control its ion transport properties. A molecular transport model is used to deconvolute the impact of Ce ions on the local solvation structure of water in the membrane from mesoscale changes of the transport pathways. The combined experimental and theoretical analysis reveals a nonlinear decrease in conductivity driven by cation solvation at the molecular level and morphological changes at longer length scales. Migration–diffusion coupling, its nonlinear dependence on ion exchange and hydration, and its overall implications for ionomer performance are also discussed. Finally, these findings have the potential to be translated into other mixed cation–ionomer systems for a wide range of energy and environmental devices.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Internal Conversion between Bright (1 1 B $^+_u$) and Dark (21 A $^-_g$) States in s- trans -Butadiene and s- trans -Hexatriene

Internal conversion (IC) between the two lowest singlet excited states, 1 1 B$_u^+$ and 2 1 A$_g^–$, of s-trans-butadiene and s-trans-hexatriene is investigated using a series of single- and multi- reference wave function and density functional theory (DFT) methodologies. Three independent types of the equation-of-motion coupled-cluster (EOMCC) theory capable of providing an accurate and balanced description of one- as well as two-electron transitions, abbreviated as δ-CR-EOMCC(2,3), DIP-EOMCC(4h2p){N o }, and DEA-EOMCC(4p2h){N u } or DEA-EOMCC(3p1h,4p2h){N u }, consistently predict that the 1 1 B$_u^+$/2 1 A$_g^–$ crossing in both molecules occurs along the bond length alternation coordinate. However, the analogous 1 1 B$_u^+$ and 2 1 A$_g^–$ potentials obtained with some multireference approaches, such as CASSCF and MRCIS(D), as well as with the linear-response formulation of time-dependent DFT (TDDFT), do not cross. Hence, caution needs to be exercised when studying the low-lying singlet excited states of polyenes with conventional multiconfigurational methods and TDDFT. The multistate many-body perturbation theory methods, such as XMCQDPT 2 , do correctly reproduce the curve crossing. Among the simplest and least expensive computational methodologies, the DFT approaches that incorporate the contributions of doubly excited configurations, abbreviated as MRSF (mixed reference spin-flip) TDDFT and SSR( 4,4 ), accurately reproduce our best EOMCC results. This is highly promising for nonadiabatic molecular dynamics simulations in larger systems.

36 MATERIALS SCIENCE↗

NASA NTRS1981

Edge delamination in angle-ply composite laminates, part 5

A theoretical method was developed for describing the edge delamination stress intensity characteristics in angle-ply composite laminates. The method is based on the theory of anisotropic elasticity. The edge delamination problem is formulated using Lekhnitskii's complex-variable stress potentials and an especially developed eigenfunction expansion method. The method predicts exact orders of the three-dimensional stress singularity in a delamination crack tip region. With the aid of boundary collocation, the method predicts the complete stress and displacement fields in a finite-dimensional, delaminated composite. Fracture mechanics parameters such as the mixed-mode stress intensity factors and associated energy release rates for edge delamination can be calculated explicity. Solutions are obtained for edge delaminated (theta/-theta theta/-theta) angle-ply composites under uniform axial extension. Effects of delamination lengths, fiber orientations, lamination and geometric variables are studied.

Wang, S. S.↗