Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “generalized linear 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 631 records · Page 35

Primordial inhomogeneities in the expanding universe. II - General features of spherical models at late times

This paper studies the density profile that forms around a spherically symmetric bound central core immersed in a homogeneous-background k = 0 or k = -1 Friedmann-Robertson-Walker cosmological model, with zero pressure. Although the density profile in the linearized regime is almost arbitrary, in the nonlinear regime certain universal features of the density profile are obtained that are independent of the details of the initial conditions. The formation of 'halos' ('holes') with densities greater than (less than) the average cosmological density is discussed. It is shown that in most regions 'halos' form, and universal values are obtained for the slope of the ln (density)-ln (radius) profile in those 'halos' at late times, independently of the shape of the initial density profile. Restrictions are derived on where it is possible for 'holes' to exist at late times and on how such 'holes' must have evolved.

Olson, D. W.↗

Geospatial modeling of near subsurface temperatures of the contiguous United States for assessment of materials degradation

Abstract Understanding subsurface temperature variations is crucial for assessing material degradation in underground structures. This study maps subsurface temperatures across the contiguous United States for depths from 50 to 3500 m, comparing linear interpolation, gradient boosting (LightGBM), neural networks, and a novel hybrid approach combining linear interpolation with LightGBM. Results reveal heterogeneous temperature patterns both horizontally and vertically. The hybrid model performed best achieving a root mean square error of 2.61 °C at shallow depths (50–350 m). Model performance generally decreased with depth, highlighting challenges in deep temperature prediction. State-level analyses emphasized the importance of considering local geological factors. This study provides valuable insights for designing efficient underground facilities and infrastructure, underscoring the need for depth-specific and region-specific modeling approaches in subsurface temperature assessment.

Science & Technology - Other Topics↗

Concepts for a theoretical and experimental study of lifting rotor random loads and vibrations. Phase 5B: Analysis of gust alleviation methods and rotor dynamic stability

The effects of various gust alleviation methods on the random blade response in flapping were studied analytically, assuming a rigid rotor support. The analytical model assumes rigid flapping blades with elastic root restraints. Linearized equations which are approximately valid at low lift conditions were used. Because of the interblade coupling from the feedback devices, the method of multiblade generalized coordinates was most convenient and was extended to include coning, differential coning, and warping of the rotor. The numerical examples cover the dynamic stability characteristics as affected by feedback gains of three- to six-bladed rotors. The number of blades had large effects on stability limits and modal time functions at these limits. The random flapping response of the blades to atmospheric turbulence was determined at 1.6 rotor advance ratio using feedback gains below the stability limit. The most effective reduction of the flapping response per unit gain was achieved with a rotor coning angle feedback.

Hohenemser, K. H.↗

Maximum likelihood identification of aircraft stability and control derivatives

Application of a generalized identification method to flight test data analysis. The method is based on the maximum likelihood (ML) criterion and includes output error and equation error methods as special cases. Both the linear and nonlinear models with and without process noise are considered. The flight test data from lateral maneuvers of HL-10 and M2/F3 lifting bodies are processed to determine the lateral stability and control derivatives, instrumentation accuracies, and biases. A comparison is made between the results of the output error method and the ML method for M2/F3 data containing gusts. It is shown that better fits to time histories are obtained by using the ML method. The nonlinear model considered corresponds to the longitudinal equations of the X-22 VTOL aircraft. The data are obtained from a computer simulation and contain both process and measurement noise. The applicability of the ML method to nonlinear models with both process and measurement noise is demonstrated.

Mehra, R. K.↗

Modelling, analyses and design of switching converters

A state-space averaging method for modelling switching dc-to-dc converters for both continuous and discontinuous conduction mode is developed. In each case the starting point is the unified state-space representation, and the end result is a complete linear circuit model, for each conduction mode, which correctly represents all essential features, namely, the input, output, and transfer properties (static dc as well as dynamic ac small-signal). While the method is generally applicable to any switching converter, it is extensively illustrated for the three common power stages (buck, boost, and buck-boost). The results for these converters are then easily tabulated owing to the fixed equivalent circuit topology of their canonical circuit model. The insights that emerge from the general state-space modelling approach lead to the design of new converter topologies through the study of generic properties of the cascade connection of basic buck and boost converters.

Cuk, S. M.↗

A multiloop generalization of the circle criterion for stability margin analysis

In order to provide a theoretical tool well suited for use in characterizing the stability margins (e.g., gain and phase margins) of multiloop feedback systems, multiloop input-output stability results generalizing the circle stability criterion are considered. Generalized conic sectors with 'centers' and 'radii' determined by linear dynamical operators are employed to enable an engineer to specify the stability margins which he desires as a frequency-dependent convex set of modeling errors (including nonlinearities, gain variations, and phase variations) which the system must be able to tolerate in each feedback loop without instability. The resulting stability criterion gives sufficient conditions for closed-loop stability in the presence of such frequency-dependent modeling errors, even when the modeling errors occur simultaneously in all loops.

Safonov, M. G.↗

The well-posedness of the Kuramoto-Sivashinsky equation

The Kuramoto-Sivashinsky equation arises in a variety of applications, among which are modeling reaction diffusion systems, flame propagation and viscous flow problems. It is considered here, as a prototype to the larger class of generalized Burgers equations: those consist of a quadratic nonlinearity and an arbitrary linear parabolic part. It is shown that such equations are well posed, thus admitting a unique smooth solution, continuously dependent on its initial data. As an attractive alternative to standard energy methods, existence and stability are derived in this case, by patching in the large short time solutions without loss of derivatives.

Tadmor, E.↗

The well-posedness of the Kuramoto-Sivashinsky equation

The Kuramoto-Sivashinsky equation arises in a variety of applications, among which are modeling reaction diffusion systems, flame propagation and viscous flow problems. It is considered here, as a prototype to the larger class of generalized Burgers equations: those consist of a quadratic nonlinearity and an arbitrary linear parabolic part. It is shown that such equations are well posed, thus admitting a unique smooth solution, continuously dependent on its initial data. As an attractive alternative to standard energy methods, existence and stability are derived in this case, by patching in the large short time solutions without 'loss of derivatives'.

Tadmor, E.↗

Deconvolution of spectra for intimate mixtures

Visible to near infrared reflectance spectra of macroscopic mixtures have been shown to be linear combinations of the reflections of the pure mineral components in the mixture. However, for microscopic mixtures the mixing systematics are in general nonlinear. The systematics may be linearized by conversion of reflectance to single scattering albedo (w), where the equations which relate reflectance to w depend on the method of data collection. Several proposed mixing models may be used to estimate mineral abundances from the reflectance spectra of intimate mixtures. These models are summarized and a revised model is presented. A noniterative (linear) least squares approach was used for curve fitting and the data, measured as bi-directional reflectance with incidence and emergence angles of 30 and 0 deg were converted to w by a simplified version of Hapke's equation for bi-directional reflectance. This model was tested with two mixture series composed of 45 to 75 micron particles: an anorthite-enstatite series and an olivine-magnetite series. The data indicate that the simplified Hapke's equation may be used to convolve reflectance spectra into mineral abundances if appropriate endmembers are known or derived from other techniques. For surfaces that contain a significant component of very low albedo material, a somewhat modified version of this technique will need to be developed. Since the abundances are calculated using a noniterative approach, the application of this method is especially efficient for large spectral data sets, such as those produced by mapping spectrometers.

Mustard, John F.↗

New parameterizations and sensitivities for simple climate models

This paper presents a reexamination of the earth radiation budget parameterization of energy balance climate models in light of data collected over the last 12 years. The study consists of three parts: (1) an examination of the infrared terrestrial radiation to space and its relationship to the surface temperature field on time scales from 1 month to 10 years; (2) an examination of the albedo of the earth with special attention to the seasonal cycle of snow and clouds; (3) solutions for the seasonal cycle using the new parameterizations with special attention to changes in sensitivity. While the infrared parameterization is not dramatically different from that used in the past, the albedo in the new data suggest that a stronger latitude dependence be employed. After retuning the diffusion coefficient the simulation results for the present climate generally show only a slight dependence on the new parameters. Also, the sensitivity parameter for the model is still about the same (1.25 C for a 1 percent increase of solar constant) for the linear models and for the nonlinear models that include a seasonal snow line albedo feedback (1.34 C). One interesting feature is that a clear-sky planet with a snow line albedo feedback has a significantly higher sensitivity (2.57 C) due to the absence of smoothing normally occurring in the presence of average cloud cover.

Graves, Charles E.↗

Coupled Attitude and Orbit Dynamics and Control in Formation Flying Systems

Formation flying systems can range from global constellations offering extended service coverage to clusters of highly coordinated vehicles that perform distributed sensing. Recently, the use of groups of micro-satellites in the areas of near Earth explorations, deep space explorations, and military applications has received considerable attention by researchers and practitioners. To date, most proposed control strategies are based on linear models (e.g., Hill-Clohessy-Wiltshire equations) or nonlinear models that are restricted to circular reference orbits. Also, all models in the literature are uncoupled between relative position and relative attitude. In this paper, a generalized dynamic model is proposed. The reference orbit is not restricted to the circular case. In this formulation, the leader or follower satellite can be in either a circular or an elliptic orbit. In addition to maintaining a specified relative position, the satellites are also required to maintain specified relative attitudes. Thus the model presented couples vehicle attitude and orbit requirements. Orbit perturbations are also included. In particular, the J(sub 2) effects are accounted in the model. Finally, a sliding mode controller is developed and used to control the relative attitude of the formation and the simulation results are presented.

Yunjun Xu↗

kL-based Linear and Nonlinear Two-Equation Turbulence Models

The development and implementation of kL-based Reynolds Average Navier-Stokes (RANS) two-equation turbulence models are reported herein. The kL is based on Abdol-Hamid's closure and Menter's modification to Rotta's two-equation model. Rotta shows that a reliable transport equation can be formed from the turbulent length scale L, and the turbulent kinetic energy k. Rotta's kL equation is well suited for term-by-term modeling and displays useful features compared to other scale formulation. One of the important differences is the inclusion of higher order velocity derivatives in the source terms of the scale equation. This can enhance the ability of RANS solvers to simulate unsteady flows in URANS mode. The present report documents the formulation of two model levels of turbulence models as implemented in the computational fluid dynamics FUN3D code. The levels are the two-equation linear k-kL and the two-equation algebraic Reynolds stress model (ARSM). Free shear, separated and corner flow cases are documented and compared with experimental, and other turbulence model data. The results show generally very good comparisons with experimental data. The results from this formulation are similar or better than results using the SST two-equation turbulence model. ARSM shows great promise with similar level of computational resources as basic two-equation turbulence models.

Abdol-Hamid, K. S.↗

Tuning Neural Network Models for Improved Prediction of Boundary Layer Transition

Boundary layer transition can strongly impact flight vehicle performance as it influences surface skin friction and aerodynamic heating, making accurate transition prediction a key to designing next generation aircraft. Artificial neural networks (ANNs) have shown promise toward predicting laminar-turbulent transition based on linear stability correlations. The computational efficiency of ANNs and the substantially reduced user involvement in relation to direct computations based on the linear stability theory (LST) makes them an attractive methodology for integrating the LST based correlations in computational fluid dynamics codes. Tollmien-Schlichting (TS) waves correspond to the dominant transition mechanism in 2D or weakly 3D subsonic boundary layers, such as those encountered in general aviation applications. Improvements to neural network model accuracy in predicting the amplification rates of TS instability waves have been investigated by leveraging recent machine learning developments in conjunction with surrogate optimization techniques and via suitable augmentation of the data used to train the networks. The optimized models trained on the modified dataset reduced the average transition location errors on different airfoils at several flow conditions by 51% of the original manually-tuned network’s errors on the same flow cases. The actual transition locations were derived from the Langley Stability and Transition Analysis Code (LASTRAC).

Machine Learning↗

Element-specific modal formulations for large-displacement multibody dynamics

Large dispacement assumed-mode modeling techniques are examined in the context of multibody elastodynamics. The range of both general and element-specific approaches are studied with the aid of examples involving beams, plates, and shells. For systems undergoing primarily structural bending and twisting with little or no membrane distortion, it is found that fully-linear, element-specific, modal formulations provide the most accurate time history solutions at the least expense. When membrane effects become dominant in structural problems due to loading and boundary conditions, one must naturally resort to a formulation involving a nonlinear stress-strain relationship in addition to nonlinear terms associated with large overall system motion. Such nonlinear models were investigated using assumed modes and found to lead to modal convergence difficulties when standard free-free structural modes are employed. A constrained mode formulation aimed at addressing the convergence problem is proposed.

Ryan, R. R.↗

Landsat test of diffuse reflectance models for aquatic suspended solids measurement

Landsat radiance data were used to test mathematical models relating diffuse reflectance to aquatic suspended solids concentration. Digital CCT data for Landsat passes over the Bay of Fundy, Nova Scotia were analyzed on a General Electric Co. Image 100 multispectral analysis system. Three data sets were studied separately and together in all combinations with and without solar angle correction. Statistical analysis and chromaticity analysis show that a nonlinear relationship between Landsat radiance and suspended solids concentration is better at curve-fitting than a linear relationship. In particular, the quasi-single-scattering diffuse reflectance model developed by Gordon and coworkers is corroborated. The Gordon model applied to 33 points of MSS 5 data combined from three dates produced r = 0.98.

Munday, J. C., Jr.↗

Predicting Flow in Fracture Networks With Quantum Algorithms

Uncertainty quantification plays a crucial role in the modeling of subsurface flow. For instance, uncertainties in the properties of geologic fracture networks significantly impact flow, requiring numerous simulations to accurately estimate quantities of interest. However, each simulation is computationally expensive because it requires solving a large linear system to capture features that involve both small and large fractures. An example is in percolation, where the interaction of many small fractures (which cumulatively can have a large surface area) with the rock matrix must be modeled precisely. Quantum computing is an emerging tool with the potential to address this issue. Quantum algorithms offer a significant speedup in solving linear systems, achieving efficiencies that are challenging to match with classical approaches. These classical approaches include direct solvers, such as LU decomposition, and iterative methods, notably preconditioned conjugate gradient, commonly used in subsurface modeling to solve large sparse systems. However, applying quantum algorithms to geologic fracture flow requires careful attention to algorithmic and problem-specific constraints to fully realize this quantum advantage. In this work we describe a quantum algorithm for generalized Monte Carlo applications with a quadratic speedup over the classical approaches which can be combined with the quantum speedup, currently under investigation, for solving quantum linear systems for subsurface flow. We show that for quantum algorithms the computational cost of estimating a quantity of interest for a statistical ensemble of networks is roughly the same as that of a single realization, essentially implying that one can get uncertainty quantification for free.

58 GEOSCIENCES↗

Analytic and numerical models of wave-CISK with conditional heating

Wave-CISK with conditional heating is investigated in the equatorial zonal-height plane by analytic and numerical techniques. For two- and three-level models, previous results are extended to give additional evidence that the most unstable mode has a single wet region of infinitesimal width. A three-level model has qualitatively similar behavior as the two-level model except that propagating solutions are possible due to coalescence of internal vertical modes. Phase speeds with conditional heating are found to be slightly greater than those for unconditional heating. The structure has one circulation cell in the vertical and is asymmetric in longitude with stronger motion on the leading edge. Growth rate is inversely proportional to the width of the single wet region. That width can be limited by second-order diffusion. A general integral relationship between growth rate, viscosity, phase speed, and heating is derived. The main conclusion is that the linear wave-CISK catastrophe is modified by conditional heating but not eliminated. The preferred mode of instability has one wet region, but it occurs on the smallest possible scale. It is likely that numerical models that use conditional heating are sensitive to resolution, especially for the commonly used spectral truncations, unless there is sufficiently strong damping at the smallest scales.

Crum, Francis X.↗

Towards a Comprehensive Model of Jet Noise Using an Acoustic Analogy and Steady RANS Solutions

An acoustic analogy is developed to predict the noise from jet flows. It contains two source models that independently predict the noise from turbulence and shock wave shear layer interactions. The acoustic analogy is based on the Euler equations and separates the sources from propagation. Propagation effects are taken into account by calculating the vector Green's function of the linearized Euler equations. The sources are modeled following the work of Tam and Auriault, Morris and Boluriaan, and Morris and Miller. A statistical model of the two-point cross-correlation of the velocity fluctuations is used to describe the turbulence. The acoustic analogy attempts to take into account the correct scaling of the sources for a wide range of nozzle pressure and temperature ratios. It does not make assumptions regarding fine- or large-scale turbulent noise sources, self- or shear-noise, or convective amplification. The acoustic analogy is partially informed by three-dimensional steady Reynolds-Averaged Navier-Stokes solutions that include the nozzle geometry. The predictions are compared with experiments of jets operating subsonically through supersonically and at unheated and heated temperatures. Predictions generally capture the scaling of both mixing noise and BBSAN for the conditions examined, but some discrepancies remain that are due to the accuracy of the steady RANS turbulence model closure, the equivalent sources, and the use of a simplified vector Green's function solver of the linearized Euler equations.

Miller, Steven A. E.↗