Search NASA⌕ Search

SEARCH · Search NASA

Results for “Two point model”

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,063 records · Page 59

Multishaker modal testing

Procedures for improving the modal modeling of structures using test data and to determine appropriate analytical models based on substructure experimental data were explored. Two related research topics were considered in modal modeling: using several independently acquired columns of frequency response data, and modal modeling using simultaneous multi-point excitation. In component mode synthesis modeling, the emphasis is on determining the best way to employ complex modes and residuals.

Craig, R. R., Jr.↗

Interfaces for high-speed fiber-optic links - Analysis and experiment

An analysis of directly and externally modulated fiberoptic links is presented. The theoretical analysis is based on the signal flow diagram of the interface circuits to the laser diode, Mach-Zehnder electro-optic modulator, and p-i-n photodiode. System parameters such as gain, noise figure, two-tone intermodulation distortion, and spurious free and compression dynamic range are expressed as a function of frequency and operating point of the laser and external modulator. Two directly and externally modulated fiber-optic links were designed and fabricated to verify the analytical models. The direct modulation FO link was developed at the Ku-band (11.6-12.4 GHz), whereas the external modulation link was at L-band (875-925 MHz). Spurious-free dynamic ranges of 95.8 dB-Hz exp 2/3 and 113 dB-Hz exp 2/3, respectively, were achieved. The predictions based on the analytical models match the measured results.

Daryoush, Afshin S.↗

Results of a zonally truncated three-dimensional model of the Venus middle atmosphere

Although the equatorial rotational speed of the solid surface of Venus is only 4 m s(exp-1), the atmospheric rotational speed reaches a maximum of approximately 100 m s(exp-1) near the equatorial cloud top level (65 to 70 km). This phenomenon, known as superrotation, is the central dynamical problem of the Venus atmosphere. We report here the results of numerical simulations aimed at clarifying the mechanism for maintaining the equatorial cloud top rotation. Maintenance of an equatorial rotational speed maximum above the surface requires waves or eddies that systematically transport angular momentum against its zonal mean gradient. The zonally symmetric Hadley circulation is driven thermally and acts to reduce the rotational speed at the equatorial cloud top level; thus wave or eddy transport must counter this tendency as well as friction. Planetary waves arising from horizontal shear instability of the zonal flow (barotropic instability) could maintain the equatorial rotation by transporting angular momentum horizontally from midlatitudes toward the equator. Alternatively, vertically propagating waves could provide the required momentum source. The relative motion between the rotating atmosphere and the pattern of solar heating, which as a maximum where solar radiation is absorbed near the cloud tops, drives diurnal and semidiurnal thermal tides that propagate vertically away from the cloud top level. The effect of this wave propagation is to transport momentum toward the cloud top level at low latitudes and accelerate the mean zonal flow there. We employ a semispectral primitive equation model with a zonal mean flow and zonal wavenumbers 1 and 2. These waves correspond to the diurnal and semidiurnal tides, but they can also be excited by barotropic or baroclinic instability. Waves of higher wavenumbers and interactions between the waves are neglected. Symmetry about the equator is assumed, so the model applies to one hemisphere and covers the altitude range 30 to 110 km. Horizontal resolution is 1.5 deg latitude, and vertical resolution is 1.5 km. Solar and thermal infrared heating, based on Venus observations and calculations drive the model flow. Dissipation is accomplished mainly by Rayleigh friction, chosen to produce strong dissipation above 85 km in order to absorb upward propagating waves and limit extreme flow velocities there, yet to give very weak Rayleigh friction below 70 km; results in the cloud layer do not appear to be sensitive to the Rayleigh friction. The model also has weak vertical diffusion, and very weak horizontal diffusion, which has a smoothing effect on the flow only at the two grid points nearest the pole.

Newman, M.↗

Research activities at the Center for Modeling of Turbulence and Transition

The main research activities at the Center for Modeling of Turbulence and Transition (CMOTT) are described. The research objective of CMOTT is to improve and/or develop turbulence and transition models for propulsion systems. The flows of interest in propulsion systems can be both compressible and incompressible, three dimensional, bounded by complex wall geometries, chemically reacting, and involve 'bypass' transition. The most relevant turbulence and transition models for the above flows are one- and two-equation eddy viscosity models, Reynolds stress algebraic- and transport-equation models, pdf models, and multiple-scale models. All these models are classified as one-point closure schemes since only one-point (in time and space) turbulent correlations, such as second moments (Reynolds stresses and turbulent heat fluxes) and third moments, are involved. In computational fluid dynamics, all turbulent quantities are one-point correlations. Therefore, the study of one-point turbulent closure schemes is the focus of our turbulence research. However, other research, such as the renormalization group theory, the direct interaction approximation method, and numerical simulations are also pursued to support the development of turbulence modeling.

Shih, Tsan-Hsing↗

Phase Sensitive X-Ray Diffraction Imaging Study of Protein Crystals

The study of defects and growth of protein crystals is of importance in providing a fundamental understanding of this important category of systems and the rationale for crystallization of better ordered crystals for structural determination and drug design. Yet, as a result of the extremely weak scattering power of x-rays in protein and other biological macromolecular crystals, the extinction lengths for those crystals are extremely large and, roughly speaking, of the order of millimeters on average compared to the scale of micrometers for most small molecular crystals. This has significant implication for x-ray diffraction and imaging study of protein crystals, and presents an interesting challenge to currently available x-ray analytical techniques. We proposed that coherence-based phase sensitive x-ray diffraction imaging could provide a way to augment defect contrast in x-ray diffraction images of weakly diffracting biological macromolecular crystals. I shall examine the principles and ideas behind this approach and compare it to other available x-ray topography and diffraction methods. I shall then present some recent experimental results in two model protein systems-cubic apofemtin and tetragonal lysozyme crystals to demonstrate the capability of the coherence-based imaging method in mapping point defects, dislocations, and the degree of perfection of biological macromolecular crystals with extreme sensitivity. While further work is under way, it is intended to show that the observed new features have yielded important information on protein crystal perfection and nucleation and growth mechanism otherwise unobtainable.

Hu, Z. W.↗

Planetary Crater Detection and Registration Using Marked Point Processes, Graph Cut Algorithms, and Wavelet Transforms

This paper addresses the problem of semi-automatic image registration on planetary images. A joint feature-based and area-based approach is proposed. Firstly, the most relevant craters are extracted from the two images to register, and then, registration is performed in two steps. The first step matches the craters extracted from the images based on a generalized Hausdorff distance. In the second step, the mutual information between the two images is maximized to achieve high registration accuracy. Craters are detected by a stochastic-geometry approach based on a marked point process model and of a multiple-birth-and-cut energy minimization algorithm. The experimental validation is carried out with 13 images for the crater extraction stage, and with 20 semi-synthetic pairs of images with ground truth and several images extracted from actual multi-temporal lunar scenes for the registration phase.

Image Processing; Pattern Recognition↗

Scaling Input Distributions for Probabilistic Models - 19472

It is important to address scaling issues when building models to support risk-based decisions for radioactive waste disposal sites. Currently, probabilistic risk assessment and probabilistic performance assessment models are run for multiple realizations where each realization uses one value from each input distribution over the entire modeled time. The input distributions represent the current state of knowledge in the mean of the input variable over the spatial region of interest and over the time period of the model run. This paper focuses on the time factor in input distribution scaling because it presents more serious conceptual and distributional challenges than the spatial factor. One of the challenges in temporal scaling is the effect of non-linearity on propagation of distribution uncertainty through to uncertainty in results. Radioactive waste disposal models tend to have both linear and non-linear components in their models. Dose assessment calculations are typically represented by linear models, while contaminant transport models have linear and non-linear components. In principle, temporal scaling of input parameter distributions for a linear model is straightforward, but scaling of input distributions for non-linear models is not. This paper uses examples to investigate the effects of input distribution temporal scaling on simple linear and non-linear models. Different scaling scenarios were investigated by comparing three analysis tools: analytical mathematical representations, statistical software programming, and GoldSim Monte Carlo simulation software. R statistical software is an open source statistical programming language often used for statistical analyses and graphical representation of data. GoldSim is a robust probabilistic dynamic modeling platform used in performance assessment modeling of radioactive waste sites. A simple linear model, a quadratic model and a multiplicative model were the cases explored for each of these tools. For the three analysis tools considered and with each of the scaling cases, 1,000 realizations were implemented from distributions of the random variables of interest. This is a sufficient number of realizations for simple models to capture the model space. An annual temporal resolution of a 100-year simulation period was chosen for all models and analysis tools. The effects of temporal scaling were demonstrated by investigating two different options for time scaling. Option 1 was set up so that the input distributions represented the expected value on an annual time scale, and a new value from that distribution was chosen every year (at every time step which equaled a year) for 100 years. In Option 2, the input distributions were scaled to represent the expected mean value of the input variable over the 100 years of simulation, and one value from that distribution was chosen for the entire elapsed time. Results show agreement across the three analysis tools. As expected, the linear models have identical results for the two different scaling Options, and the non-linear models diverge in the results for different time scale implementation. The implications of this work point to the need for careful consideration of input distribution temporal scaling. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Modelling of Collision Induced Absorption Spectra Of H2-H2 Pairs for the Planetary Atmospheres Structure: The Second Overtone Band

The main objective of the proposal was to model the collision induced, second overtone band of gaseous hydrogen at low temperatures. The aim of this work is to assist planetary scientists in their investigation of planetary atmospheres, mainly those of Uranus and Neptune. The recently completed extended database of collision induced dipole moments of hydrogen pairs allowed us, for the first time, to obtain dipole moment matrix elements responsible for the roto-vibrational collision induced absorption spectra of H2-H2 in the second overtone band. Despite our numerous attempts to publish those data, the enormous volume of the database did not allow us to do this. Instead, we deposited the data on a www site. The final part of this work has been partially supported by NASA, Division for Planetary Atmospheres. In order to use our new data for modelling purpose, we first needed to test how well we can reproduce the existing experimental data from theory, when using our new input data. Two papers resulted from this work. The obtained agreement between theoretical results and the measurements appeared to be within 10-30%. The obviously poorer agreement than observed for the first H2 overtone, the fundamental, and the rototranslational bands can be attributed to the fact that dipole moments responsible for the second overtone are much weaker, therefore susceptible to larger numerical uncertainties. At the same time, the intensity of the second overtone band is much weaker and therefore it is much harder to be measured accurately in the laboratory. We need to point out that until now, no dependable model of the 2nd overtone band was available for modelling of the planetary atmospheres. The only one, often referred to in previous works on Uranian and Neptune's atmospheres, uses only one lineshape, with one (or two) parameter(s) deduced at the effective temperature of Uranus (by fitting the planetary observation). After that, the parameter(s) was(were) made temperature dependent according to some very simple relation. Summarizing, no reliable temperature-dependent model has been available yet. Our approach was a bit different from similar attempts done earlier, on account of the poorer agreement of theory with experiment. We needed to resort to some semi-empirical procedure. While we were in a favourable position to be able to rely on the physical input data, these, apparently, did not supply the most dependable predictions (simply because the results did not agree well enough with experimental data). On the other hand, the relative deviations between the theory and experiment were comparable at 77 and at 298 K. That fact indicated that theory is capable of predicting the temperature dependence of the absorption spectra well. We have thus chosen the "middle way". We have fitted the existing measurements with many 3- parameter lineshapes, in order to achieve the closest fit.

Borysow, Aleksandra↗

Fission gas diffusion and release for Cr 2 O 3 -doped UO 2 : From the atomic to the engineering scale

Here, the anticipated benefits of large grains in Cr 2 O 3 -doped UO 2 pellets include improved mechanical and fission gas retention properties. To support the assessment of fission gas release (FGR) from doped pellets, the impact of doping on fission gas diffusivity for in-reactor conditions must be understood. In this work, we tackle this issue by informing the fission gas model within the BISON fuel performance code using material models developed at the atomic scale. The investigation of intra-granular fission gas diffusivity in Cr 2 O 3 -doped UO 2 is carried out by adapting a cluster dynamics model that, accounting for UO 2 thermochemistry, is capable of describing Xe diffusion under irradiation in undoped UO 2 as the starting point. Using a thermodynamic analysis, it is shown that in stoichiometric UO 2 with additions of Cr 2 O 3 the oxygen potential is defined by the Cr-Cr 2 O 3 two-phase equilibrium. Using the cluster dynamics model, the predicted Xe diffusivity in doped UO 2 was significantly increased in both the intrinsic and irradiation-enhanced regimes compared to undoped UO 2 as a result of higher concentrations of uranium and oxygen vacancies, respectively. This is a consequence of the more oxidizing conditions at high temperature, and more reducing conditions at low temperature, as a result of doping. Arrhenius functions have been fitted to the cluster dynamics results to enable implementation of the new diffusivities in the BISON fission gas behavior model. BISON simulations were carried out, showing the competing effects of the enlarged grains and the new fission gas diffusivity model, which act to suppress and enhance fission gas release, respectively. The new physics-informed model was validated against in-reactor experimental measurements under normal operation. Additionally, benchmarking was carried out for power ramp conditions. The predicted fission gas release agreed well with the experimental data, showing noticeable improvements over the standard UO 2 model.

36 MATERIALS SCIENCE↗

Multiscale spatial analysis of fracture nodes in two dimensions

Spatial arrangement of fractures as a function of scale is an important component of fracture quantification for inferential and predictive modeling. Available methods that analyze fracture spatial arrangement are based on one-dimensional spacing data; therefore, they are limited to semi-parallel fractures. Such methods cannot be applied to fracture networks in higher dimensions, particularly when fractures have different orientations. Here to characterize fracture arrangements in two dimensions, we propose using Ripley’s K-function, as a method of point pattern analysis, to quantify spatial arrangement of fracture nodes. Fracture nodes, such as barycenters, intersection points, and tips, are point-based representations of fracture locations and connectivity within the fracture network. We introduce formulations for isotropic as well as directional analyses of spatial arrangement. In addition, we derive formulations for edge correction in circular and rectangular study domains. Finally, we demonstrate applications of Ripley’s K-function on two natural fracture datasets. Our proposed method supports quantification and characterization of fracture spatial arrangements that allow practitioners to build representative models of fractures in the subsurface.

02 PETROLEUM↗

Instantons in chiral magnets

In this work, we exhaustively construct instanton solutions and elucidate their properties in one-dimensional antiferromagnetic chiral magnets based on the $\textit{O}(3)$ nonlinear sigma model description of spin chains with the Dzyaloshinskii-Moriya (DM) interaction. By introducing an easy-axis potential and a staggered magnetic field, we obtain a phase diagram consisting of ground-state phases with two points (or one point) in the easy-axis dominant cases, a helical modulation at a fixed latitude of the sphere, and a tricritical point allowing helical modulations at an arbitrary latitude. We find that instantons (or skyrmions in two-dimensional Euclidean space) appear as composite solitons in different fashions in these phases: temporal domain walls or wall-antiwall pairs (bions) in the easy-axis dominant cases, dislocations (or phase slips called merons) with fractional instanton numbers in the helical state, and isolated instantons and calorons living on the top of the helical modulation at the tricritical point. We also show that the models with DM interaction and an easy-plane potential can be mapped into those without them, providing a useful tool to investigate the model with the DM interaction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Chronology of surface units on the icy satellites of Saturn

Because the impact flux history in the Saturnian system has been similar in its general form to that for the inner solar system and Jupiter, it is possible to develop an absolute calibration model for the relative ages derived from the density of impact craters. This has been done by identifying the density at which the crater populations change characteristics from that during post-accretion heavy bombardment to a later population of low but continuous flux. The transition between these two populations occurred at a known time throughout the solar system, 3.9 by, and provides an absolute calibration point on which to construct a calibration curve. Results based on calibrating relative age data to this absolute flux model suggest that most surfaces are ancient, greater than 3.9 by (or within several 100 million years). Only Enceladus and Dione have had prolonged activity that has resulted in resurfacing. The activity on Enceladus was probably driven by tidal forces, whereas the mechanism for driving resurfacing on Dione has yet to be identified.

Boyce, J. M.↗

A dimensionality reducing model for distributed filtering.

An approach is made to filtering of two-dimensional steady-state problems based on the notion of nearest neighbor interaction, i.e., at a given point in the spatial grid, the value at the variable of interest can be assumed to depend only on the values at adjacent grid points. It is shown that for linear steady-state problems significant dimensionality reductions can be accomplished. It was possible to achieve the desired results using a small amount of computer time and without getting into stability difficulties.

Angel, E.↗

Dynamics modeling of multibody flexible systems with traveling multi-point joints

This paper presents a general methodology for modeling the dynamics of a flexible structure which moves on the surface of another flexible structure. The relative motion of the two bodies need not follow a straight line, nor need the contact surfaces be planar in their undeformed state. A model reduction approach is also developed, which makes the methodology applicable to large structural systems. The Space Station Freedom Mobile Transporter represents one such system, when attached to the Shuttle Orbiter. Numerical examples are provided.

Messac, Achille↗

Uncovering Non-Fermi-Liquid Behavior in Hund Metals: Conformal Field Theory Analysis of an SU(2)×SU(3) Spin-Orbital Kondo Model

Hund metals have attracted attention in recent years due to their unconventional superconductivity, which supposedly originates from non-Fermi-liquid (NFL) properties of the normal state. When studying Hund metals using dynamical mean-field theory, one arrives at a self-consistent “Hund impurity problem” involving a multiorbital quantum impurity with nonzero Hund coupling interacting with a metallic bath. If its spin and orbital degrees of freedom are screened at different energy scales, T sp < T orb , the intermediate energy window is governed by a novel NFL fixed point, whose nature had not yet been clarified. We resolve this problem by providing an analytical solution of a paradigmatic example of a Hund impurity problem, involving two spin and three orbital degrees of freedom. To this end, we combine a state-of-the-art implementation of the numerical renormalization group, capable of exploiting non-Abelian symmetries, with a generalization of Affleck and Ludwig’s conformal field theory (CFT) approach for multichannel Kondo models. We characterize the NFL fixed point of Hund metals in detail for a Kondo model with an impurity forming an SU ( 2 ) × SU ( 3 ) spin-orbital multiplet, tuned such that the NFL energy window is very wide. The impurity’s spin and orbital susceptibilities then exhibit striking power-law behavior, which we explain using CFT arguments. We find excellent agreement between CFT predictions and numerical renormalization group results. Our main physical conclusion is that the regime of spin-orbital separation, where orbital degrees of freedom have been screened but spin degrees of freedom have not, features anomalously strong local spin fluctuations: the impurity susceptibility increases as χ sp imp ~ ω – γ , with γ > 1 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Implementation of a Parallel Kalman Filter for Stratospheric Chemical Tracer Assimilation

A Kalman filter for the assimilation of long-lived atmospheric chemical constituents has been developed for two-dimensional transport models on isentropic surfaces over the globe. An important attribute of the Kalman filter is that it calculates error covariances of the constituent fields using the tracer dynamics. Consequently, the current Kalman-filter assimilation is a five-dimensional problem (coordinates of two points and time), and it can only be handled on computers with large memory and high floating point speed. In this paper, an implementation of the Kalman filter for distributed-memory, message-passing parallel computers is discussed. Two approaches were studied: an operator decomposition and a covariance decomposition. The latter was found to be more scalable than the former, and it possesses the property that the dynamical model does not need to be parallelized, which is of considerable practical advantage. This code is currently used to assimilate constituent data retrieved by limb sounders on the Upper Atmosphere Research Satellite. Tests of the code examined the variance transport and observability properties. Aspects of the parallel implementation, some timing results, and a brief discussion of the physical results will be presented.

Chang, Lang-Ping↗

Improving the Performance of NEML2 with Modern Graph Compilation Backends

NEML2 vectorizes constitutive-model evaluation for large-scale multiphysics simulation, using PyTorch as its tensor backend so that a batch of material-point updates runs on CPU or GPU through a single implementation. In the two prior reports in this series it was a C++-native library, deployed through TorchScript tracing and just-in-time (JIT) compilation; it has since been rewritten from the ground up into a Python-native library deployed through Ahead-of-Time Inductor (AOTInductor), a modern PyTorch graph-compilation backend. The rewrite is driven by a persistent tension, not a language preference: NEML2 composes constitutive models at runtime from a registry of small, independently-authored pieces, and that flexibility is difficult to reconcile with the compile-time knowledge an efficient GPU kernel needs. This report documents the rewrite and the investment that accompanied it: the AOTInductor export pipeline that turns a Python-authored model into a portable, Python-free compiled artifact loadable from pure C++; the eager and compiled runtimes and the new implicit solver layer built on them; a head-to-head benchmark of legacy JIT against AOTInductor; the physics-model catalog and its worked examples; the developer tooling; and the corresponding overhaul of MOOSE’s NEML2 integration that lets MOOSE consume it. A central objective is to examine whether modern PyTorch graph-compilation backends are effective for MOOSE GPU integration. The benchmark answers directly: AOTInductor outperforms legacy JIT on every GPU scenario measured, by 1.0–4.5×. Modern graph-compilation backends are effective for MOOSE GPU integration, and AOTInductor specifically – not compilation in the abstract – is why.

Hu, Gary (Tianchen) [Argonne National Laboratory (↗

Prediction of an axisymmetric combusting flow

A numerical model for turbulent, recirculating combusting flow is developed and applied to a research combustor. The model is based on the time-averaged Navier-Stokes equations with k-epsilon turbulence closure and the compositional fluctuations at each point are given probabilistically in terms of the mixture fraction. The probability density function is derived from transport equations for its first two moments along with an assumption regarding its shape. The resulting equations are solved using a standard line relaxation algorithm. It is found that the predictions of the model are in good agreement with data from an axisymmetric, bluff-body stabilized research combustor, while the major discrepancies are similar to those found in isothermal flow comparisons. The peculiar features of this flow which contribute to the errors are examined. The agreement between the theory and data deteriorates as the central jet velocity is increased which indicates an enhanced role for unsteady effects.

Correa, S. M.↗