Data-Driven Mori-Zwanzig: Reduced Order Modeling of Sparse Sensor Measurements for Boundary Layer Transition
Explore the source record for details and available documents.
Reduced Models, Data-Driven, Mori-Zwanzig, Hyperso↗
SEARCH · Engineering Papers
Results for “Transition 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.
Simulating neutron stars with a flexible enthalpy-based equation of state parametrization in sp ECTRE
Numerical simulations of neutron star mergers represent an essential step toward interpreting the full complexity of multimessenger observations and constraining the properties of supranuclear matter. Currently, simulations are limited by an array of factors, including computational performance and input physics uncertainties, such as the neutron star equation of state. In this work, we expand the range of nuclear phenomenology efficiently available to simulations by introducing a new analytic parametrization of cold, beta-equilibrated matter that is based on the relativistic enthalpy. Here we show that the new enthalpy parametrization can capture a range of nuclear behavior, including strong phase transitions. We implement the enthalpy parametrization in the spectre code, simulate isolated neutron stars, and compare performance to the commonly used spectral and polytropic parametrizations. We find comparable computational performance for nuclear models that are well represented by either parametrization, such as simple hadronic equations of state. We show that the enthalpy parametrization further allows us to simulate more complicated hadronic models or models with phase transitions that are inaccessible to current parametrizations.
79 ASTRONOMY AND ASTROPHYSICS↗
DEMographic MicrOSimulation (DEMOS) v2.0
DEMOS is an agent-based simulation framework used to evolve population demographic characteristics or lifecycle events such as education, marital status etc. DEMOS modules are designed to capture the interdependencies of short-term and long-term lifecycle events often influential in downstream transportation and land use modeling. An important facet of DEMOS is that it is dynamic, meaning it captures the impact of an agent's demographic characteristics in year 't' on their demographic status in year 't+1'. This has important consequences on medium- and long-term downstream modeling in different domains, such as transportation where decisions such as household vehicle transactions (i.e., buying, selling, or replacing a vehicle) depend current or past household vehicle holdings as well as household transitions. Other downstream modeling that can be enabled by DEMOS include building technology adoption or retrofit decisions, and more. DEMOS modules include: individual age progression year by year, individual mortality events, household restructuring events including marriage match-making; household birth events; individual education participation; individual labor force participation; child leaving household; household in-migration and out-migration for the region; household residence location choice and household mandatory destination choice (work and school). This version of DEMOS is tailored to be integrated with the land use simulation software UrbanSim.
Caicedo, Juan [UrbanSim, Inc.]↗
Signatures of deconfined quark phases in binary neutron star mergers
We investigate the quark deconfinement phase transition in the context of binary neutron star (BNS) mergers. We employ a new finite-temperature composition-dependent equation of state (EOS) with a first order phase transition between hadrons and deconfined quarks to perform numerical relativity simulations of BNS mergers. Here, the softening of the EOS due to the phase transition causes the merger remnants to be more compact and to collapse to a black hole (BH) at earlier times. The phase transition is imprinted on the postmerger gravitational wave (GW) signal duration, amplitude, and peak frequency. However, this imprint is only detectable for binaries with sufficiently long-lived remnants. Moreover, the phase transition does not result in significant deviations from quasi-universal relations for the postmerger GW peak frequency. We also study the impact of the phase transition on dynamical ejecta, remnant accretion disk masses, r-process nucleosynthetic yields and associated electromagnetic (EM) counterparts. While there are differences in the EM counterparts and nucleosynthesis yields between the purely hadronic models and the models with phase transitions, these can be primarily ascribed to the difference in remnant collapse time between the two. An exception is the non-thermal afterglow caused by the interaction of the fastest component of the dynamical ejecta and the interstellar medium, which is systematically boosted in the binaries with phase transition as a consequence of the more violent merger they experience.
79 ASTRONOMY AND ASTROPHYSICS↗
Thermomechanical modeling of the Colorado Plateau-Basin and range transition zone
The Colorado Plateau (CP) basin and range (B & R) boundary is marked by a transition zone on the order of 75 to 150 km in width. As one moves westward across this transition from the CP interior to the B & R there is a variation in the surface topography, surface heat flow, Bouguer gravity, seismicity, and crustal structure. This transition extends eastward into the western CP from the Wastach-Hurricane fault line and is largely coincident with the high plateaus of Utah and the Wasatch Mountains. It has been suggested that this transition zone marks a thermal and tectonic encroachment of the CP by the B & R. A simple two dimensional numerical model of the thermal regime for the transition zone was constructed to test the hypothesis that the observed geophysical signatures across the transition are due to lateral heat conduction from steady state uniform extension within the B & R lithosphere. Surface heat flow, uplift due to flexure from thermal buoyant loading, and regional Bouguer gravity are computed for various extension rates, crustal structures, and compensation depths.
Londe, M. D.↗
Numerical modeling of plasma assisted deflagration to detonation transition in a microscale channel
Here, this work numerically studies the plasma assisted deflagration to detonation transition (DDT) of H 2 /O 2 mixtures in a microscale channel with detailed chemistry and transport. The results show that the DDT onset time is non-monotonically dependent on the discharge pulse number. The DDT is accelerated with small pulse numbers, whereas retarded with large ones. Two different DDT regimes, respectively at a small and large plasma discharge number, via acoustic choking of the burned gas and plasma-enhanced reactivity gradient without acoustic choking, are observed. Without plasma discharge, pronounced pressure and temperature gradients in front of the flame are generated by acoustic compression after the choking of the burned gas, triggering DDT via autoignition. With small plasma pulse numbers, the plasma-generated species enhance the ignition kinetics and lead to an increased reactivity in the boundary layer. After the choking of the burned gas, the plasma-enhanced reactivity advances the sequence of autoignition near the wall, strengthens ignition-shock wave coupling, and accelerates DDT. However, with a large discharge pulse number, a direct autoignition initiating DDT can occur without the acoustic choking of the burned gas due to the strongly accelerated reactivity and elevated temperature. In this case, DDT onset is retarded because the elevated temperature increases sonic velocity and the increased reactivity accelerates fuel oxidation in front of the flame, decelerating the formation of a leading shock and subsequent pressure buildup ahead of the flame. The present modeling reveals that no matter with or without plasma discharge, DDT is initiated by autoignition in thermal, pressure, and reactivity gradient fields via the Zel'dovich gradient mechanism. The acoustic choking of the burned gas may not be the necessary condition of DDT with strong plasma-enhanced reactivity gradient. This work provides an answer to the experimentally observed non-monotonic DDT onset time by plasma, which provides guidance to control DDT in advanced detonation engines and fire safety of hydrogen-fueled catalytic reactors in microchannels by non-equilibrium plasma discharge.
33 ADVANCED PROPULSION SYSTEMS↗
Thermodynamic Control of Activity Patterns in Cytoskeletal Networks
Biological materials, such as the actin cytoskeleton, exhibit remarkable structural adaptability to various external stimuli by consuming different amounts of energy. In this Letter, we use methods from large deviation theory to identify a thermodynamic control principle for structural transitions in a model cytoskeletal network. Specifically, we demonstrate that biasing the dynamics with respect to the work done by nonequilibrium components effectively renormalizes the interaction strength between such components, which can eventually result in a morphological transition. Further, our work demonstrates how a thermodynamic quantity can be used to renormalize effective interactions, which in turn can tune structure in a predictable manner, suggesting a thermodynamic principle for the control of cytoskeletal structure and dynamics.
59 BASIC BIOLOGICAL SCIENCES↗
Effects of mass addition on blunt-body boundary-layer transition and heat transfer
The model bodies tested at Mach number 7.32 were hemispheres, blunt cones, and spherical segments. The mass addition consisted of air ejected through porous forward surfaces of the models. The experimental data consisted of heat transfer measurements from which boundary layer transitions were deduced. The data verified various applicable boundary layer codes in the laminar and transitional flow regimes. Empirical heating rate data correlations were developed for the laminar and turbulent flow regimes.
Kaattari, G. E.↗
Automated theoretical chemical kinetics: Predicting the kinetics for the initial stages of pyrolysis
Large scale implementation of high level computational theoretical chemical kinetics offers the prospect for dramatically improving the fidelity of combustion chemical modeling. To facilitate such efforts, we have developed a suite of codes, collectively referred to as AutoMech, that allow for the automatic prediction of the kinetics for large sets of reactions via ab initio transition-state-theory based master-equation calculations. The primary input is simply the mechanism, a dictionary relating chemically identifiable species descriptors (e.g., SMILES or InChIs) to species labels in the mechanism, and a specification of the electronic structure and transition state theory models to be implemented. Here we illustrate the current utility of AutoMech through a study of the initial stages of pyrolysis for 3 sets of fuels: sequences of alkanes, alcohols, and aldehydes. For simplicity, the analysis focuses on abstractions from the fuel by H, CH 3 , and OH, and the decomposition of the resulting radicals. Altogether, there are a total of 166 input channels in these sets (more than 363 forward reactions when expanded to the full set of elementary reactions). The code successfully produces high quality rate estimates (with apparent uncertainties less than a factor of two in limited comparisons with experiment) for > 95% of these. For the radical decomposition reactions, the analysis includes predictions for the pressure dependence of the kinetics. This wide-ranging exploration illustrates (i) the effect of different levels of prediction on the expected accuracy, (ii) the branching between abstractions at different sites for different abstractors, (iii) the dependence of the rates on the chemical structure, and (iv) the variation in radical stabilities across chemical families. These results, as well as the demonstrated feasibility of the methodology, should find further utility in the development of accurate rate expressions for arbitrary fuels. (c) 2020 The Combustion Institute. Published by Elsevier Inc. All rights reserved.
Elliott, Sarah N.↗
Empirical Characterization and Modeling of Cohesive – to – Adhesive Shear Fracture Mode Transition due to Increased Adhesive Layer Thicknesses of Fiber Reinforced Composite Single – Lap Joints
Here, to ensure a strong adhesive bond, most standards and adhesive manufacturers specify a maximum adhesive gap of 1 mm when bonding fiber reinforced composite structures. In manufacturing large components, such as joining two halves of wind turbine blades, meeting this gap tolerance specification is impractical; gaps larger than 10 mm are common in large adhesively bonded composite structures using state-of-the-art manufacturing techniques. Currently, there is a lack of fundamental understanding of the failure mechanics of adhesive gaps larger than 3 mm. To create such understanding, glass fiber - acrylic thermoplastic composite panels bonded using different epoxy adhesives within single-lap joint samples with adhesive thicknesses of 0.1 mm, 0.3 mm, 1 mm, 3 mm, 5 mm, and 10 mm were sheared to failure. A transition from cohesive to adhesive failure was observed to occur about 1 mm to 3 mm joint thicknesses. Plotting the shear stress normalized by the ratio of the joint width to thickness as a function of the joint thickness normalized by the joint length is shown to result in the ability to fit simple empirically derived models of the cohesive-to-adhesive failure transition, regardless of the adhesive. Furthermore, using these normalized variables, all the observed cohesively failed specimens collapse to a single master curve, as do the adhesively failed specimens.
36 MATERIALS SCIENCE↗
A New Objective Technique for Verifying Mesoscale Numerical Weather Prediction Models
This report presents a new objective technique to verify predictions of the sea-breeze phenomenon over east-central Florida by the Regional Atmospheric Modeling System (RAMS) mesoscale numerical weather prediction (NWP) model. The Contour Error Map (CEM) technique identifies sea-breeze transition times in objectively-analyzed grids of observed and forecast wind, verifies the forecast sea-breeze transition times against the observed times, and computes the mean post-sea breeze wind direction and speed to compare the observed and forecast winds behind the sea-breeze front. The CEM technique is superior to traditional objective verification techniques and previously-used subjective verification methodologies because: It is automated, requiring little manual intervention, It accounts for both spatial and temporal scales and variations, It accurately identifies and verifies the sea-breeze transition times, and It provides verification contour maps and simple statistical parameters for easy interpretation. The CEM uses a parallel lowpass boxcar filter and a high-order bandpass filter to identify the sea-breeze transition times in the observed and model grid points. Once the transition times are identified, CEM fits a Gaussian histogram function to the actual histogram of transition time differences between the model and observations. The fitted parameters of the Gaussian function subsequently explain the timing bias and variance of the timing differences across the valid comparison domain. Once the transition times are all identified at each grid point, the CEM computes the mean wind direction and speed during the remainder of the day for all times and grid points after the sea-breeze transition time. The CEM technique performed quite well when compared to independent meteorological assessments of the sea-breeze transition times and results from a previously published subjective evaluation. The algorithm correctly identified a forecast or observed sea-breeze occurrence or absence 93% of the time during the two- month evaluation period from July and August 2000. Nearly all failures in CEM were the result of complex precipitation features (observed or forecast) that contaminated the wind field, resulting in a false identification of a sea-breeze transition. A qualitative comparison between the CEM timing errors and the subjectively determined observed and forecast transition times indicate that the algorithm performed very well overall. Most discrepancies between the CEM results and the subjective analysis were again caused by observed or forecast areas of precipitation that led to complex wind patterns. The CEM also failed on a day when the observed sea- breeze transition affected only a very small portion of the verification domain. Based on the results of CEM, the RAMS tended to predict the onset and movement of the sea-breeze transition too early and/or quickly. The domain-wide timing biases provided by CEM indicated an early bias on 30 out of 37 days when both an observed and forecast sea breeze occurred over the same portions of the analysis domain. These results are consistent with previous subjective verifications of the RAMS sea breeze predictions. A comparison of the mean post-sea breeze winds indicate that RAMS has a positive wind-speed bias for .all days, which is also consistent with the early bias in the sea-breeze transition time since the higher wind speeds resulted in a faster inland penetration of the sea breeze compared to reality.
Case, Jonathan L.↗
Shear and bulk viscosity for a pure glue theory using an effective matrix model
At nonzero temperatures, the deconfining phase transition can be analyzed using an effective matrix model to characterize the change in holonomy. The model includes gluons and two-dimensional ghost fields in the adjoint representation, or “teens.” As ghosts, the teen fields are responsible for the decrease of the pressure as 𝑇 →𝑇 𝑑 , with 𝑇 𝑑 the transition temperature for deconfinement. Using the solution of this matrix model for a large number of colors, the parameters of the teen fields are adjusted so that the expectation value of the Polyakov loop is close to the values from the lattice. The shear, 𝜂, and bulk, 𝜁, viscosities are computed at nonzero holonomy to leading logarithmic order in weak coupling. In the pure glue theory, the value of the Polyakov loop is relatively large in the deconfined phase, ≈1/2 at 𝑇 𝑑 . Consequently, if 𝑠 is the entropy density, while 𝜂/𝑠 decreases as 𝑇 →𝑇 𝑑 , it is still well above the conformal bound. In contrast, 𝜁/𝑠 is largest at 𝑇 𝑑 , comparable to 𝜂/𝑠, then falls off rapidly with increasing temperature and is negligible by ∼2𝑇 𝑑 .
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Lessons learned in transitioning to an open systems environment
Software development organizations, both commercial and governmental, are undergoing rapid change spurred by developments in the computing industry. To stay competitive, these organizations must adopt new technologies, skills, and practices quickly. Yet even for an organization with a well-developed set of software engineering models and processes, transitioning to a new technology can be expensive and risky. Current industry trends are leading away from traditional mainframe environments and toward the workstation-based, open systems world. This paper presents the experiences of software engineers on three recent projects that pioneered open systems development for NASA's Flight Dynamics Division of the Goddard Space Flight Center (GSFC).
Boland, Dillard E.↗
Sign Problem in Tensor-Network Contraction
We investigate how the computational difficulty of contracting tensor networks depends on the sign structure of the tensor entries. Using results from computational complexity, we observe that the approximate contraction of tensor networks with only positive entries has lower computational complexity as compared to tensor networks with general real or complex entries. This raises the question of how this transition in computational complexity manifests itself in the hardness of different tensor-network-contraction schemes. We pursue this question by studying random tensor networks with varying bias toward positive entries. First, we consider contraction via Monte Carlo sampling and find that the transition from hard to easy occurs when the tensor entries become predominantly positive; this can be understood as a tensor-network manifestation of the well-known negative-sign problem in quantum Monte Carlo. Second, we analyze the commonly used contraction based on boundary tensor networks. The performance of this scheme is governed by the number of correlations in contiguous parts of the tensor network (which by analogy can be thought of as entanglement). Remarkably, we find that the transition from hard to easy—i.e., from a volume-law to a boundary-law scaling of entanglement—already occurs for a slight bias of the tensor entries toward a positive mean, scaling inversely with the bond dimension D , and thus the problem becomes easy the earlier the larger D occurs. This is in contrast both to expectations and to the behavior found in Monte Carlo contraction, where the hardness at fixed bias increases with the bond dimension. To provide insight into this early breakdown of computational hardness and the accompanying entanglement transition, we construct an effective classical statistical-mechanical model that predicts a transition at a bias of the tensor entries of 1 / D , confirming our observations. We conclude by investigating the computational difficulty of computing expectation values of tensor-network wave functions (projected entangled-pair states, PEPSs) and find that in this setting, the complexity of entanglement-based contraction always remains low. We explain this by providing a local transformation that maps PEPS expectation values to a positive-valued tensor network. This not only provides insight into the origin of the observed boundary-law entanglement scaling but also suggests new approaches toward PEPS contraction based on positive decompositions. Published by the American Physical Society 2025
Chen, Jielun (ORCID:0000000178411545)↗
The structure and heating of the chromosphere-corona transition region
The structure and heating (or energy balance) of the transition region and the role of the transition region in the structure and heating of the solar atmosphere as a whole are investigated. The features of the structure of the atmosphere and radiative energy losses of the atmosphere are summarized. A static, planar model of the solar temperature which has a temperature profile representative of the actual solar atmosphere is considered. Then a static, planar model of the transition region which is heated by thermal conduction from the corona and cooled by radiative losses is developed. A general conclusion is that the temperature profile of the transition region and lower corona results primarily from the energy balance of the corona, while the number density is determined by the energy balance of the transition region.
Moore, R. L.↗
An Experimental Investigation of Transition as Applied to Low Pressure Turbine Suction Surface Flows
Results of an experimental study of flow separation and transition in either attached boundary layers or separated shear layers over the suction surface of a simulation of a low-pressure turbine airfoil flow are presented. Detailed velocity profiles were measured with the hot-wire technique. Static pressure distributions are also presented. Flow transition is documented using measured intermittency distributions in the boundary layer and the separated shear layer. Cases for Reynolds numbers of 50,000, 100,000, 200,000 and 300,000 are reported. These Reynolds numbers are based on suction surface length and exit velocity. Three Free Stream Turbulence Intensity values, 0.5%, 2.5% and 10%, are represented. Flow separation is observed for all the low-FSTI cases. Of these, the lowest Reynolds number case was not able to complete transition of the shear layer and the separation bubble persisted over the entire blade surface. For the other low-FSTI cases, transition is observed in the shear layer over the separation bubble. This transition proceeded quickly, spreading rapidly toward the wall. Elevated FSTI drives an earlier transition than in the low-FSTI cases and the separation bubbles are smaller. For the highest Reynolds number cases with 2.5% and 10% FSTI, transition is of the attached boundary layer and no separation exists. Flow separation with shear flow transition is observed for the lower-Re cases. Models for intermittency and transition length and location from the modern literature are assessed.
Qiu, Songgang↗
Coupling of the coronal helium abundance to the solar wind
Models of the transition region-corona-solar wind system are investigated in order to find the coronal helium abundance and to study the role played by coronal helium in controlling the solar wind proton flux. The thermal force on alpha-particles in the transition region sets the flow of helium into the corona. The frictional coupling between alpha-particles and protons and/or the electric polarization field determines the proton flux in the solar wind as well as the fate of the coronal helium content. The models are constructed by solving the time-dependent population and momentum equations for all species of hydrogen and helium in an atmosphere with a given temperature profile. Several temperature profiles are considered in order to very the roles of frictional coupling and electric polarization field in the solar wind, and the thermal force in the transition region. Steady-state solutions are found for coronae with a hydrogen flux at 1 AU of 1.0 x 10(exp 9)/cm(exp 2)/sec or larger. For coronae with lower hydrogen fluxes, the helium flux into the corona is larger than the flux 'pulled out' by the solar wind protons, and solutions with increasing coronal helium content are found. The timescale for forming a helium-filled corona, that may allow for a steady outflow, is long compared to the mixing time for the corona.
Hansteen, Viggo H.↗
Status of the TESS Science Processing Operations Center
The Transiting Exoplanet Survey Satellite (TESS) science pipeline is being developed by the Science Processing Operations Center (SPOC) at NASA Ames Research Center based on the highly successful Kepler Mission science pipeline. Like the Kepler pipeline, the TESS science pipeline will provide calibrated pixels, simple and systematic error-corrected aperture photometry, and centroid locations for all 200,000+ target stars, observed over the 2-year mission, along with associated uncertainties. The pixel and light curve products are modeled on the Kepler archive products and will be archived to the Mikulski Archive for Space Telescopes (MAST). In addition to the nominal science data, the 30-minute Full Frame Images (FFIs) simultaneously collected by TESS will also be calibrated by the SPOC and archived at MAST. The TESS pipeline will search through all light curves for evidence of transits that occur when a planet crosses the disk of its host star. The Data Validation pipeline will generate a suite of diagnostic metrics for each transit-like signature discovered, and extract planetary parameters by fitting a limb-darkened transit model to each potential planetary signature. The results of the transit search will be modeled on the Kepler transit search products (tabulated numerical results, time series products, and pdf reports) all of which will be archived to MAST.
high performance computing↗