SEARCH · Engineering Papers

Results for “linearized distribution”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 595 records · Page 33

All sources NASA NTRS DOE OSTI

NASA NTRS1988

Approximation theory for LQG (Linear-Quadratic-Gaussian) optimal control of flexible structures

An approximation theory is presented for the LQG (Linear-Quadratic-Gaussian) optimal control problem for flexible structures whose distributed models have bounded input and output operators. The main purpose of the theory is to guide the design of finite dimensional compensators that approximate closely the optimal compensator. The optimal LQG problem separates into an optimal linear-quadratic regulator problem and an optimal state estimation problem. The solution of the former problem lies in the solution to an infinite dimensional Riccati operator equation. The approximation scheme approximates the infinite dimensional LQG problem with a sequence of finite dimensional LQG problems defined for a sequence of finite dimensional, usually finite element or modal, approximations of the distributed model of the structure. Two Riccati matrix equations determine the solution to each approximating problem. The finite dimensional equations for numerical approximation are developed, including formulas for converting matrix control and estimator gains to their functional representation to allow comparison of gains based on different orders of approximation. Convergence of the approximating control and estimator gains and of the corresponding finite dimensional compensators is studied. Also, convergence and stability of the closed-loop systems produced with the finite dimensional compensators are discussed. The convergence theory is based on the convergence of the solutions of the finite dimensional Riccati equations to the solutions of the infinite dimensional Riccati equations. A numerical example with a flexible beam, a rotating rigid body, and a lumped mass is given.

Gibson, J. S.↗

DOE OSTI2020

Machine learning approaches for structural and thermodynamic properties of a Lennard-Jones fluid

Predicting the functional properties of many molecular systems relies on understanding how atomistic interactions give rise to macroscale observables. However, current attempts to develop predictive models for the structural and thermodynamic properties of condensed-phase systems often rely on extensive parameter fitting to empirically selected functional forms whose effectiveness is limited to a narrow range of physical conditions. Here, we illustrate how these traditional fitting paradigms can be superseded using machine learning. Specifically, we use the results of molecular dynamics simulations to train machine learning protocols that are able to produce the radial distribution function, pressure, and internal energy of a Lennard-Jones fluid with increased accuracy in comparison to previous theoretical methods. The radial distribution function is determined using a variant of the segmented linear regression with the multivariate function decomposition approach developed by Craven et al. [J. Phys. Chem. Lett. 11, 4372 (2020)]. The pressure and internal energy are determined using expressions containing the learned radial distribution function and also a kernel ridge regression process that is trained directly on thermodynamic properties measured in simulation. The presented results suggest that the structural and thermodynamic properties of fluids may be determined more accurately through machine learning than through human-guided functional forms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1972

Kullback-Leibler information function and the sequential selection of experiments to discriminate among several linear models

The error variance of the process prior multivariate normal distributions of the parameters of the models are assumed to be specified, prior probabilities of the models being correct. A rule for termination of sampling is proposed. Upon termination, the model with the largest posterior probability is chosen as correct. If sampling is not terminated, posterior probabilities of the models and posterior distributions of the parameters are computed. An experiment was chosen to maximize the expected Kullback-Leibler information function. Monte Carlo simulation experiments were performed to investigate large and small sample behavior of the sequential adaptive procedure.

Sidik, S. M.↗

DOE OSTI2024

Minimum entropy filtering for a single output non-Gaussian stochastic system using state transformation

This paper presents a novel filter design for the single-output stochastic non-linear systems subjected to non-Gaussian noises and the proposed assumptions. Based on a state transformation, the unmeasurable states of the systems can be estimated where non-linear terms in the systems have been eliminated. It has been shown that the estimation error is linearly dynamical regarding to the presented vector-valued filter gain which can be optimised by minimising the entropy-based performance criterion. In addition, the convergence of the presented algorithm is analysed in mean-square sense and a numerical example is given to verify the effectiveness of the presented filtering algorithm. Meanwhile, the extended Kalman filter, unscented particle filter and minimum entropy filter are given for the comparisons of the filtering performance. Following the presented framework, some extensions of the presented filtering algorithm are discussed to indicate the flexibility of the filter design. The contribution of this paper can be summarised as establishing a novel minimum entropy filtering framework which consists of model transformation, entropy optimisation and convergence analysis.

42 ENGINEERING↗

NASA NTRS1995

Linear operating region in the ozone dial photon counting system

Ozone is a relatively unstable molecule found in Earth's atmosphere. An ozone molecule is made up of three atoms of oxygen. Depending on where ozone resides, it can protect or harm life on Earth. High in the atmosphere, about 15 miles up, ozone acts as a shield to protect Earth's surface from the sun's harmful ultraviolet radiation. Without this shield, we would be more susceptible to skin cancer, cataracts, and impaired immune systems. Closer to Earth, in the air we breathe, ozone is a harmful pollutant that causes damage to lung tissue and plants. Since the early 1980's, airborne lidar systems have been used for making measurements of ozone. The differential absorption lidar (DIAL) technique is used in the remote measurement of O3. This system allows the O3 to be measured as function of the range in the atmosphere. Two frequency-doubled Nd:YAG lasers are used to pump tunable dye lasers. The lasers are operating at 289 nm for the DIAL on-line wavelength of O3, and the other one is operated at 300 nm for the off-line wavelength. The DIAL wavelengths are produced in sequential laser pulses with a time separation of 300 micro s. The backscattered laser energy is collected by telescopes and measured using photon counting systems. The photon counting system measures the light signal by making use of the photon nature of light. The output pulse from the Photo-Multiplier Tube (PE), caused by a photon striking the PMT photo-cathode, is amplified and passed to a pulse height discriminator. The peak value of the pulse is compared to a reference voltage (discrimination level). If the pulse amplitude exceeds the discrimination level, the discriminator generates a standard pulse which is counted by the digital counter. Non-linearity in the system is caused by the overlapping of pulses and the finite response time of the electronics. At low count rates one expects the system to register one event for each output pulse from the PMT corresponding to a photon incident upon the photocathode, however, at higher rates the limitations of the discrimination/counting system will cause the observed count rate to be non-linear with respect to the true count rate. Depending on the pulse height distribution and the discriminator level, the overlapping of pulses (pulse pile-up) can cause count loss or even an additional apparent count gain as the signal levels increase. Characterization of the system, including the pulse height distribution, the signal to noise ratio, and the effect of the discriminator threshold level, is critical in maximizing the linear operating region of the system, thus greatly increasing the useful dynamic range of the system.

Andrawis, Madeleine↗

NASA NTRS2019

Evaluation of Wing Load Calibration and Sensing Methods Using Conventional Strain Gages and a Fiber Optic Sensing System Installed on a Straight Tapered Wing

State-of-the-art instrumentation techniques have provided an opportunity to obtain greater insight into structural wing loads during flight. An investigation was undertaken to research new wing instrumentation and load sensing techniques for measuring accurate in-flight spanwise load distributions on wing structures. A straight tapered wing was instrumented with both conventional foil strain gages at five spanwise wing stations and fiber optic strain sensors at every half inch along the entire wing span. Thirty-nine unique load cases were applied to the wing lower surface using hydraulic actuators to obtain various shear, bending moment, and torque load distributions on the wing. This paper will highlight three load calibration approaches. Conventional linear regression calibration methods were applied to foil strain gages providing a single wing station vertical shear, bending moment, and torque load. Linear regression methods were applied to a fiber optic sensing system to provide bending moment and torque spanwise load distributions. A load sensing scheme using strain derived wing shape information derived from a single load case provided vertical shear and bending moment spanwise load distribution information. Aspects of the three different approaches will be compared and contrasted to inform the reader of the benefits or disadvantages of each. Instrument installation, sensor characteristics, test execution aspects, and recommended calibration techniques will be discussed.

structural mechanics↗

DOE OSTI2020

Rupture Dynamics Simulations of Shallow Crustal Earthquakes on Reverse Slip Faults

The main goal of this study was to perform dynamic rupture simulations in order constrain the characterization of kinematic rupture models implemented in techniques for simulating strong ground motion for crustal earthquakes. First, we investigated the rupture process of the 2008 Iwate-Miyagi Nairiku earthquake by developing a dynamic rupture model on a reverse fault using a trial-and-error technique that produced a slip distribution and near-fault ground motion that matched the recorded ones. The simulations were performed in the frequency range 0-2 Hz, using a 3D staggered grid finite-difference method and a linear slip weakening friction law. Constrained by the observed slip distribution and consistent with dynamic rupture models, the derived kinematic rupture model of the Iwate-Miyagi Nairiku earthquake contains areas with relatively large slip rate, representing strong motion generation areas, set against lower amplitude heterogeneous background slip, and a relatively low slip rate in the weak zone of the top-most crust (upper 3km). Second, we performed rupture dynamics modeling to constrain shallow slip characterization in rupture models for strike-slip crustal earthquakes. The objective was to establish general rules about the characterization of slip rate function and slip at shallow depths and in the strong motion generation areas (SMGAs). The simulations of spontaneous rupture were performed in the frequency range 0-2.0Hz, using a 3D staggered-grid finite-difference method and a layer over half-space 1D crustal velocity model with a minimum shear-wave velocity of 2.8 km/s. In order to account for changes in material ductility and reduction of stress drop, observed in the shallow crust (upper 3-5 km), and the transition from ductile state to brittle state in the upper seismogenic zone, in our stress models we included a shallow weak zone (<4km). In this zone the stress drop was set to zero at the free surface and gradually increased with depth, while the slip weakening distance was set to 75 cm at the free surface and decreased to 50 cm at the base of the weak zone. As in the case of reverse faulting during the Iwate-Miyagi Nairiku earthquake, from these computations we found a systematic change in the shape of the slip-rate function, from Kostrov-type in the deeper part of the fault to a more symmetric cosine-type in the upper few kilometers, near the free-surface. Moreover, the average slip duration in the weak zone, with respect to slip duration in the deeper parts of the fault, increases by at most a factor of 1.5. We found a systematic and gradual change in the shape of the slip-rate function from Kostrov-type in the asperity areas (SMGAs), to more symmetric cosine-type in the upper few km near the free-surface, on the long period motion generation area (LMGA). The effective rise time in the LMGA, with respect to that in the SMGAs, increases by at most a factor of 2. Effective rise time is the time difference between the time at which the slip rate drops to a level that is equal to 25% of its peak and the onset time of the signal. In addition, the slip in the LMGA located above the SMGAs is almost the same as the one in the SMGAs, and about 1.5 times larger than the average slip.

58 GEOSCIENCES↗

NASA NTRS1992

Design, test, and evaluation of three active flutter suppression controllers

Three control law design techniques for flutter suppression are presented. Each technique uses multiple control surfaces and/or sensors. The first method uses traditional tools (such as pole/zero loci and Nyquist diagrams) for producing a controller that has minimal complexity and which is sufficiently robust to handle plant uncertainty. The second procedure uses linear combinations of several accelerometer signals and dynamic compensation to synthesize the model rate of the critical mode for feedback to the distributed control surfaces. The third technique starts with a minimum-energy linear quadratic Gaussian controller, iteratively modifies intensity matrices corresponding to input and output noise, and applies controller order reduction to achieve a low-order, robust controller. The resulting designs were implemented digitally and tested subsonically on the active flexible wing wind-tunnel model in the Langley Transonic Dynamics Tunnel. Only the traditional pole/zero loci design was sufficiently robust to errors in the nominal plant to successfully suppress flutter during the test. The traditional pole/zero loci design provided simultaneous suppression of symmetric and antisymmetric flutter with a 24-percent increase in attainable dynamic pressure. Posttest analyses are shown which illustrate the problems encountered with the other laws.

Adams, William M., Jr.↗

NASA NTRS1976

Determination of statistics for any rotation of axes of a bivariate normal elliptical distribution

Transformation of statistics from a dimensional set to another dimensional set involves linear functions of the original set of statistics. Similarly, linear functions will transform statistics within a dimensional set such that the new statistics are relevant to a new set of coordinate axes. A restricted case of the latter is the rotation of axes in a coordinate system involving any two correlated random variables. A special case is the transformation for horizontal wind distributions. Wind statistics are usually provided in terms of wind speed and direction (measured clockwise from north) or in east-west and north-south components. A direct application of this technique allows the determination of appropriate wind statistics parallel and normal to any preselected flight path of a space vehicle. Among the constraints for launching space vehicles are critical values selected from the distribution of the expected winds parallel to and normal to the flight path. These procedures are applied to space vehicle launches at Cape Kennedy, Florida.

Falls, L. W.↗

DOE OSTI2026

Probing gluon saturation with forward di-hadron correlations in proton-nucleus collisions

We present a detailed numerical investigation of semi-inclusive forward di-hadron production in proton–nucleus collisions employing the Color Glass Condensate effective theory. We focus on the regime where di-hadrons are produced nearly back-to-back in the transverse plane, thereby justifying a transverse-momentum-dependent factorization approach in terms of small-x gluon distributions. Our computation integrates several key elements: i) non-linear rapidity evolution via the Balitsky–Kovchegov equation with running coupling, ii) both perturbative and non-perturbative Sudakov resummation, and iii) a phenomenologically constrained model for the initial conditions for small-x gluon distributions. We compare this phenomenological framework to experimental data from the STAR Collaboration on azimuthal correlations in forward di-pion production in both proton–proton and proton–gold collisions. We analyze the systematic theoretical uncertainties associated with the saturation scales of nuclei at the initial scale for rapidity evolution and with those associated with the hadronization process. Finally, we make predictions for the kinematics anticipated to be covered by the ALICE Forward Calorimeter (FoCal) upgrade at the Large Hadron Collider.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1976

Dispersion of ensembles of non-interacting particles

The dynamics of an ensemble of noninteracting particles dispersing from a common origin and moving in a common force field with an initial distribution of momenta is analyzed using an approach where the particles are considered as a continuum described by a phase-space distribution function. General solutions are obtained for both the distribution function and the associated spatial density function. The linear case of small departures from circular orbits in an axisymmetric gravitational field is treated along with the specific case of particle dispersion from an object in a circular orbit in the same type of field. Numerical results are presented for the latter case, and consideration is given to the inverse problem of determining the initial time and velocity distribution from knowledge of the ensemble structure at a later time. Explicit results are provided for the case of an ellipsoidal distribution of initial momenta, and a numerical procedure is indicated for treating more general cases.

Heard, W. B.↗

NASA NTRS1984

Cosmic Ray Exposure & Linear-Energy-Transfer Evaluations for the COBE Mission

Magnetospherically attenuated, orbit integrated, surface incident cosmic ray fluxes of galactic origin were determined for the COBE mission. All heavy ions up to Nickel (z=28) were considered in this, evaluation. In order to provide worst case approximations, estimates were based on solar minimum conditions. Transport and material shielding calculations were then performed for these vehicle encountered particles, for a simple spherical 3-D aluminum geometry, and the materially attenuated cosmic ray distributions emerging behind selected shield thicknesses were obtained. Finally, Linear-Energy-Transfer (LET) spectra were evaluated for the energy spectra of each ion specie and these were, in turn, integrated into the sum-total LET distribution contained in this report. The results are presented in tabular and graphical form for eight (8) different shield thicknesses from.01 to 10.0 gm/cm2 of aluminum. They show most conclusively, and not unexpectedly, that material shielding has virtually no 'effect on the final LET spectrum, at least not for aluminum shields. This was predictable. Please note that in the LET graphs, two curves were plotted. Of these, please ignore the one labeled "grazing" and use only the one labeled "normal". Again, it should be remembered that the end results contain contributions from all elements up to Nickel (z=28).

Stassinopoulos, E. G.↗

NASA NTRS1986

Sensitivity of the Lidar ratio to changes in size distribution and index of refraction

In order to invert lidar signals to obtain reliable extinction coefficients, sigma, a relationship between sigma and the backscatter coefficient, beta, must be given. These two coefficients are linearly related if the complex index of refraction, m, particle shape size distribution, N, does not change along the path illuminated by the laser beam. This, however, is generally not the case. An extensive Mie computation of the lidar ratio R = beta/sigma and the sensitivity of R to the changes in a parametric space defined by N and m were examined.

Evans, B. T. N.↗

NASA NTRS1987

Surface flaw reliability analysis of ceramic components with the SCARE finite element postprocessor program

The SCARE (Structural Ceramics Analysis and Reliability Evaluation) computer program on statistical fast fracture reliability analysis with quadratic elements for volume distributed imperfections is enhanced to include the use of linear finite elements and the capability of designing against concurrent surface flaw induced ceramic component failure. The SCARE code is presently coupled as a postprocessor to the MSC/NASTRAN general purpose, finite element analysis program. The improved version now includes the Weibull and Batdorf statistical failure theories for both surface and volume flaw based reliability analysis. The program uses the two-parameter Weibull fracture strength cumulative failure probability distribution model with the principle of independent action for poly-axial stress states, and Batdorf's shear-sensitive as well as shear-insensitive statistical theories. The shear-sensitive surface crack configurations include the Griffith crack and Griffith notch geometries, using the total critical coplanar strain energy release rate criterion to predict mixed-mode fracture. Weibull material parameters based on both surface and volume flaw induced fracture can also be calculated from modulus of rupture bar tests, using the least squares method with known specimen geometry and grouped fracture data. The statistical fast fracture theories for surface flaw induced failure, along with selected input and output formats and options, are summarized. An example problem to demonstrate various features of the program is included.

Gyekenyesi, John P.↗

NASA NTRS1987

Surface flaw reliability analysis of ceramic components with the SCARE finite element postprocessor program

Gyekenyesi, J. P.↗

NASA NTRS1934

On the theory of laminar boundary layers involving separation

This paper presents a mathematical discussion of the laminar boundary layer, which was developed with a view of facilitating the investigation of those boundary layers in particular for which the phenomenon of separation occurs. The treatment starts with a slight modification of the form of the boundary layer equation first published by Von Mises. Two approximate solutions of this equation are found, one of which is exact at the outer edge of the boundary layer while the other is exact at the wall. The final solution is obtained by joining these two solutions at the inflection points of the velocity profiles. The final solution is given in terms of a series of universal functions for a fairly broad class of potential velocity distributions outside of the boundary layer. Detailed calculations of the boundary layer characteristics are worked out for the case in which the potential velocity is a linear function of the distance from the upstream stagnation point. Finally, the complete separation point characteristics are determined for the boundary layer associated with a potential velocity distribution made up of two linear functions of the distance from the stagnation point. It appears that extensions of the detailed calculations to more complex potential flows can be fairly easily carried out by using the explicit formulae given in the paper. (author)

Von Karman, TH↗

DOE OSTI2022

Charged particle motion in spherically symmetric distributions of magnetic monopoles

The classical equations of motion of a charged particle in a spherically symmetric distribution of magnetic monopoles can be transformed into a system of linear equations, thereby providing a type of integrability. In the case of a single monopole, the solution was given long ago by Poincaré. In the case of a uniform distribution of monopoles, the solution can be expressed in terms of parabolic cylinder functions (essentially the eigenfunctions of an inverted harmonic oscillator). Further, this solution is relevant to recent studies of nonassociative star products, symplectic lifts of twisted Poisson structures, and fluids and plasmas of electric and magnetic charges.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1991

Light scattering by lunar-like particle size distributions

A fundamental input to models of light scattering from planetary regoliths is the mean phase function of the regolith particles. Using the known size distribution for typical lunar soils, the mean phase function and mean linear polarization for a regolith volume element of spherical particles of any composition were calculated from Mie theory. The two contour plots given here summarize the changes in the mean phase function and linear polarization with changes in the real part of the complex index of refraction, n - ik, for k equals 0.01, the visible wavelength 0.55 micrometers, and the particle size distribution of the typical mature lunar soil 72141. A second figure is a similar index-phase surface, except with k equals 0.1. The index-phase surfaces from this survey are a first order description of scattering by lunar-like regoliths of spherical particles of arbitrary composition. They form the basis of functions that span a large range of parameter-space.

Goguen, Jay D.↗