Search NASA⌕ Search

SEARCH · Search NASA

Results for “random fields”

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,387 records · Page 77

Vortex simulation of the pressure field of a jet

Fluctuations of the pressure field of a jet are simulated numerically by a flow model consisting of axisymmetric vortex rings with viscous cores submerged in an inviscid uniform stream. Vortex shedding time intervals, randomly created to imitate the time-history characteristics of the pressure signals of a jet, are generated based on a probability distribution of the intervals between successive pressure peaks obtained from experiments. It is found that, up to five diameters downstream of the jet exit, the characteristics of the pressure fluctuations and the most probable time intervals between experimental and numerical results show good qualitative agreements. The role played by the axisymmetric vortex model in pressure field as well as extensions of the model is also discussed.

Fung, Y. T.↗

Effects of internal swelling on residual elasticity of a quasi-brittle material through a composite sphere model

This work describes the development of a micromechanical-based constitutive model accounting for the effect of internal expansion on the residual elasticity of a Hashin composite material. This material is made of spherical inclusions that are subjected to gradual swelling within a quasi-brittle matrix. The main focus of this work is to describe and analyze the material mechanical response, with an additional focus on the internal swelling’s effect on the stress–strain response and residual elasticity. Microstructural features and parameters of major importance for the mechanical responses were identified. The innovative characteristics of the proposed approach are summarized as follows: (1) a full determination of the physics of a complete-damage problem throughout the whole process of inclusion swelling with upscaling techniques, which transfers the microcrack-related properties from the lower scale to upper scale; and (2) an evolution of the mechanical fields and the corresponding residual elasticity for various inclusion swelling levels. Concerning the matrix–inclusion composite, it was hypothesized that only the matrix was susceptible to cracking, with varied degrees of damage, whereas the inclusions behave elastically and the elastic modulus of the expanding inclusions remains constant. It is to emphasize that the gradual swelling of inclusions is modeled by an increasing strain in the current micromechanical-based constitutive model, and the microcracks are represented by a set of randomly oriented penny-shaped microcracks with identical radii (namely, closed cracks). The main contribution of the current research is to establish the exact mathematical solutions for the mechanical fields (stress, strain) caused by the swelling of the inclusions (mechanical loading) and derive the effective residual elasticity of a composite (structural response) subject to internal expansion and quasi-brittle damage. Based on the assumption of closed cracks that could be extended to open cracks in the upcoming work, the results proposed by the present paper help to understand the non-linear mechanical behavior of quasi-brittle materials subject to microcracking and provide a theoretical framework to be used as academical benchmark for numerical simulations.

36 MATERIALS SCIENCE↗

Overlap Properties of Clouds Generated by a Cloud Resolving Model

In order for General Circulation Models (GCMs), one of our most important tools to predict future climate, to correctly describe the propagation of solar and thermal radiation through the cloudy atmosphere a realistic description of the vertical distribution of cloud amount is needed. Actually, one needs not only the cloud amounts at different levels of the atmosphere, but also how these cloud amounts are related, in other words, how they overlap. Currently GCMs make some idealized assumptions about cloud overlap, for example that contiguous cloud layers overlap maximally and non-contiguous cloud layers overlap in a random fashion. Since there are difficulties in obtaining the vertical profile of cloud amount from observations, the realism of the overlap assumptions made in GCMs has not been yet rigorously investigated. Recently however, cloud observations from a relatively new type of ground radar have been used to examine the vertical distribution of cloudiness. These observations suggest that the GCM overlap assumptions are dubious. Our study uses cloud fields from sophisticated models dedicated to simulate cloud formation, maintenance, and dissipation called Cloud Resolving Models . These models are generally considered capable of producing realistic three-dimensional representation of cloudiness. Using numerous cloud fields produced by such a CRM we show that the degree of overlap between cloud layers is a function of their separation distance, and is in general described by a combination of the maximum and random overlap assumption, with random overlap dominating as separation distances increase. We show that it is possible to parameterize this behavior in a way that can eventually be incorporated in GCMs. Our results seem to have a significant resemblance to the results from the radar observations despite the completely different nature of the datasets. This consistency is encouraging and will promote development of new radiative transfer codes that will estimate the radiation effects of multi-layer cloud fields more accurately.

Oreopoulos, L.↗

A 3d disordered superconformal fixed point

We initiate the study of a three dimensional disordered supersymmetric field theory. Specifically, we consider a $\mathcal{N}$= 2 large N Wess-Zumino like model with cubic superpotential involving couplings drawn from a Gaussian random ensemble. Taking inspiration from analyses of lower dimensional SYK like models we demonstrate that the theory flows to a strongly coupled superconformal fixed point in the infra-red. In particular, we obtain leading large N spectral data and operator product coefficients at the critical point. Moreover, the analytic control accorded by the model allows us to compare our results against those derived in the conformal bootstrap program and demonstrate consistency with general expectations.

1/N expansion↗

An application of survival analysis to disruption prediction via Random Forests</strong

We report one of the most pressing challenges facing the fusion community is adequately mitigating or, even better, avoiding disruptions of tokamak plasmas. However, before this can be done, disruptions must first be predicted with sufficient warning time to actuate a response. The established field of survival analysis provides a convenient statistical framework for time-to-event (i.e. time-to-disruption) studies. This paper demonstrates the integration of an existing disruption prediction machine learning algorithm with the Kaplan–Meier estimator of survival probability. Specifically discussed are the implied warning times from binary classification of disruption databases and the interpretation of output signals from Random Forest algorithms trained and tested on these databases. This survival analysis approach is applied to both smooth and noisy test data to highlight important features of the survival and hazard functions. In addition, this method is applied to three Alcator C-Mod plasma discharges and compared to a threshold-based scheme for triggering alarms. In one case, both techniques successfully predict the disruption; although, in another, neither warns of the impending disruption with enough time to mitigate. For the final discharge, the survival analysis approach could avoid the false alarm triggered by the threshold method. Limitations of this analysis and opportunities for future work are also presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Non-Wiener Dynamics of Open Systems for Nonzero Number Density of Surrounding Photons

We propose a collective decay model for a localized ensemble of identical atoms in a vacuum broadband electromagnetic field with a nonzero photon number density considering second-order terms in the atom–field interaction constant. The model is based on the stochastic differential equation for the evolution operator for the atomic ensemble and surroundings, which is controlled by the main quantum random processes and independent classical Wiener processes. It is shown that in some cases, nonzero photon number density does not influence the suppression of collective emission of atoms, and its existence is manifested as an additional mechanism of redistribution of atoms over collective sublevels.

74 ATOMIC AND MOLECULAR PHYSICS↗

Compressible pairwise interaction extended point-particle model for force prediction of shock-particle bed interaction

We propose a pairwise influence framework for the complex unsteady compressible particle-laden flow problem by accounting for the scattered hydrodynamic waves emitting from neighboring particles in a Euler-Lagrange simulation. It has been observed from particle-resolved (PR) simulations of randomly dispersed particle beds under a loading shock that the compressible pseudoturbulence dominates the flow system even after the primary shock has passed, which causes fluctuations observed in the forces experienced by the particle. Moreover, the fact that each particle exists in the vicinity of a random arrangement of other particles modifies the time history of the drag force experienced by each particle during and after the passage of the shock. First, the scattering flow field due to an incoming shock interacting with a single sphere is constructed using an analysis of the flow in the acoustic limit. Then we examine the validity of the compressible Maxey-Riley-Gatignol force model by comparing the force prediction against a PR simulation of two interacting particles for various particle arrangements and incoming shock strength. Subsequently, the neighboring influences are stored as a library of maps that can be used readily in the calculation of the perturbation force. Lastly, the pairwise interaction assumption is evaluated by comparing the force predicted with the model with PR simulations of a randomly packed particle bed of 10% volume fraction for both water and air as the fluid medium for an incoming shock Mach number 1.22. With a considerably lower cost for the implementation of the model compared to PR simulations, it is verified that the model is reasonably accurate in pinpointing particles whose peak force is significantly larger or smaller than the mean drag but also to capture the prolonged fluctuations after the initial shock.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Determination of permeability for hydrocarbon release due to excavation-induced stress redistribution in rock salt

Due to a stress redistribution after the excavation of an underground tunnel for radioactive waste disposal, an Ed/DZ (excavation disturbed/damaged zone) will be generated in the near field of the opening, resulting in significant changes in the hydraulic and mechanical properties of the rock mass in the zone. Initially more or less randomly distributed hydrocarbons at grain boundaries in rock salt, which sometimes can only be observed with ultraviolet light, can then be mobilised and migrate at a potentially significant rate towards the excavation. Within the international cooperative project DECOVALEX 2019, the migration mechanism of such fluid inclusions in rock salt is being studied intensively. Specifically, a multi-scale modelling strategy has been developed. A macroscale coupled hydro-mechanical modelling of an underground excavation was performed to determine hydraulic and time-dependent deviatoric stress conditions, by taking into account the rock salt creep behaviour. Under the obtained macro-scale constraints, micro-scale modelling of a pathway dilation along halite grain boundaries was performed using different model strategies: a) coupled hydromechanical modelling with a consideration of hydraulic pressure-induced dilatant deformation, b) nonlinear dynamic model taking account of fluid migration, stress-dependent grain boundary wetting and shear-induced dilatancy of salt, and c) phase-field modelling of flow pathway propagation. The permeability increase resulting from the pathway dilation is estimated to be as high as two orders of magnitude. Based on the permeability determined, a series of pressure buildups measured from a borehole with a high hydrocarbon release rate, a total of 430 build-ups within a monitoring time of 938 days, can be simulated with a macro-scale compressible flow model accounting for different zones around the opening.

58 GEOSCIENCES↗

Parameter estimation for X-ray scattering analysis with Hamiltonian Markov Chain Monte Carlo

Bayesian-inference-based approaches, in particular the random-walk Markov Chain Monte Carlo (MCMC) method, have received much attention recently for X-ray scattering analysis. Hamiltonian MCMC, a state-of-the-art development in the field of MCMC, has become popular in recent years. It utilizes Hamiltonian dynamics for indirect but much more efficient drawings of the model parameters. We described the principle of the Hamiltonian MCMC for inversion problems in X-ray scattering analysis by estimating high-dimensional models for several motivating scenarios in small-angle X-ray scattering, reflectivity, and X-ray fluorescence holography. Hamiltonian MCMC with appropriate preconditioning can deliver superior performance over the random-walk MCMC, and thus can be used as an efficient tool for the statistical analysis of the parameter distributions, as well as model predictions and confidence analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Strong spin-orbit interaction and quadratic temperature dependence of electron-phonon scattering in disorder V{sub 75}X{sub 25} (X = Pd, Al) alloys at low temperature

The results of a comprehensive study of weak electron localization (WEL) and electron-electron interaction (EEI) effects in disordered V{sub 75}X{sub 25} (X = Pd, Al) alloys has been reported. The resistivity in absence of magnetic field shows a minimum at temperature T = T{sub m} and follows T{sup 1/2} law within the temperature range 5 K ≤ T ≤ T{sub m}, which suggests predominant EEI effect. Magnetoresistivity is positive due to strong spin-orbit interaction. The dephasing scattering time is dominated by the electron-phonon scattering. The electron-phonon scattering rate shows quadratic temperature dependence behavior, which is explained by the theory of incomplete dragging at the random scattering potential by phonons. The zero temperature scattering time strongly depends on the disorder and its magnitude decreases with increasing disorder.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Moisture availability mediates the relationship between terrestrial gross primary production and solar-induced chlorophyll fluorescence: Insights from global-scale variations

Effective use of solar-induced chlorophyll fluorescence (SIF) to estimate and monitor gross primary production (GPP) in terrestrial ecosystems requires a comprehensive understanding and quantification of the relationship between SIF and GPP. To date, this understanding is incomplete and somewhat controversial in the literature. Here we derived the GPP/SIF ratio from multiple data sources as a diagnostic metric to explore its global-scale patterns of spatial variation and potential climatic dependence. We found that the growing season GPP/SIF ratio varied substantially across global land surfaces, with the highest ratios consistently found in boreal regions. Spatial variation in GPP/SIF was strongly modulated by climate variables. The most striking pattern was a consistent decrease in GPP/SIF from cold-and-wet climates to hot-and-dry climates. We propose that the reduction in GPP/SIF with decreasing moisture availability may be related to stomatal responses to aridity. Furthermore, we show that GPP/SIF can be empirically modeled from climate variables using a machine learning (random forest) framework, which can improve the modeling of ecosystem production and quantify its uncertainty in global terrestrial biosphere models. Finally, our results point to the need for targeted field and experimental studies to better understand the patterns observed and to improve the modeling of the relationship between SIF and GPP over broad scales.

59 BASIC BIOLOGICAL SCIENCES↗

The Mechanism for the Energy Buildup Driving Solar Eruptive Events

The underlying origin of solar eruptive events (SEEs), ranging from giant coronal mass ejections to small coronalhole jets, is that the lowest-lying magnetic flux in the Sun's corona undergoes continual buildup of stress and free energy. This magnetic stress has long been observed as the phenomenon of "filament channels:" strongly sheared magnetic field localized around photospheric polarity inversion lines. However, the mechanism for the stress buildup-formation of filament channels-is still debated. We present magnetohydrodynamic simulations of a coronal volume that is driven by transient, cellular boundary flows designed to model the processes by which the photosphere drives the corona. The key feature of our simulations is that they accurately preserve magnetic helicity, the topological quantity that is conserved even in the presence of ubiquitous magnetic reconnection. Although small-scale random stress is injected everywhere at the photosphere, driving stochastic reconnection throughout the corona, the net result of the magnetic evolution is a coherent shearing of the lowest-lying field lines. This highly counterintuitive result-magnetic stress builds up locally rather than spreading out to attain a minimum energy state-explains the formation of filament channels and is the fundamental mechanism underlying SEEs. Furthermore, this process is likely to be relevant to other astrophysical and laboratory plasmas.

Knizhnik, K. J.↗

Astrometry with Hubble Space Telescope Fine Guidance Sensor number 3: Position-mode stability and precision

We report results from a test exploring the long- and short-term astrometric stability of Hubble Space Telescope Fine Guidance Sensor (FGS) #3. A test field was observed 40 times over 522 days to determine the precision and accuracy of FGS astrometry and to measure the character and magnitude of possible secular scale changes. We examine the astrometric data and the associated guide-star data to determine random errors. These data are also explored to find sources of systematic error. After correcting for some systematic effects we obtain a precision of 0.002 arcsec (2 mas) per observation (RSS of x and y). This is relative astrometry within a central 2.5 arcmin FGS field of view for any orientation. We find that the scale varies over time and confirm the sense of the trend with independent data. From the 40 observation sets we produce a catalog of an astrometry test field containing eight stars whose relative positions are known to an average 0.7 and 0.9 mas in x and y. One reference star has a relative parallax of 3.1 plus or minus 0.5 mas. Finally, we report that eleven observation sets acquired over 387 days produce parallaxes and relative positions with 1-mas precision.

Benedict, G. F.↗

On the replica structure of Sachdev-Ye-Kitaev model

We investigate existence of replica off-diagonal solutions in the field-theoretical description of Sachdev-Ye-Kitaev model. To this end we evaluate a set of local and non-local dynamic correlation functions in the long time limit. We argue that the structure of the soft-mode Schwarzian action is qualitatively different in replica-diagonal vs. replica-off-diagonal scenarios, leading to distinct long-time predictions for the correlation functions. We then evaluate the corresponding correlation functions numerically and compare the simulations with analytical predictions of replica-diagonal and replica-off-diagonal calculations. We conclude that all our numerical results are in a quantitative agreement with the theory based on the replica-diagonal saddle point plus Schwarzian and massive Gaussian fluctuations (the latter do contain replica off-diagonal components). This seems to exclude any contributions from replica-off-diagonal saddle points, at least on the time scales shorter than the inverse many-body level spacing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Out-of-time-order correlators and Lyapunov exponents in sparse SYK

We use a combination of analytical and numerical methods to study out-of-time order correlators (OTOCs) in the sparse Sachdev-Ye-Kitaev (SYK) model. We find that at a given order of N, the standard result for the q-local, all-to-all SYK, obtained through the sum over ladder diagrams, is corrected by a series in the sparsity parameter, k. We present an algorithm to sum the diagrams at any given order of 1/(kq) n . We also study OTOCs numerically as a function of the sparsity parameter and determine the Lyapunov exponent. We find that numerical stability when extracting the Lyapunov exponent requires averaging over a massive number of realizations. This trade-off between the efficiency of the sparse model and consistent behavior at finite N becomes more significant for larger values of N.

2D gravity↗

Emergence of high piezoelectricity from competing local polar order-disorder in relaxor ferroelectrics

Relaxor ferroelectrics are known for outstanding piezoelectric properties, finding a broad range of applications in advanced electromechanical devices. Decoding the origins of the enhanced properties, however, have long been complicated by the heterogeneous local structures. Here, we employ the advanced big-box refinement method by fitting neutron-, X-ray-based total scattering, and X-ray absorption spectrum simultaneously, to extract local atomic polar displacements and construct 3D polar configurations in the classical relaxor ferroelectric Pb(Mg 1/3 Nb 2/3 )O 3 –PbTiO 3 . Our results demonstrate that prevailing order-disorder character accompanied by the continuous rotation of local polar displacements commands the composition-driven global structure evolution. The omnidirectional local polar disordering appears as an indication of macroscopic relaxor characteristics. Combined with phase-field simulations, it demonstrates that the competing local polar order-disorder between different states with balanced local polar length and direction randomness leads to a flattening free-energy profile over a wide polar length, thus giving rise to high piezoelectricity. Our work clarifies that the critical structural feature required for high piezoelectricity is the competition states of local polar rather than relaxor.

42 ENGINEERING↗