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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Using system identification in modeling the yaw response of tail fins for small wind turbines with bearing friction

Here, we describe three main developments of our previous study of the nonlinear yaw dynamics of the tail fins for small wind turbines [Khedr et al., J. Renewable Sustainable Energy 16, 053305 (2024)]. First, the model constants derived from archived computational and experimental studies are adjusted by employing system identification (SI) to maximize the model's agreement with wind tunnel tests. This adjustment was done for high wind speeds, where yaw bearing friction can be ignored. When starting a turbine at low wind speed, however, friction can become important. Our second development is to implement a model for the frictional resistive torque and use SI to maximize its accuracy. These developments used wind tunnel experiments on generic delta, elliptical, and rectangular planforms that were described by Khedr et al. [J. Renewable Sustainable Energy 16, 053305 (2024)]. Since the aerodynamic and friction models employ a large number of constants, we describe ways to constrain the values using linearized solutions of the response equations for small and large yaw angles. Third, we test the generality of the aerodynamic and friction modeling using the complex planform from a commercial small turbine, for which limited theoretical and computational guidance is available in selecting the model constants. The model implemented with SI is shown to provide an accurate description of the yaw response of the complex planform. Guidelines are given for the use of wind tunnel tests to determine the model constants for tail fins of any planform.

17 WIND ENERGY↗
DOE OSTI2024

FARM User Guidance and Instructions

The FARM (Feasible Actuator Range Modifier) module is a component of the RAVEN-based FORCE framework for the analysis of Integrated Energy Systems (IES). FARM aids HERON in the solution of the power dispatch problem by evaluating feasible set-point signals to be issued to the control systems of the different IES unit components. Set-points need to satisfy limits on both production variables (i.e., the variables to be optimized such as the electrical power, etc.) and process variables tied to the service life of equipment (e.g., steam flowrate, vessel pressure, turbine firing temperature, etc.). To enforce all these limits, a two-stage approach is adopted. First, the power dispatcher algorithm in HERON module estimates set-points that meet the constraints on the production variables, e.g., power levels and power ramp rate limits. These constraints are called explicit constraints. Then, if necessary, FARM adjusts these set-points to ensure the respect of the limits on the process variables of interest, given the knowledge of the system dynamics acquired through machine learning algorithms. These constraints are called implicit constraints. From this standpoint, FARM constitutes a bridge between the HERON power dispatcher that adopts a simplified description of the IES unit (low-resolution physics) and the HYBRID high-fidelity models (high-resolution physics). In this report, the general workflow and the software requirements of FARM module are summarized, and the detailed instructions for installing FARM software, running built-in example cases, deriving Linear Parameter-Varying (LPV) state-space models, and using FARM for user-defined power dispatch problems are provided.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Linear Time-Invariant Models of a Large Cumulus Ensemble

Abstract Methods in system identification are used to obtain linear time-invariant state-space models that describe how horizontal averages of temperature and humidity of a large cumulus ensemble evolve with time under small forcing. The cumulus ensemble studied here is simulated with cloud-system-resolving models in radiative–convective equilibrium. The identified models extend steady-state linear response functions used in past studies and provide accurate descriptions of the transfer function, the noise model, and the behavior of cumulus convection when coupled with two-dimensional gravity waves. A novel procedure is developed to convert the state-space models into an interpretable form, which is used to elucidate and quantify memory in cumulus convection. The linear problem studied here serves as a useful reference point for more general efforts to obtain data-driven and interpretable parameterizations of cumulus convection.

Meteorology & Atmospheric Sciences↗
NASA NTRS1985

Variability simulations with a steady, linearized primitive equations model

Solutions of the steady, primitive equations on a sphere, linearized about a zonally symmetric basic state are computed for the purpose of simulating monthly mean variability in the troposphere. The basic states are observed, winter monthly mean, zonal means of zontal and meridional velocities, temperatures and surface pressures computed from the 15 year NMC time series. A least squares fit to a series of Legendre polynomials is used to compute the basic states between 20 H and the equator, and the hemispheres are assumed symmetric. The model is spectral in the zonal direction, and centered differences are employed in the meridional and vertical directions. Since the model is steady and linear, the solution is obtained by inversion of a block, pente-diagonal matrix. The model simulates the climatology of the GFDL nine level, spectral general circulation model quite closely, particularly in middle latitudes above the boundary layer. This experiment is an extension of that simulation to examine variability of the steady, linear solution.

Kinter, J. L., III↗
DOE OSTI2025

Next-Level Energy Management in Manufacturing: Facility-Level Energy Digital Twin Framework Based on Machine Learning and Automated Data Collection

This research introduces an energy prediction framework at the facility level supported by automated data collection and machine learning models. It investigates whether reducing the prediction time scale allows for applying more complex machine learning techniques and if those techniques improve the prediction accuracy. The primary advantages of this framework lie in its automation of the energy prediction process and its provision of real-time energy data suitable for use in energy dashboards or digital twins. A sitewide dataset was created by combining 15 min energy and daily production data of five shops—assembly, battery, body (electric), body (gas), and paint—from a globally recognized electric vehicle manufacturer. Various machine learning models were evaluated on daily, weekly, and monthly datasets, including, in increasingly complex order: naïve, simple linear regression, net regularized generalized linear regression, principal component regression, k-nearest neighbor, random forest, and Bayesian regularized neural network. Compared to the current state-of-the-art energy consumption prediction for the industrial facility level, this research investigates more complex models and smaller time intervals for higher accuracy. The findings revealed that the more complex monthly models require a minimum of a year and a half of data to operate, while weekly models demand a year of data to achieve improved accuracy. Daily models can operate with only six months of data but exhibit poor performance due to reduced prediction accuracy of production. Key challenges identified include access to reliable, high-quality energy and production data and the initial demand for human labor.

digital twin↗
DOE OSTI2021

Topographical shifts in fine root lifespan in a mixed, mesic temperate forest

Root lifespan, often is estimated in landscape- and ecosystem-level carbon models using linear approximations. In water manipulation experiments, fine root lifespan can vary with soil water content. Soil water content is generally structured by complex topography, which is largely unaccounted for in landscape- and ecosystem-scale carbon models. Topography governs the range of soil water content experienced by roots which may impact their lifespan. We hypothesized that root lifespan varied nonlinearly across a temperate, mesic, forested catchment due to differences in soil water content associated with topographic position. We expected regions of the landscape that were too wet or too dry would have soils that were not optimal for roots and thus result in shorter root lifespans. Specifically, we hypothesized that root lifespan would be longest in areas that consistently had soil water content in the middle of the soil water content spectrum, while in soils at either very low or very high soil water content, root lifespan would be relatively short. We tested this hypothesis by collecting and analyzing two years of minirhizotron and soil moisture data in plots widely distributed in the Shale Hills catchment of the Susquehanna-Shale Hills Critical Zone Observatory in Pennsylvania. We found that fine root lifespans were longer in traditionally wetter topographic regions, but detected no short term (biweekly) effect of soil moisture on root lifespan. Additionally, depth in soil, soil series, slope face orientation, and season of birth strongly affected root lifespans across the catchment. In contrast, lifespan was unaffected by root diameter or mycorrhizal association. Failure to account for these variables could result in erroneous estimates of fine root lifespan and, consequentially, carbon flux in temperate forested regions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Nature of size effects in compact models of field effect transistors

Investigations have shown that in the local approximation (for sizes L < 100 μm), AlGaN/GaN high electron mobility transistor (HEMT) structures satisfy to all properties of chaotic systems and can be described in the language of fractal geometry of fractional dimensions. For such objects, values of their electrophysical characteristics depend on the linear sizes of the examined regions, which explain the presence of the so-called size effects—dependences of the electrophysical and instrumental characteristics on the linear sizes of the active elements of semiconductor devices. In the present work, a relationship has been established for the linear model parameters of the equivalent circuit elements of internal transistors with fractal geometry of the heteroepitaxial structure manifested through a dependence of its relative electrophysical characteristics on the linear sizes of the examined surface areas. For the HEMTs, this implies dependences of their relative static (A/mm, mA/V/mm, Ω/mm, etc.) and microwave characteristics (W/mm) on the width d of the sink-source channel and on the number of sections n that leads to a nonlinear dependence of the retrieved parameter values of equivalent circuit elements of linear internal transistor models on n and d. Thus, it has been demonstrated that the size effects in semiconductors determined by the fractal geometry must be taken into account when investigating the properties of semiconductor objects on the levels less than the local approximation limit and designing and manufacturing field effect transistors. In general, the suggested approach allows a complex of problems to be solved on designing, optimizing, and retrieving the parameters of equivalent circuits of linear and nonlinear models of not only field effect transistors but also any arbitrary semiconductor devices with nonlinear instrumental characteristics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Characterization and Modeling of Polymeric Foam Under Multi-Axial Static and Dynamic Loading

A polymeric foam commonly used in composite sandwich structures was characterized under multi-axial loading at strain rates varying from quasi-static to dynamic. Additionally, tests were conducted under uniaxial compression, tension, pure shear and combinations of normal and shear stresses. Quasi-static and intermediate strain rate tests were conducted in a servo-hydraulic testing machine. High strain rate tests were conducted using a split Hopkinson pressure bar (Kolsky bar) system made of polycarbonate bars having an impedance compatible to that of the foam material. The typical compressive stress-strain behavior of the polymeric foam exhibits a linear elastic region up to a yield point, a nonlinear elastic-plastic region up to an initial peak or “critical stress” corresponding to collapse initiation of the cells, followed by strain softening up to a local minimum (plateau or saddle point stress) and finally, a strain hardening region up to densification of the foam. The characteristic stresses of the stress-strain behavior vary linearly with the logarithm of strain rate. A general three-dimensional elastic-viscoplastic model, formulated in strain space, was proposed. The model expresses the multi-axial state of stress in terms of an effective stress, incorporates strain rate effects and includes the large deformation region. Stress-strain curves obtained under multi-axial loading at different strain rates were used to develop and validate the elastic-viscoplastic constitutive model. Lastly, excellent agreement was shown between model predictions and experimental results.

36 MATERIALS SCIENCE↗
DOE OSTI2010

GENERALIZATION OF THE LAST SCATTERING APPROXIMATION FOR THE SECOND SOLAR SPECTRUM MODELING: THE Ca I 4227 A LINE AS A CASE STUDY

To model the second solar spectrum (the linearly polarized spectrum of the Sun that is due to coherent scattering processes), one needs to solve the polarized radiative transfer (RT) equation. For strong resonance lines, partial frequency redistribution (PRD) effects must be accounted for, which make the problem computationally demanding. The 'last scattering approximation' (LSA) is a concept that has been introduced to make this highly complex problem more tractable. An earlier application of a simple LSA version could successfully model the wings of the strong Ca I 4227 A resonance line in Stokes Q/I (fractional linear polarization), but completely failed to reproduce the observed Q/I peak in the line core. Since the magnetic field signatures from the Hanle effect only occur in the line core, we need to generalize the existing LSA approach if it is to be useful for the diagnostics of chromospheric and turbulent magnetic fields. In this paper, we explore three different approximation levels for LSA and compare each of them with the benchmark represented by the solution of the full polarized RT, including PRD effects. The simplest approximation level is LSA-1, which uses the observed center-to-limb variation of the intensity profile to obtain the anisotropy of the radiation field at the surface, without solving any transfer equation. In contrast, the next two approximation levels use the solution of the unpolarized transfer equation to derive the anisotropy of the incident radiation field and use it as an input. In the case of LSA-2, the anisotropy at level {tau}{sub {lambda}} = {mu}, the atmospheric level from which an observed photon is most likely to originate, is used. LSA-3, on the other hand, makes use of the full depth dependence of the radiation anisotropy. The Q/I formula for LSA-3 is obtained by keeping the first term in a series expansion of the Q-source function in powers of the mean number of scattering events. Computationally, LSA-1 is 21 times faster than LSA-2, which is 5 times faster than the more general LSA-3, which itself is 8 times faster than the polarized RT approach. A comparison of the calculated Q/I spectra with the RT benchmark shows excellent agreement for LSA-3, including good modeling of the Q/I core region with its PRD effects. In contrast, both LSA-1 and LSA-2 fail to model the core region. The RT and LSA-3 approaches are then applied to model the recently observed Q/I profile of the Ca I 4227 A line in quiet regions of the Sun. Apart from a global scale factor both give a very good fit to the Q/I spectra for all the wavelengths, including the core peak and blend line depolarizations. We conclude that LSA-3 is an excellent substitute for the full polarized RT and can be used to interpret the second solar spectrum, including the Hanle effect with PRD. It also allows the techniques developed for unpolarized three-dimensional RT to be applied to the modeling of the second solar spectrum.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

A note on the S-matrix bootstrap for the 2d O(N) bosonic model

In this work we apply the S-matrix bootstrap maximization program to the 2d bosonic O ( N ) integrable model which has N species of scalar particles of mass m and no bound states. Since in previous studies theories were defined by maximizing the coupling between particles and their bound states, the main problem appears to be to find what other functional can be used to define this model. Instead, we argue that the defining property of this integrable model is that it resides at a vertex of the convex space determined by the unitarity and crossing constraints. Thus, the integrable model can be found by maximizing any linear functional whose gradient points in the general direction of the vertex, namely within a cone determined by the normals to the faces intersecting at the vertex. This is a standard problem in applied mathematics, related to semi-definite programming and solvable by fast available numerical algorithms. The information provided by the numerical solution is enough to reproduce the known analytical solution without using integrability, namely the Yang-Baxter equation. This situation seems quite generic so we expect that other theories without continuous parameters can also be found by maximizing linear functionals in the convex space of allowed S-matrices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Inflection-point inflation in a hyper-charge oriented U ( 1 ) X model

Inflection-point inflation is an interesting possibility to realize a successful slow-roll inflation when inflation is driven by a single scalar field with its value during inflation below the Planck mass (ΦI≲M Pl ). In order for a renormalization group (RG) improved effective λΦ 4 potential to develop an inflection-point, the running quartic coupling λ(Φ) must exhibit a minimum with an almost vanishing value in its RG evolution, namely λ(Φ I )≃0 and β λ (ΦI)≃0, where β λ is the beta-function of the quartic coupling. Here in this paper, we consider the inflection-point inflation in the context of the minimal gauged U(1) X extended Standard Model (SM), which is a generalization of the minimal U(1) B$-$L model, and is constructed as a linear combination of the SM U(1) Y and U(1) B$-$L gauge symmetries. We identify the U(1) X Higgs field with the inflaton field. For a successful inflection-point inflation to be consistent with the current cosmological observations, the mass ratios among the U(1) X gauge boson, the right-handed neutrinos and the U(1) X Higgs boson are fixed. Focusing on the case that the extra U(1) X gauge symmetry is mostly aligned along the SM U(1) Y direction, we investigate a consistency between the inflationary predictions and the latest LHC Run-2 results on the search for a narrow resonance with the di-lepton final state.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2011

Coherent Vortices in Strongly Coupled Liquids

Strongly coupled liquids are ubiquitous in both nature and laboratory plasma experiments. They are unique in the sense that their average potential energy per particle dominates over the average kinetic energy. Using ''first principles'' molecular dynamics (MD) simulations, we report for the first time the emergence of isolated coherent tripolar vortices from the evolution of axisymmetric flows in a prototype two-dimensional (2D) strongly coupled liquid, namely, the Yukawa liquid. Linear growth rates directly obtained from MD simulations are compared with a generalized hydrodynamic model. Our MD simulations reveal that the tripolar vortices persist over several turn over times and hence may be observed in strongly coupled liquids such as complex plasma, liquid metals and astrophysical systems such as white dwarfs and giant planetary interiors, thereby making the phenomenon universal.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Multiresonant Feedback Control of a Three-Degree-of-Freedom Wave Energy Converter

For a three-degree-of-freedom wave energy converter (heave, pitch, and surge), the equations of motion could be coupled depending on the buoy shape. This paper presents a multiresonant feedback control, in a general framework, for this type of a wave energy converter that is modeled by linear time invariant dynamic systems. The proposed control strategy finds the optimal control in the sense that it computes the control based on the complex conjugate criteria. This control strategy is relatively easy to implement since it is a feedback control in the time domain that requires only measurements of the buoy motion. Numerical tests are presented for two different buoy shapes: a sphere and a cylinder. Regular, Bretschnieder, and Ochi-Hubble waves are tested. Simulation results show that the proposed controller harvests energy in the pitch-surge-heave modes that is about three times the energy that can be harvested using a heave-only device. In conclusion, this multiresonant control can also be used to shift the energy harvesting between the coupled modes, which can be exploited to eliminate one of the actuators while maintaining about the same level of energy harvesting.

42 ENGINEERING↗
DOE OSTI2019

The allocation of a variable series reactor considering AC constraints and contingencies

The Variable Series Reactors (VSRs) can efficiently control the power flow through the adjustment of the line reactance. When they are appropriately allocated in the power network, the transmission congestion and generation cost can be reduced. Our work proposes a planning model to optimally allocate VSRs considering AC constraints and multi-scenarios including base case and contingencies. The planning model is originally a non-convex large scale mixed integer nonlinear program (MINLP), which is generally intractable. The proposed Benders approach decomposes the MINLP model into a mixed integer linear program (MILP) master problem and a number of nonlinear subproblems. Numerical case studies based on IEEE 118-bus demonstrate the high performance of the proposed approach.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2011

A Brief Review of Viscosity Models for Slag in Coal Gasification

Many researchers have defined the phenomenon of 'slagging' as the deposition of ash in the radiative section of a boiler, while 'fouling' refers to the deposition of ash in the convective-pass region. Among the important parameters affecting ash deposition that need to be studied are ash chemistry, its transport, deposit growth, and strength development; removability of the ash deposit; heat transfer mechanisms; and the mode of operation for boilers. The heat transfer at the walls of a combustor depends on many parameters including ash deposition. This depends on the processes or parameters controlling the impact efficiency and the sticking efficiency. For a slagging combustor or furnace, however, the temperatures are so high that much of the coal particles are melted and the molten layer, in turn, captures more particles as it flows. The main problems with ash deposition are reduced heat transfer in the boiler and corrosion of the tubes. Common ways of dealing with these issues are soot blowing and wall blowing on a routine basis; however, unexpected or uncontrolled depositions can also complicate the situation, and there are always locations inaccessible to the use of such techniques. Studies have indicated that slag viscosity must be within a certain range of temperatures for tapping and the membrane wall to be accessible, for example, between 1300 C and 1500 C, the viscosity is approximately 25 Pa {center_dot} s. As the operating temperature decreases, the slag cools and solid crystals begin to form. In such cases the slag should be regarded as a non-Newtonian suspension, consisting of liquid silicate and crystals. A better understanding of the rheological properties of the slag, such as yield stress and shear-thinning, are critical in determining the optimum operating conditions. To develop an accurate heat transfer model in any type of coal combustion or gasification process, the heat transfer and to some extent the rheological properties of ash and slag, especially in high-temperature environments need to be understood and properly modeled. The viscosity of slag and the thermal conductivity of ash deposits are among two of the most important constitutive parameters that need to be studied. The accurate formulation or representations of the (transport) properties of coal (and biomass for co-firing cases) present a special challenge of modeling efforts in computational fluid dynamics applications. In this report, we first provide a brief review of the various approaches taken by different researchers in formulating or obtaining a slag viscosity model. In general, these models are based on experiments. Since slag behaves as a non-linear fluid, we discuss the constitutive modeling of slag and the important parameters that must be studied.

01 COAL, LIGNITE, AND PEAT↗
NASA NTRS1968

Survey of fracture toughness test methods

Comprehensive survey presents current methods of fracture toughness testing that are based on linear elastic fracture mechanics. General principles of the basic two dimensional crack stress field model are discussed in relation to real three dimensional specimens. Methods of test instrumentation and procedure are described.

Brown, W. F., Jr.↗
NASA NTRS1979

Identification evaluation methods

Methods for airplane parameter estimation, the equation error method, output error method, and two advanced methods are presented and their basic properties described. The advanced methods include the maximum likelihood and extended Kalman filter method. For a better understanding of the estimation techniques a first-order scalar differential equation is used as a model of the system under test. Application of the methods to a general multivariable linear system is briefly outlined. A note on the parameter estimation in the frequency domain is also presented. Numerical examples along with the comparison of results from various methods are given.

Klein, V.↗
DOE OSTI2015

Stable static structures in models with higher-order derivatives

We investigate the presence of static solutions in generalized models described by a real scalar field in four-dimensional space–time. We study models in which the scalar field engenders higher-order derivatives and spontaneous symmetry breaking, inducing the presence of domain walls. Despite the presence of higher-order derivatives, the models keep to equations of motion second-order differential equations, so we focus on the presence of first-order equations that help us to obtain analytical solutions and investigate linear stability on general grounds. We then illustrate the general results with some specific examples, showing that the domain wall may become compact and that the zero mode may split. Moreover, if the model is further generalized to include k-field behavior, it may contribute to split the static structure itself.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗