Results for “linearized distribution” — Page 96

NASA NTRS1986

A three-dimensional dynamic analysis and libration study of a tethered satellites system

By a three-dimensional Lagrangian formulation, an exact set of ten coupled nonlinear second-order differential equations has been derived for a system with an extensible tether connecting two end satellites of distributed mass. The effects of tether mass, small orbital eccentricity, central body oblateness, aerodynamic drag force, and solar radiation pressure are also included in the formulation. By linearizing the exact differential equations, the in-plane (orbital plane) differential equations are found to be decoupled from the out-of-plane ones. The characteristic equation of the in-plane differential equations is derived and some associated stability constraints are shown.

Tang, C. C. H.↗

NASA NTRS1987

A Comparison of Experimental and Theoretical Results for Labyrinth Gas Seals

The basic equations are derived for a two control volume model for compressible flow in a labyrinth seal. The flow is assumed to be completely turbulent and isoenergetic. The wall friction factors are determined using the Blasius formula. Jet flow theory is used for the calculation of the recirculation velocity in the cavity. Linearized zeroth and first order perturbation equations are developed for small motion about a centered position by an expansion in the eccentricity ratio. The zeroth order pressure distribution is found by satisfying the leakage equation. The circumferential velocity distribution is determined by satisfying the momentum equations. The first order equations are solved by a separation of variable solution. Integration of the resultant pressure distribution along and around the seal defines the reaction force developed by the seal and the corresponding dynamic coefficients. The results of this analysis are compared to experimental test results.

Scharrer, Joseph Kirk↗

NASA NTRS1987

New methods and results for quantification of lightning-aircraft electrodynamics

The NASA F-106 collected data on the rates of change of electromagnetic parameters on the aircraft surface during over 700 direct lightning strikes while penetrating thunderstorms at altitudes from 15,000 t0 40,000 ft (4,570 to 12,190 m). These in situ measurements provided the basis for the first statistical quantification of the lightning electromagnetic threat to aircraft appropriate for determining indirect lightning effects on aircraft. These data are used to update previous lightning criteria and standards developed over the years from ground-based measurements. The proposed standards will be the first which reflect actual aircraft responses measured at flight altitudes. Nonparametric maximum likelihood estimates of the distribution of the peak electromagnetic rates of change for consideration in the new standards are obtained based on peak recorder data for multiple-strike flights. The linear and nonlinear modeling techniques developed provide means to interpret and understand the direct-strike electromagnetic data acquired on the F-106. The reasonable results obtained with the models, compared with measured responses, provide increased confidence that the models may be credibly applied to other aircraft.

Pitts, Felix L.↗

NASA NTRS1986

Numerical solution of random singular integral equation appearing in crack problems

The solution of several elasticity problems, and particularly crack problems, can be reduced to the solution of one-dimensional singular integral equations with a Cauchy-type kernel or to a system of uncoupled singular integral equations. Here a method for the numerical solution of random singular integral equations of Cauchy type is presented. The solution technique involves a Chebyshev series approximation, the coefficients of which are the solutions of a system of random linear equations. This method is applied to the problem of periodic array of straight cracks inside an infinite isotropic elastic medium and subjected to a nonuniform pressure distribution along the crack edges. The statistical properties of the random solution are evaluated numerically, and the random solution is used to determine the values of the stress-intensity factors at the crack tips. The error, expressed as the difference between the mean of the random solution and the deterministic solution, is established. Values of stress-intensity factors at the crack tip for different random input functions are presented.

Sambandham, M.↗

NASA NTRS1987

Stochastic estimation of organized structures in turbulent channel flow

Linear stochastic estimation is used to approximate the conditional vector fields associated with high Reynolds stress producing events in numerically simulated turbulent channel flow. Joint probability density distributions of u-v are presented for y(+) between 1.35 and 180, and at each position the u-v values that make the greatest contributions to the average uv3 are used to define conditional ejection- and sweep-type events. In the y-z plane the conditional fields appear to be pairs of counter-rotating vortices. When the conditional event defined by the weighted probability density function is specified at a point close to the wall, the spacing and size of these eddies are consistent with the accepted wall streak spacing. As this point moves outwards, the eddy size increases significantly. An abrupt change in the flow angle occurs in the buffer layer and may indicate transition from streamwise oriented wall layer structures to hairpin vortices characterizing the outer layer.

Moin, P.↗

NASA NTRS1988

Load balancing and closed chain multiple arm control

The authors give the general dynamical equations for several rigid link manipulators rigidly grasping a commonly held rigid object. It is shown that the number of arm-configuration degrees of freedom lost due to imposing the closed-loop kinematic constraints is the same as the number of degrees of freedom gained for controlling the internal forces of the closed-chain system. This number is equal to the dimension of the kernel of the Jacobian operator which transforms contact forces to the net forces acting on the held object, and it is shown that this kernel can be identified with the subspace of controllable internal forces of the closed-chain system. Control of these forces makes it possible to regulate the grasping forces imparted to the held object or to control the load taken by each arm. It is shown that the internal forces can be influenced without affecting the control of the configuration degrees of freedom. Control laws of the feedback linearization type are shown to be useful for controlling the location and attitude of a frame fixed with respect to the held object, while simultaneously controlling the internal forces of the closed-chain system. Force feedback can be used to linearize and control the system even when the held object has unknown mass properties. If saturation effects are ignored, an unconstrained quadratic optimization can be performed to distribute the load optimally among the joint actuators.

Kreutz, Kenneth↗

NASA NTRS1991

Photon interaction cross sections for aggregations of finite-dimensional leaves

The Ross (1981) plate turbid-medium theory, which abstracts the leaf canopy as a gaseous mixture of oriented planar nondimensional plates that are randomly distributed in the configuration space, does not account for the 'hot spot' effect typical of all layered media with finite-dimensional scattering centers. The concept of 'particle distribution functions' is presently used to derive the interaction cross sections. This formalism is statistical-mechanical in nature, since the radiative transfer equation is essentially a linearized form of the Boltzmann equation for N-body systems.

Myneni, Ranga B.↗

NASA NTRS1994

Planetary-scale circulations in the presence of climatological and wave-induced heating

Interaction between the large-scale circulation and the convective pattern is investigated in a coupled system governed by the linearized primitive equations. Convection is represented in terms of two components of heating: A 'climatological component' is prescribed stochastically to represent convection that is maintained by fixed distributions of land and sea and sea surface temperature (SST). An 'induced component' is defined in terms of the column-integrated moisture flux convergence to represent convection that is produced through feedback with the circulation. Each component describes the envelope organizing mesoscale convective activity. As SST on the equator is increased, induced heating amplifies in the gravest zonal wavenumbers at eastward frequencies, where positive feedback offsets dissipation. Under barotropic stratification, a critical SST of 29.5 C results in positive feedback exactly cancelling dissipation in wavenumber 1 for an eastward phase speed of 6 m/s. Sympathetic interaction between the circulation and the induced heating is the basis for 'frictional wave-Conditional Instability of the Second Kind (CISK)', which is distinguished from classical wave-CISK by rendering the gravest zonal dimensions most unstable. Under baroclinic stratification, the coupled system exhibits similar behavior. The critical SST is only 26.5 C for conditions representative of equinox, but in excess of 30 C for conditions representative of solstice. Having the form of an unsteady Walker circulation, the disturbance produced by frictional wave-CISK compares favorably with the observed life cycle of the Madden-Julian oscillation (MJO). SST above the critical value produces an amplifying disturbance in which enhanced convection coincides with upper-tropospheric westerlies and is positively correlated with temperature and surface convergence. Conversely, SST below the critical value produces a decaying disturbance in which enhanced convection coincides with upper-tropospheric easterlies and is nearly in quadrature with temperature and surface convergence. While sharing essential features with the MJO in the Eastern Hemisphere, frictional wave-CISK does not explain observed behavior in the Western Hemisphere, where the convective signal is largely absent. Comprised of Kelvin structure with the same frequency, observed behavior in the Western Hemisphere can be understood as a propagating response that is excited in and radiates away from the fluctuation of convection in the Eastern Hemisphere.

Salby, Murry L↗

NASA NTRS2005

Using Remote Sensing Data to Evaluate Surface Soil Properties in Alabama Ultisols

Evaluation of surface soil properties via remote sensing could facilitate soil survey mapping, erosion prediction and allocation of agrochemicals for precision management. The objective of this study was to evaluate the relationship between soil spectral signature and surface soil properties in conventionally managed row crop systems. High-resolution RS data were acquired over bare fields in the Coastal Plain, Appalachian Plateau, and Ridge and Valley provinces of Alabama using the Airborne Terrestrial Applications Sensor multispectral scanner. Soils ranged from sandy Kandiudults to fine textured Rhodudults. Surface soil samples (0-1 cm) were collected from 163 sampling points for soil organic carbon, particle size distribution, and citrate dithionite extractable iron content. Surface roughness, soil water content, and crusting were also measured during sampling. Two methods of analysis were evaluated: 1) multiple linear regression using common spectral band ratios, and 2) partial least squares regression. Our data show that thermal infrared spectra are highly, linearly related to soil organic carbon, sand and clay content. Soil organic carbon content was the most difficult to quantify in these highly weathered systems, where soil organic carbon was generally less than 1.2%. Estimates of sand and clay content were best using partial least squares regression at the Valley site, explaining 42-59% of the variability. In the Coastal Plain, sandy surfaces prone to crusting limited estimates of sand and clay content via partial least squares and regression with common band ratios. Estimates of iron oxide content were a function of mineralogy and best accomplished using specific band ratios, with regression explaining 36-65% of the variability at the Valley and Coastal Plain sites, respectively.

Sullivan, Dana G.↗

NASA NTRS2010

Time-Dependent Selection of an Optimal Set of Sources to Define a Stable Celestial Reference Frame

Temporal statistical position stability is required for VLBI sources to define a stable Celestial Reference Frame (CRF) and has been studied in many recent papers. This study analyzes the sources from the latest realization of the International Celestial Reference Frame (ICRF2) with the Allan variance, in addition to taking into account the apparent linear motions of the sources. Focusing on the 295 defining sources shows how they are a good compromise of different criteria, such as statistical stability and sky distribution, as well as having a sufficient number of sources, despite the fact that the most stable sources of the entire ICRF2 are mostly in the Northern Hemisphere. Nevertheless, the selection of a stable set is not unique: studying different solutions (GSF005a and AUG24 from GSFC and OPA from the Paris Observatory) over different time periods (1989.5 to 2009.5 and 1999.5 to 2009.5) leads to selections that can differ in up to 20% of the sources. Observing, recording, and network improvement are some of the causes, showing better stability for the CRF over the last decade than the last twenty years. But this may also be explained by the assumption of stationarity that is not necessarily right for some sources.

Le Bail, Karine↗

NASA NTRS2011

A Pure Marine Aerosol Model, for Use in Remote Sensing Applications

Retrievals of aerosol optical depth (AOD) and related parameters from satellite measurements typically involve prescribed models of aerosol size and composition, and are therefore dependent on how well these models are able to represent the radiative behaviour of real aerosols, This study uses aerosol volume size distributions retrieved from Sun-photometer measurements at 11 Aerosol Robotic Network (AERONET) island sites, spread throughout the world's oceans, as a basis to define such a model for unpolluted maritime aerosols. Size distributions are observed to be bimodal and approximately lognormal, although the coarse mode is skewed with a long tail on the low-radius end, The relationship of AOD and size distribution parameters to meteorological conditions is also examined, As wind speed increases, so do coarse-mode volume and radius, The AOD and Angstrom exponent (alpha) show linear relationships with wind speed, although there is considerable scatter in all these relationships, limiting their predictive power. Links between aerosol properties and near-surface relative humidity, columnar water vapor, and sea surface temperature are also explored. A recommended bimodal maritime model, which is able to reconstruct the AERONET AOD with accuracy of order 0.01-0.02, is presented for use in aerosol remote sensing applications. This accuracy holds at most sites and for wavelengths between 340 nm and 1020 nm. Calculated lidar ratios are also provided, and differ significantly from those currently used in Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) processing.

Sayer, A. M.↗

NASA NTRS2015

An Experimental Investigation of the Flow Physics Associated With End Wall Losses and Large Rotor Tip Clearances as Found in the Rear Stages of a High Pressure Compressor

The focus of this work was to characterize the fundamental flow physics and the overall performance effects due to increased rotor tip clearance heights in axial compressors. Data have been collected in the three-stage axial research compressor at Purdue University with a specific focus on analyzing the multistage effects resulting from the tip leakage flow. Three separate rotor tip clearance heights were studied with nominal tip clearance heights of 1.5%, 3.0%, and 4.0% based on a constant annulus height. Overall compressor performance was investigated at four corrected speedlines (100%, 90%, 80%, and 68%) for each of the three tip clearance configurations using total pressure and total temperature rakes distributed throughout the compressor. The results have confirmed results from previous authors showing a decrease of total pressure rise, isentropic efficiency, and stall margin which is approximately linear with increasing tip clearance height. The stall inception mechanisms have also been evaluated at the same corrected speeds for each of the tip clearance configurations. Detailed flow field measurements have been collected at two loading conditions, nominal loading (NL) and high loading (HL), on the 100% corrected speedline for the smallest and largest tip clearance heights (1.5% and 4.0%). Steady detailed radial traverses of total pressure at the exit of each stator row have been supported by flow visualization techniques to identify regions of flow recirculation and separation. Furthermore, detailed radial traverses of time-resolved total pressures at the exit of each rotor row have been measured with a fast-response pressure probe. These data have helped to quantify the size of the leakage flow at the exit of each rotor. Thermal anemometry has also been implemented to evaluate the time-resolved three-dimensional components of velocity throughout the compressor and calculate blockage due to the rotor tip leakage flow throughout the compressor. These measurements have also been used to calculate streamwise vorticity. Time-resolved static pressure measurements have been collected over the rotor tips for all rotors with each of the three tip clearance configurations for up to five loading conditions along the 100% corrected speedline using fast-response piezoresistive pressure sensors. These time-resolved static pressure measurements, as well as the time-resolved total pressures and velocities have helped to reveal a profound influence of the upstream stator vane on the size and shape of the rotor tip leakage flow. Finally, a novel particle image velocimetry (PIV) technique has been developed as a proof-of- concept. In contrast to PIV methods that have been typically been utilized for turbomachinery applications in the past, the method used for this study introduced the laser light through the same access window that was also used to image the flow. This new method addresses potential concerns related to the intrusive laser-introducing techniques that have typically been utilized by other authors in the past. Ultimately, the data collected for this project represent a unique data set which contributes to build a better understanding of the tip leakage flow field and its associated loss mechanisms. These data will facilitate future engine design goals leading to small blade heights in the rear stages of high pressure compressors and aid in the development of new blade designs which are desensitized to the performance penalties attributed to rotor tip leakage flows.

turbomachinery↗

NASA NTRS2015

An Experimental Investigation of the Flow Physics Associated With End Wall Losses and Large Rotor Tip Clearances as Found in the Rear Stages of a High Pressure Compressor

The focus of this work was to characterize the fundamental flow physics and the overall performance effects due to increased rotor tip clearance heights in axial compressors. Data have been collected in the three-stage axial research compressor at Purdue University with a specific focus on analyzing the multistage effects resulting from the tip leakage flow. Three separate rotor tip clearance heights were studied with nominal tip clearance heights of 1.5%, 3.0%, and 4.0% based on a constant annulus height. Overall compressor performance was investigated at four corrected speedlines (100%, 90%, 80%, and 68%) for each of the three tip clearance configurations using total pressure and total temperature rakes distributed throughout the compressor. The results have confirmed results from previous authors showing a decrease of total pressure rise, isentropic efficiency, and stall margin which is approximately linear with increasing tip clearance height. The stall inception mechanisms have also been evaluated at the same corrected speeds for each of the tip clearance configurations. Detailed flow field measurements have been collected at two loading conditions, nominal loading (NL) and high loading (HL), on the 100% corrected speedline for the smallest and largest tip clearance heights (1.5% and 4.0%). Steady detailed radial traverses of total pressure at the exit of each stator row have been supported by flow visualization techniques to identify regions of flow recirculation and separation. Furthermore, detailed radial traverses of time-resolved total pressures at the exit of each rotor row have been measured with a fast-response pressure probe. These data have helped to quantify the size of the leakage flow at the exit of each rotor. Thermal anemometry has also been implemented to evaluate the time-resolved three-dimensional components of velocity throughout the compressor and calculate blockage due to the rotor tip leakage flow throughout the compressor. These measurements have also been used to calculate streamwise vorticity. Time-resolved static pressure measurements have been collected over the rotor tips for all rotors with each of the three tip clearance configurations for up to five loading conditions along the 100% corrected speedline using fast-response piezoresistive pressure sensors. These time-resolved static pressure measurements, as well as the time-resolved total pressures and velocities have helped to reveal a profound influence of the upstream stator vane on the size and shape of the rotor tip leakage flow. Finally, a novel particle image velocimetry (PIV) technique has been developed as a proof-of- concept. In contrast to PIV methods that have been typically been utilized for turbomachinery applications in the past, the method used for this study introduced the laser light through the same access window that was also used to image the flow. This new method addresses potential concerns related to the intrusive laser-introducing techniques that have typically been utilized by other authors in the past. Ultimately, the data collected for this project represent a unique data set which contributes to build a better understanding of the tip leakage flow field and its associated loss mechanisms. These data will facilitate future engine design goals leading to small blade heights in the rear stages of high pressure compressors and aid in the development of new blade designs which are desensitized to the performance penalties attributed to rotor tip leakage flows.

Berdanier, Reid A.↗

DOE OSTI2024

Measurement of the t-dependence of the Beam Asymmetry of Photoproduced eta at the GlueX Experiment at Jefferson Lab

We report on the photoproduction of eta mesons studied at the GlueX experiment at Thomas Jefferson National Laboratory in Newport News, VA. A linearly polarized photon beam at energies between 8.2 and 8.8 GeV, incident on a liquid hydrogen target, produces eta dominantly through t-channel exchanges. Azimuthal angular distributions of the eta, for the direction of the polarized photon, facilitate the extraction of the beam asymmetry Sigma for the reaction gamma p -> eta p. We can understand the production mechanisms of the eta by analyzing Sigma_eta, which is sensitive to the particles exchanged. Compared to previous GlueX results [[4],[17]], we have access to five times the statistics, thereby allowing us to extend these measurements to values beyond the previous limitation of -t lte 1.1 (GeV)2. Results will be shown for events reconstructed from the decay of eta -> gamma gamma and eta -> pi^+ pi^- pi^0 Systematic checks are also made for the neutral mode decay to confirm the stability of the extracted results. We find through data binned up to a momentum transfer of 2.0 (GeV)2 that the photoproduction of eta at these beam energies is largely carried out by natural processes, with a preliminary comparison with the photoproduction of eta' suggesting possible strangeness exchanges at wider kinematic regions.

Erbora, Tolga↗

DOE OSTI2025

Sub-microsecond Transformers for Jet Tagging on FPGAs

We present the first sub-microsecond transformer implementation on an FPGA achieving competitive performance for state-of-the-art high-energy physics benchmarks. Transformers have shown exceptional performance on multiple tasks in modern machine learning applications, including jet tagging at the CERN Large Hadron Collider (LHC). However, their computational complexity prohibits use in real-time applications, such as the hardware trigger system of the collider experiments up until now. In this work, we demonstrate the first application of transformers for jet tagging on FPGAs, achieving $\mathcal{O}(100)$ nanosecond latency with superior performance compared to alternative baseline models. We leverage high-granularity quantization and distributed arithmetic optimization to fit the entire transformer model on a single FPGA, achieving the required throughput and latency. Furthermore, we add multi-head attention and linear attention support to hls4ml, making our work accessible to the broader fast machine learning community. This work advances the next-generation trigger systems for the High Luminosity LHC, enabling the use of transformers for real-time applications in high-energy physics and beyond.

Laatu, Lauri [Imperial Coll., London]↗

DOE OSTI2021

Framework for optimization of long-term, multi-period investment planning of integrated urban energy systems

In order to achieve stringent greenhouse gas emission reductions, a transition of our entire energy system from fossil to renewable resources needs to be designed. Such an energy transition brings two main challenges: most renewables generate variable electric energy, yet most demand is currently not electric (carrier mismatch) and does not always manifest at the same time as supply (temporal mismatch). Integrating multiple energy infrastructures can address both challenges by using the synergy between different energy carriers; building on existing infrastructure, while allowing a robust and flexible integration of the new. This paper proposes an optimization framework for long-term, multi-period investment planning of urban energy systems in an integrated manner. We formulate it as a mixed-integer linear program, combining a capacitated facility location with a multi-dimensional, capacitated network design problem. It includes generation and network expansion planning as well as interconnections between networks and storage infrastructure for each energy system. It can incorporate pathway effects like techno-economic developments, policy measures, and weather variations. The intended use is to support urban decision makers with long-term investment planning, though it can be tailored to fit other geographical or temporal scales. We demonstrate the model using two cases based on an average city in The Netherlands, which wants to reduce its CO 2 -emissions with 95% by 2050. In the first case, we include explicit carbon-emission constraints to study the effects of the carrier mismatch. In the second case, we implement interannual weather variations to analyze the temporal mismatch. The results give valuable insights into the energy transition design strategy for urban decision makers. They also show the future potential, as well as the computational challenges of the optimization framework.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2024

Revealing Decision Conservativeness Through Inverse Distributionally Robust Optimization

This paper introduces Inverse Distributionally Robust Optimization (I-DRO) as a method to infer the conservativeness level of a decision-maker, represented by the size of a Wasserstein metric-based ambiguity set, from the optimal decisions made using Forward Distributionally Robust Optimization (F-DRO). By leveraging the Karush-Kuhn-Tucker (KKT) conditions of the convex F-DRO model, we formulate I-DRO as a bi-linear program, which can be solved using off-the-shelf optimization solvers. Additionally, this formulation exhibits several advantageous properties. We demonstrate that I-DRO not only guarantees the existence and uniqueness of an optimal solution but also establishes the necessary and sufficient conditions for this optimal solution to accurately match the actual conservativeness level in F-DRO. Furthermore, we identify three extreme scenarios that may impact I-DRO effectiveness. Our case study applies F-DRO for power system scheduling under uncertainty and employs I-DRO to recover the conservativeness level of system operators. Numerical experiments based on an IEEE 5-bus system and a realistic NYISO 11-zone system demonstrate I-DRO performance in both normal and extreme scenarios. An extended version of this paper with additional analyses is available at li2024revealing.

distributionally robust optimization↗

NASA NTRS1943

Heat Transfer Through Turbulent Friction Layers

The "general Prandtl number" Pr(exp 1) - A(sub q)/A Pr, aside from the Reynolds number determines the ratio of turbulent to molecular heat transfer, and the temperature distribution in turbulent friction layers. A(sub q) = exchange coefficient for heat; A = exchange coefficient for momentum transfer. A formula is derived from the equation defining the general Prandtl number which describes the temperature as a function of the velocity. For fully developed thermal boundary layers all questions relating to heat transfer to and from incompressible fluids can be treated in a simple manner if the ratio of the turbulent shear stress to the total stress T(sub t)/T in the layers near the wall is known, and if the A(sub q)/A can be regarded as independent of the distance from the wall. The velocity distribution across a flat smooth channel and deep into the laminar sublayer was measured for isothermal flow to establish the shear stress ratio T(sub t)/T and to extend the universal wall friction law. The values of T(sub t)/T which resulted from these measurements can be approximately represented by a linear function of the velocity in the laminar-turbulent transition zone. The effect of the temperature relationship of the material values on the flow near the wall is briefly analyzed. It was found that the velocity at the laminar boundary (in contrast to the thickness of the laminar layer) is approximately independent of the temperature distribution. The temperature gradient at the wall and the distribution of temperature and heat flow in the turbulent friction layers were calculated on the basis of the data under two equations. The derived formulas and the figures reveal the effects of the Prandtl number, the Reynolds number, the exchange quantities and the temperature relationship of the material values.

Reichardt, H.↗