Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “statistical mechanics models”

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 901 records · Page 50

Computational thermodynamics and its applications

Thermodynamics is a science concerning the state of a system, whether it is stable, metastable or un- stable, when interacting with the surroundings. In this overview, fundamentals of thermodynamics are briefly reviewed through the combination of first and second laws of thermodynamics for open and nonequilibrium systems to demonstrate that the reversible equilibrium and irreversible nonequilibrium thermodynamics can be integrated to enhance the power and utilities of thermodynamics. The recent progresses in computational thermodynamics, the remaining challenges, and potential impacts in broad scientific fields are discussed here. It is shown that computational thermodynamics enables the modeling of thermodynamics of a state as a function of both external and internal variables and quantitative calculations of a broad range of properties of a multicomponent system in terms of first and second derivatives of energy, including not only equilibrium states when there are no driving forces for any internal processes and but also non-equilibrium states with driving forces for internal processes. Consequently, external constraints such as fixed strain and internal degree of freedoms such as ordering and defects can be described in a coherent framework and applied to materials design. Furthermore, two important but largely overlooked aspects in thermodynamics will be discussed, i.e. the rigorous application of statistical thermodynamics with the probability of configurations and their contributions to system properties, and the applications of second derivatives of energy with respect to either two extensive variables or two potentials or a mixture of them in terms of understanding and predicting emergent behaviors, critical phenomena, kinetic coefficients, and mechanical properties.

36 MATERIALS SCIENCE↗

Tea consumption and risk of bladder cancer in the Bladder Cancer Epidemiology and Nutritional Determinants (BLEND) Study: Pooled analysis of 12 international cohort studies

Tea has been shown to be associated with reduced risk of several diseases including cardiovascular diseases, stroke, metabolic syndrome, and obesity. However, the results on the relationship between tea consumption and bladder cancer are conflicting. This research aimed to assess the association between tea consumption and risk of bladder cancer using a pooled analysis of prospective cohort data. Individual data from 532,949 participants in 12 cohort studies, were pooled for analyses. Cox regression models stratified by study centre was used to estimate hazard ratios (HR) and corresponding 95% CIs. Fractional polynomial regression models were used to examine the dose–response relationship. A higher level of tea consumption was associated with lower risk of bladder cancer incidence (compared with no tea consumption: HR = 0.87, 95% C.I. = 0.77–0.98 for low consumption; HR = 0.86, 95% C.I. = 0.77–0.96 for moderate consumption; HR = 0.84, 95% C.I. = 0.75–0.95 for high consumption). When stratified by sex and smoking status, this reduced risk was statistically significant among men and current and former smokers. In addition, dose–response analyses showed a lower bladder cancer risk with increment of 100 ml of tea consumption per day (HR-increment = 0.97; 95% CI = 0.96–0.98). A similar inverse association was found among males, current and former smokers while never smokers and females showed non-significant results, suggesting potential sex-dependent effect. Higher consumption of tea is associated with reduced risk of bladder cancer with potential interaction with sex and smoking status. Further studies are needed to clarify the mechanisms for a protective effect of tea (e.g. inhibition of the survival and proliferation of cancer cells and anti-inflammatory mechanisms) and its interaction with smoking and sex.

60 APPLIED LIFE SCIENCES↗

Wide-Field Hubble Space Telescope Observations of the Globular Cluster System in NGC 1399*

We present a comprehensive high spatial resolution imaging study of globular clusters (GCs) in NGC 1399, thecentral giant elliptical cD galaxy in the Fornax galaxy cluster, conducted with the Advanced Camera for Surveys(ACS) aboard theHubble Space Telescope(HST).Using a novel technique to construct drizzled point-spreadfunction libraries for HSTACS data, we accurately determine the fidelity of GC structural parameter measurementsfrom detailed artificial star cluster experiments and show the superior robustness of the GC half-light radius,rh,compared with other GC structural parameters, such as King core and tidal radius. The measurement ofrhfor themajor fraction of the NGC 1399 GC system reveals a trend of increasingrhversus galactocentric distance,Rgal,out to about 10 kpc and a flat relation beyond. This trend is very similar for blue and red GCs, which are found tohave a mean size ratio ofrh,redrh,blue0.820.11 at all galactocentric radii from the core regions of the galaxyout to40 kpc. This suggests that the size difference between blue and red GCs is due to internal mechanismsrelated to the evolution of their constituent stellar populations. Modeling the mass density profile of NGC 1399shows that additional external dynamical mechanisms are required to limit the GC size in the galaxy halo regionstorh2 pc. We suggest that this may be realized by an exotic GC orbit distribution function, an extended darkmatter halo, andor tidal stress induced by the increased stochasticity in the dwarf halo substructure at largergalactocentric distances. We compare our results with the GCrhdistribution functions in various galaxies and findthat the fraction of extended GCs withrh5 pc is systematically larger in late-type galaxies compared with GCsystems in early-type galaxies. This is likely due to the dynamically more violent evolution of early-type galaxies.We match our GCrhmeasurements with radial velocity data from the literature and split the resulting sample at themedianrhvalue into compact and extended GCs. We find that compact GCs show a significantly smaller line-of-sight velocity dispersion,cmp22525 km s1, than their extended counterparts,ext31721 km s1.Considering the weaker statistical correlation in the GCrhcolor and the GCrhRgalrelations, the more significantGC sizedynamics relation appears to be astrophysically more relevant and hints at the dominant influence of theGC orbit distribution function on the evolution of GC structural parameters.

Hubble↗

Why it is Unfortunate that Linear Machine Learning “Works” so well in Electromechanical Switching of Ferroelectric Thin Films

Machine learning (ML) is relied on for materials spectroscopy. It is challenging to make ML models fail because statistical correlations can mimic the physics without causality. Here, using a benchmark band-excitation piezoresponse force microscopy polarization spectroscopy (BEPS) dataset the pitfalls of the so-called “better”, “faster”, and “less-biased” ML of electromechanical switching are demonstrated and overcome. Using a toy and real experimental dataset, it is demonstrated how linear nontemporal ML methods result in physically reasonable embedding (eigenvalues) while producing nonsensical eigenvectors and generated spectra, promoting misleading interpretations. A new method of unsupervised multimodal hyperspectral analysis of BEPS is demonstrated using long-short-term memory (LSTM) β-variational autoencoders (β-VAEs) . By including LSTM neurons, the ordinal nature of ferroelectric switching is considered. Further, to improve the interpretability of the latent space, a variational Kullback–Leibler-divergency regularization is imposed . Finally, regularization scheduling of β as a disentanglement metric is leveraged to reduce user bias. Combining these experiment-inspired modifications enables the automated detection of ferroelectric switching mechanisms, including a complex two-step, three-state one. Ultimately, this work provides a robust ML method for the rapid discovery of electromechanical switching mechanisms in ferroelectrics and is applicable to other multimodal hyperspectral materials spectroscopies.

36 MATERIALS SCIENCE↗

Screening of metal complexes and organic solvents using the COSMOSAC-LANL model to enhance the energy density in a non-aqueous redox flow cell: an insight into the solubility

In this paper, we have proposed a first-principles methodology to screen transition metal complexes against a particular organic solvent and organic solvents against a particular transition metal complex based on their solubility information without the knowledge of heat of fusion and melting temperature. The energy density of a non-aqueous redox flow cell directly depends on the solubility of the redox active species in the non-aqueous medium. Here, we have used the “COSMOSAC-LANL” activity coefficient model which is based on first-principles COSMO calculations where the microscopic information is passed to the macroscopic world via a dielectric continuum solvation model, followed by a post-statistical thermodynamic treatment of the self-consistent properties of the solute particle to calculate the solubility. To model the activity coefficient at infinite dilution for the binary mixtures, a 3-suffix Margules (3sM) function is introduced for the quantitative estimation of the asymmetric interactions and, for the combinatorial term, the Staverman–Guggenheim (SG) form is used. The new activity coefficient model is separately called the “LANL” activity coefficient model. The metal complex and the organic solvent have been treated as a simple binary mixture. The present model has been applied to a set of 14 different organic solvents and 16 different transition metal complexes. Using the new LANL activity coefficient model in combination with the ADF-COSMOSAC-2013 model, we have shown how one can improve the solubility of a transition metal complex in an organic solvent. We applied our model to screen 84 binary mixtures to predict the compatible pair of redox active species and organic solvent to increase the energy density. The solvation mechanism of the transition metal complexes in the organic solvents was obtained using the new model. The results have been compared with the experimental and theoretical results where they are available.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Advanced Laboratory and Field Arrays: Debris Modeling, Detection,& Mitigation (Task 1)

The statement of project objectives for this task was: develop tools, methods and models to assess, and mitigate the risk of damage to MHK infrastructure from woody debris. Develop the capability to detect woody debris using sonar and/or physical methods for purposes of characterizing debris statistics in river (at UAF’s Tanana River Test Site) and near-shore wave (at Yakutat, AK) environments and to activate debris mitigation measures. Develop debris impact risk maps and tables using statistics on debris size, geometry, type, prevalence, mobility and location. Improve and apply the COUPi discrete element method (DEM) to develop models of debris movement and impact on MHK infrastructure to evaluate risk of damage, and interference, to operations from debris. The proposed final deliverable for the task was a set of tools or techniques for providing estimates of the probability of debris impact, and resulting impact forces, on MHK infrastructure as a function of debris size, type, wave regime, and current velocity. Such estimates are required to assess damage risks to operational MHK infrastructure.

13 HYDRO ENERGY↗

Detailed study of oxidation/wear mechanism in lox turbopump bearings

Wear of 440C angular contact ball bearings of the phase 2 high pressure oxygen turbopump (HPOTP) of the space shuttle main engine (SSME) has been studied by means of various advanced nondestructive techniques (NDT) and modeled with reference to all known material, design, and operation variables. Three modes dominating the wear scenario were found to be the adhesive/sheer peeling (ASP), oxidation, and abrasion. Bearing wear was modeled in terms of the three modes. Lacking a comprehensive theory of rolling contact wear to date, each mode is modeled after well-established theories of sliding wear, while sliding velocity and distance are related to microsliding in ball-to-ring contacts. Microsliding, stress, temperature, and other contact variables are evaluated with analytical software packages of SHABERTH(TM)/SINDA(TM) and ADORE(TM). Empirical constants for the models are derived from NIST experiments by applying the models to the NIST wear data. The bearing wear model so established precisely predicts quite well the average ball wear rate for the HPOTP bearings. The wear rate has been statistically determined for the entire population of flight and development bearings based on Rocketdyne records to date. Numerous illustrations are given.

Chase, T. J.↗

Experimental measurements of fluid flow in an 84-pin hexagonal rod bundle with spacer grid for a gas-cooled fast modular reactor

A 50 MW e helium-cooled fast modular reactor (FMR) is under development, and as part of the Department of Energy Integrated Research Project (IRP), Texas A&M University is conducting the thermal-hydraulic characterization of its baseline core configuration. Here, we experimentally investigated the axial and cross flow fields characteristics in the hexagonal fuel rod bundle composed of 84 rods, a central rod, and spacer grids at a Reynolds number of 12,000. A fully transparent experimental facility resembling of one unit of the fuel assembly was constructed. Time-resolved particle image velocimetry (TR-PIV) measurements were performed to characterize hydraulic behavior downstream the spacer grid. Velocity measurements were conducted in the axial and radial direction of the rod bundle. From the PIV velocity vector fields, the full-field flow statistics were computed for the mean velocity, vorticity, and Reynolds stresses. An analysis of the energy decay downstream of the spacer grid was conducted with calculations of secondary-flow intensities, turbulent kinetic energy, power spectrum, and spatial-temporal velocity cross correlations. The vorticity field was obtained and pairs of counter-rotating vortexes were identified using a 3D reconstruction of several measurement planes. The data from this experimental campaign will be used to validate numerical models based on Computational Fluid Dynamics and other codes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Time reversal symmetry breaking and odd viscosity in active fluids: Green–Kubo and NEMD results

Active fluids, which are driven at the microscale by non-conservative forces, are known to exhibit novel transport phenomena due to the breaking of time reversal symmetry. Recently, Epstein and Mandadapu [arXiv:1907.10041 (2019)] obtained Green-Kubo relations for the full set of viscous coefficients governing isotropic chiral active fluids, including the so-called odd viscosity, invoking Onsager's regression hypothesis for the decay of fluctuations in active non-equilibrium steady states. In this Communication, we test these Green-Kubo relations using molecular dynamics simulations of a canonical model system consisting of actively torqued dumbbells. We find the resulting odd and shear viscosity values from the Green-Kubo relations to be in good agreement with values measured independently through non-equilibrium molecular dynamics flow simulations. This provides a test of the Green-Kubo relations and lends support to the application of the Onsager regression hypothesis in relation to viscous behaviors of active matter systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Data-Driven Discovery of Active Nematic Hydrodynamics

Active nematics can be modeled using phenomenological continuum theories that account for the dynamics of the nematic director and fluid velocity through partial differential equations (PDEs). While these models provide a statistical description of the experiments, the relevant terms in the PDEs and their parameters are usually identified indirectly. We adapt a recently developed method to automatically identify optimal continuum models for active nematics directly from spatio-temporal data, via sparse regression of the coarse-grained fields onto generic low order PDEs. After extensive benchmarking, we apply the method to experiments with microtubule-based active nematics, finding a surprisingly minimal description of the system. Furthermore, our approach can be generalized to gain insights into active gels, microswimmers, and diverse other experimental active matter systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

CO Signatures in Subtropical Convective Clouds and Anvils during CRYSTAL-FACE: An Analysis of Convective Transport and Entrainment using Observations and a Cloud-Resolving Model

Convective systems are an important mechanism in the transport of boundary layer air into the upper troposphere. The Cirrus Regional Study of Tropical Anvils and Cirrus Layers-Florida Area Cirrus Experiment (CRYSTAL-FACE) campaign, in July 2002, was developed as a comprehensive atmospheric mission to improve knowledge of subtropical cirrus systems and their roles in regional and global climate. In situ measurements of carbon monoxide (CO), water vapor (H2Ov), and total water (H2Ot) aboard NASA's WB-57F aircraft and CO aboard the U.S. Navy's Twin Otter aircraft were obtained to study the role of convective transport. Three flights sampled convective outflow on 11, 16 and 29 July found varying degrees of CO enhancement relative to the free troposphere. A cloud-resolving model used the in situ observations and meteorological fields to study these three systems. Several methods of filtering the observations were devised here using ice water content, relative humidity with respect to ice, and particle number concentration as a means to statistically sample the model results to represent the flight tracks. A weighted histogram based on ice water content observations was then used to sample the simulations for the three flights. In addition, because the observations occurred in the convective outflow cirrus and not in the storm cores, the model was used to estimate the maximum CO within the convective systems. In general, anvil-level air parcels contained an estimated 20-40% boundary layer air in the analyzed storms.

Lopez, Jimena P.↗

CO Signatures in Subtropical Convective Clouds and Anvils During CRYSTAL-FACE: An Analysis of Convective Transport and Entertainment Using Observations and a Cloud-Resolving Model

Convective systems are an important mechanism in the transport of boundary layer air into the upper troposphere. The Cirrus Regional Study of Tropical Anvils and Cirrus Layers-Florida Area Cirrus Experiment (CRYSTAL-FACE) campaign, in July 2002, was developed as a comprehensive atmospheric mission to improve knowledge of subtropical cirrus systems and their roles in regional and global climate. In situ measurements of carbon monoxide (CO), water vapor (H20v), and total water (H20t) aboard NASA's . WB-57F aircraft and CO aboard the U.S. Navy's Twin Otter aircraft were obtained to study the role of convective transport. Three flights sampled convective outflow on 11, 16 and 29 July found varying degrees of CO enhancement relative to the fiee troposphere. A cloud-resolving model used the in situ observations and meteorological fields to study these three systems. Several methods of filtering the observations were devised here using ice water content, relative humidity with respect to ice, and particle number concentration as a means to statistically sample the model results to represent the flight tracks. A weighted histogram based on ice water content observations was then used to sample the simulations for the three flights. In addition, because the observations occurred in the convective outflow cirrus and not in the storm cores, the model was used to estimate the maximum CO within the convective systems. In general, anvil-level air parcels contained an estimated 20-40% boundary layer air in the analyzed storms.

Lopez, Jimena P.↗

Multiple-time scales analysis of physiological time series under neural control

We discuss multiple-time scale properties of neurophysiological control mechanisms, using heart rate and gait regulation as model systems. We find that scaling exponents can be used as prognostic indicators. Furthermore, detection of more subtle degradation of scaling properties may provide a novel early warning system in subjects with a variety of pathologies including those at high risk of sudden death.

Non-NASA Center↗

An engineering model of woven composites based on micromechanics

Composites with three-dimensional woven architectures exhibit large strains to failure when compared to composites made up of the same materials but not with three-dimensional interlocking tows. The fracture mechanics of such three-dimensional architectures is a subject requiring substantial investigation and experimental testing. Classical fracture mechanics concepts (for instance, an isolated defect in a homogeneous body) will not be applicable to the woven fracture test specimen. The use of an isolated singularity to characterize an entire specimen is inadequate when the density of defects is considerable and the material is heterogeneous. Modelling of such a complex system requires a great deal of insight and consideration as well as prudent choices of model sizes to make numerical schemes feasible. The purpose of this manuscript is to review our recently acquired knowledge of damage accumulation in woven composites and to describe a practicable model of the macroscopic behavior in these and other complex composite architectures based on such knowledge. In this manuscript, discussion will be limited to uniaxial compressive loading; considerations of general loading (monotonic and cyclic) will appear in a subsequent manuscript. Our modelling efforts may be briefly described as follows: the composite is subdivided into microstructural elements (microelements) in which the micromechanical modelling is either understood rigorously or can be represented adequately by statistical parameters. There can be microstructural elements for many different types of composite components, such as the various types of warp and weft and matrix for three-dimensional woven composites. The physical dimensions of microelements are made as large as possible while the response within the element can still be represented by a single micromechatlical calculation. The various elements are linked together(sometimes by associating distinct corners and edges, sometimes by superposition) in a pattern which resembles a particular weave architecture. The model can then be loaded in any manner and the linear and nonlinear elastic responses of representative weaves can be calculated. After the elastic regime, the fracture response is determined by monitoring the damage accumulation.

Carter, W. C.↗

Turbulent Concentration of mm-Size Particles in the Protoplanetary Nebula: Scale-Dependent Cascades

The initial accretion of primitive bodies (here, asteroids in particular) from freely-floating nebula particles remains problematic. Traditional growth-by-sticking models encounter a formidable "meter-size barrier" (or even a mm-to-cm-size barrier) in turbulent nebulae, making the preconditions for so-called "streaming instabilities" difficult to achieve even for so-called "lucky" particles. Even if growth by sticking could somehow breach the meter size barrier, turbulent nebulae present further obstacles through the 1-10km size range. On the other hand, nonturbulent nebulae form large asteroids too quickly to explain long spreads in formation times, or the dearth of melted asteroids. Theoretical understanding of nebula turbulence is itself in flux; recent models of MRI (magnetically-driven) turbulence favor low-or- no-turbulence environments, but purely hydrodynamic turbulence is making a comeback, with two recently discovered mechanisms generating robust turbulence which do not rely on magnetic fields at all. An important clue regarding planetesimal formation is an apparent 100km diameter peak in the pre-depletion, pre-erosion mass distribution of asteroids; scenarios leading directly from independent nebula particulates to large objects of this size, which avoid the problematic m-km size range, could be called "leapfrog" scenarios. The leapfrog scenario we have studied in detail involves formation of dense clumps of aerodynamically selected, typically mm-size particles in turbulence, which can under certain conditions shrink inexorably on 100-1000 orbit timescales and form 10-100km diameter sandpile planetesimals. There is evidence that at least the ordinary chondrite parent bodies were initially composed entirely of a homogeneous mix of such particles. Thus, while they are arcane, turbulent concentration models acting directly on chondrule size particles are worthy of deeper study. The typical sizes of planetesimals and the rate of their formation can be estimated using a statistical model with properties inferred from large numerical simulations of turbulence. Nebula turbulence is described by its Reynolds number Re = (L/eta)(exp 4/3), where L = H alpha(exp 1/2) is the largest eddy scale, H is the nebula gas vertical scale height, alpha the turbulent viscosity parameter, and eta is the Kolmogorov or smallest scale in turbulence (typically about 1km), with eddy turnover time t(sub eta). In the nebula, Re is far larger than any numerical simulation can handle, so some physical arguments are needed to extend the results of numerical simulations to nebula conditions. In this paper, we report new physics to be incorporated into our statistical models.

mm-size particles↗

Fundamental limitations on photoisomerization from thermodynamic resource theories

Small, out-of-equilibrium, and quantum systems defy simple thermodynamic expressions. Such systems are exemplified by molecular switches, which exchange heat with a bath. These molecules can photoisomerize, or change conformation, or switch, on absorbing light. The photoisomerization probability depends on kinetic details that couple the molecule’s energetics to its dissipation. Therefore, a simple, general, thermodynamic-style bound on the photoisomerization probability seems out of reach. We derive such a bound using a resource theory. The resource-theory framework is a set of mathematical tools, developed in quantum information theory, used to generalize thermodynamics to small and quantum settings. From this toolkit has been derived a generalization of the second law, the thermomajorization preorder. We use thermomajorization to upper-bound the photoisomerization probability. Then, we compare the bound with an equilibrium prediction and with a Lindbladian model. We identify a realistic parameter regime in which the Lindbladian evolution saturates the thermomajorization bound. We also quantify the energy coherence in the electronic degree of freedom, and we argue that this coherence cannot promote photoisomerization. Overall, this work illustrates how quantum-information-theoretic thermodynamics can elucidate complex quantum processes in nature, experiments, and synthetics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Tracking atomic structure evolution during directed electron beam induced Si-atom motion in graphene via deep machine learning

Using electron beam manipulation, we enable deterministic motion of individual Si atoms in graphene along predefined trajectories. Structural evolution during the dopant motion was explored, providing information on changes of the Si atom neighborhood during atomic motion and providing statistical information of possible defect configurations. The combination of a Gaussian mixture model and principal component analysis applied to the deep learning-processed experimental data allowed disentangling of the atomic distortions for two different graphene sublattices. This approach demonstrates the potential of e-beam manipulation to create defect libraries of multiple realizations of the same defect and explore the potential of symmetry breaking physics. The rapid image analytics enabled via a deep learning network further empowers instrumentation for e-beam controlled atom-by-atom fabrication. Here, the analysis described in the paper can be reproduced via an interactive Jupyter notebook at https://git.io/JJ3Bx.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Statistical tests of additional plate boundaries from plate motion inversions

The application of the F-ratio test, a standard statistical technique, to the results of relative plate motion inversions has been investigated. The method tests whether the improvement in fit of the model to the data resulting from the addition of another plate to the model is greater than that expected purely by chance. This approach appears to be useful in determining whether additional plate boundaries are justified. Previous results have been confirmed favoring separate North American and South American plates with a boundary located beween 30 N and the equator. Using Chase's global relative motion data, it is shown that in addition to separate West African and Somalian plates, separate West Indian and Australian plates, with a best-fitting boundary between 70 E and 90 E, can be resolved. These results are generally consistent with the observation that the Indian plate's internal deformation extends somewhat westward of the Ninetyeast Ridge. The relative motion pole is similar to Minster and Jordan's and predicts the NW-SE compression observed in earthquake mechanisms near the Ninetyeast Ridge.

Stein, S.↗