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 667 records · Page 37

Verification of a Byzantine-Fault-Tolerant Self-stabilizing Protocol for Clock Synchronization

This paper presents the mechanical verification of a simplified model of a rapid Byzantine-fault-tolerant self-stabilizing protocol for distributed clock synchronization systems. This protocol does not rely on any assumptions about the initial state of the system except for the presence of sufficient good nodes, thus making the weakest possible assumptions and producing the strongest results. This protocol tolerates bursts of transient failures, and deterministically converges within a time bound that is a linear function of the self-stabilization period. A simplified model of the protocol is verified using the Symbolic Model Verifier (SMV). The system under study consists of 4 nodes, where at most one of the nodes is assumed to be Byzantine faulty. The model checking effort is focused on verifying correctness of the simplified model of the protocol in the presence of a permanent Byzantine fault as well as confirmation of claims of determinism and linear convergence with respect to the self-stabilization period. Although model checking results of the simplified model of the protocol confirm the theoretical predictions, these results do not necessarily confirm that the protocol solves the general case of this problem. Modeling challenges of the protocol and the system are addressed. A number of abstractions are utilized in order to reduce the state space.

Malekpour, Mahyar R.↗

Generalized Predictive Control of Dynamic Systems with Rigid-Body Modes

Numerical simulations to assess the effectiveness of Generalized Predictive Control (GPC) for active control of dynamic systems having rigid-body modes are presented. GPC is a linear, time-invariant, multi-input/multi-output predictive control method that uses an ARX model to characterize the system and to design the controller. Although the method can accommodate both embedded (implicit) and explicit feedforward paths for incorporation of disturbance effects, only the case of embedded feedforward in which the disturbances are assumed to be unknown is considered here. Results from numerical simulations using mathematical models of both a free-free three-degree-of-freedom mass-spring-dashpot system and the XV-15 tiltrotor research aircraft are presented. In regulation mode operation, which calls for zero system response in the presence of disturbances, the simulations showed reductions of nearly 100%. In tracking mode operations, where the system is commanded to follow a specified path, the GPC controllers produced the desired responses, even in the presence of disturbances.

Kvaternik, Raymond G.↗

Radiative transfer calculated from a Markov chain formalism

The theory of Markov chains is used to formulate the radiative transport problem in a general way by modeling the successive interactions of a photon as a stochastic process. Under the minimal requirement that the stochastic process is a Markov chain, the determination of the diffuse reflection or transmission from a scattering atmosphere is equivalent to the solution of a system of linear equations. This treatment is mathematically equivalent to, and thus has many of the advantages of, Monte Carlo methods, but can be considerably more rapid than Monte Carlo algorithms for numerical calculations in particular applications. We have verified the speed and accuracy of this formalism for the standard problem of finding the intensity of scattered light from a homogeneous plane-parallel atmosphere with an arbitrary phase function for scattering. Accurate results over a wide range of parameters were obtained with computation times comparable to those of a standard 'doubling' routine. The generality of this formalism thus allows fast, direct solutions to problems that were previously soluble only by Monte Carlo methods. Some comparisons are made with respect to integral equation methods.

Esposito, L. W.↗

Revisiting matrix-based inversion of scanning mobility particle sizer (SMPS) and humidified tandem differential mobility analyzer (HTDMA) data

Abstract. Tikhonov regularization is a tool for reducing noise amplification during data inversion. This work introduces RegularizationTools.jl, a general-purpose software package for applying Tikhonov regularization to data. The package implements well-established numerical algorithms and is suitable for systems of up to ∼ 1000 equations. Included is an abstraction to systematically categorize specific inversion configurations and their associated hyperparameters. A generic interface translates arbitrary linear forward models defined by a computer function into the corresponding design matrix. This obviates the need to explicitly write out and discretize the Fredholm integral equation, thus facilitating fast prototyping of new regularization schemes associated with measurement techniques. Example applications include the inversion involving data from scanning mobility particle sizers (SMPSs) and humidified tandem differential mobility analyzers (HTDMAs). Inversion of SMPS size distributions reported in this work builds upon the freely available software DifferentialMobilityAnalyzers.jl. The speed of inversion is improved by a factor of ∼ 200, now requiring between 2 and 5 ms per SMPS scan when using 120 size bins. Previously reported occasional failure to converge to a valid solution is reduced by switching from the L-curve method to generalized cross-validation as the metric to search for the optimal regularization parameter. Higher-order inversions resulting in smooth, denoised reconstructions of size distributions are now included in DifferentialMobilityAnalyzers.jl. This work also demonstrates that an SMPS-style matrix-based inversion can be applied to find the growth factor frequency distribution from raw HTDMA data while also accounting for multiply charged particles. The outcome of the aerosol-related inversion methods is showcased by inverting multi-week SMPS and HTDMA datasets from ground-based observations, including SMPS data obtained at Bodega Marine Laboratory during the CalWater 2/ACAPEX campaign and co-located SMPS and HTDMA data collected at the US Department of Energy observatory located at the Southern Great Plains site in Oklahoma, USA. Results show that the proposed approaches are suitable for unsupervised, nonparametric inversion of large-scale datasets as well as inversion in real time during data acquisition on low-cost reduced-instruction-set architectures used in single-board computers. The included software implementation of Tikhonov regularization is freely available, general, and domain-independent and thus can be applied to many other inverse problems arising in atmospheric measurement techniques and beyond.

54 ENVIRONMENTAL SCIENCES↗

Collisional excitation of interstellar molecules - Linear molecules CO, CS, OCS, and HC3N

Rates for excitation of CO, CS, OCS, and HC3N by collisions with He atoms and H2 molecules are presented. These have been obtained from extensive theoretical calculations using methods of known reliability. In general, rates necessary for modeling interstellar molecular clouds with kinetic temperatures to 100 K have been included. Methods for fitting these rates in more compact form are discussed.

Green, S.↗

The Total Synthesis Problem of linear multivariable control. II - Unity feedback and the design morphism

Zames (1981) has observed that there is, in general, no 'separation principle' to guarantee optimality of a division between control law design and filtering of plant uncertainty. Peczkowski and Sain (1978) have solved a model matching problem using transfer functions. Taking into consideration this investigation, Peczkowski et al. (1979) proposed the Total Synthesis Problem (TSP), wherein both the command/output-response and command/control-response are to be synthesized, subject to the plant constraint. The TSP concept can be subdivided into a Nominal Design Problem (NDP), which is not dependent upon specific controller structures, and a Feedback Synthesis Problem (FSP), which is. Gejji (1980) found that NDP was characterized in terms of the plant structural matrices and a single, 'good' transfer function matrix. Sain et al. (1981) have extended this NDP work. The present investigation is concerned with a study of FSP for the unity feedback case. NDP, together with feedback synthesis, is understood as a Total Synthesis Problem.

Sain, M. K.↗

PROTEUS two-dimensional Navier-Stokes computer code, version 1.0. Volume 1: Analysis description

A new computer code was developed to solve the two-dimensional or axisymmetric, Reynolds averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The thin-layer or Euler equations may also be solved. Turbulence is modeled using an algebraic eddy viscosity model. The objective was to develop a code for aerospace applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The equations are written in nonorthogonal body-fitted coordinates, and solved by marching in time using a fully-coupled alternating direction-implicit procedure with generalized first- or second-order time differencing. All terms are linearized using second-order Taylor series. The boundary conditions are treated implicitly, and may be steady, unsteady, or spatially periodic. Simple Cartesian or polar grids may be generated internally by the program. More complex geometries require an externally generated computational coordinate system. The documentation is divided into three volumes. Volume 1 is the Analysis Description, and describes in detail the governing equations, the turbulence model, the linearization of the equations and boundary conditions, the time and space differencing formulas, the ADI solution procedure, and the artificial viscosity models.

Towne, Charles E.↗

A Coupled Aeroelastic Model for Launch Vehicle Stability Analysis

A technique for incorporating distributed aerodynamic normal forces and aeroelastic coupling effects into a stability analysis model of a launch vehicle is presented. The formulation augments the linear state-space launch vehicle plant dynamics that are compactly derived as a system of coupled linear differential equations representing small angular and translational perturbations of the rigid body, nozzle, and sloshing propellant coupled with normal vibration of a set of orthogonal modes. The interaction of generalized forces due to aeroelastic coupling and thrust can be expressed as a set of augmenting non-diagonal stiffness and damping matrices in modal coordinates with no penalty on system order. While the eigenvalues of the structural response in the presence of thrust and aeroelastic forcing can be predicted at a given flight condition independent of the remaining degrees of freedom, the coupled model provides confidence in closed-loop stability in the presence of rigid-body, slosh, and actuator dynamics. Simulation results are presented that characterize the coupled dynamic response of the Ares I launch vehicle and the impact of aeroelasticity on control system stability margins.

Orr, Jeb S.↗

Lattice strains and diffraction elastic constants of cubic polycrystals

In situ synchrotron X-ray and neutron diffraction experiments provide a powerful approach to measure lattice strains in bulk polycrystalline materials. They are being increasingly used for quantitative characterization of microscale deformation within and between grains and phases. Here we use a self-consistent micromechanics model to obtain a general analytic solution of the grain-level lattice strains and diffraction elastic constants for a broad class of elastically isotropic polycrystals with cubic crystal symmetry. This analytic solution reveals a direct linear relationship between the reciprocal of the elastic diffraction constant and the orientation index parameter along the direction of any diffraction vector, including tensile loading and transverse directions. The straightforward numerical implementation of this solution provides diffraction elastic constants for 26 representative cubic polycrystals. Finally, analytic solutions of this kind can serve to benchmark in situ diffraction measurements of lattice strains and also facilitate high-throughput studies of microscale stresses and diffraction elastic constants in polycrystalline materials.

36 MATERIALS SCIENCE↗

Drag reduction and the Vogel exponent of a flexible beam in transient shear flows

Interactions between a flexible beam and a fluid in a channel are of great relevance to biological hairy surfaces, aquatic vegetation, marine life (e.g., fish gills), and many industrial systems alike. While steady state response of a beam to such flows is fairly well-explored, their behavior in the transient regime is not fully understood. Here, a series of numerical simulations are performed to study the laminar Couette flow of an incompressible viscous fluid past an elastic beam in a two-dimensional channel. The flexible beam is perpendicular to the direction of flow, and its base is fixed to the stationary bottom of the channel. We measure the evolution of the Vogel exponent, drag reduction, and reconfiguration number during the transient and steady-state response of the fluid–structure system for different geometrical and physical properties. Our benchmark shows a good agreement between numerical and experimental observations. Our results show that the system's steady-state response at different bulk-fluid velocities can be reproduced by investigating the shear flow response during the transient regime. We define a new variable that characterizes the evolution of the local velocity profile in the proximity of the free end of the beam and use that to characterize the transient-regime response. The analysis yields insight into the competing effects of elasticity of the beam and non-linear flow response.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Multiscale Approach to Fast ModSim for Laser Processing of Metals for Future Nuclear Deterrence Environments

Predicting performance of parts produced using laser-metal processing remains an out- standing challenge. While many computational models exist, they are generally too computationally expensive to simulate the build of an engineering-scale part. This work develops a reduced order thermal model of a laser-metal system using analytical Green's function solutions to the linear heat equation, representing a step towards achieving a full part performance prediction in an "overnight" time frame. The developed model is able to calculate a thermal history for an example problem 72 times faster than a traditional FEM method. The model parameters are calibrated using a non-linear solution and microstructures and residual stresses calculated and compared to a non-linear case. The calibrated model shows promising agreement with a non-linear solution.

36 MATERIALS SCIENCE↗

On the synthesis of multivariable systems.

Description of a general synthesis procedure for the compensation of linear multivariable systems through the combined use of dynamic feed-forward compensation and linear state variable feedback. Applications of the synthesis algorithm presented to problems of decoupling and exact model matching illustrate its use.

Wolovich, W. A.↗

A sequential linear optimization approach for controller design

A linear optimization approach with a simple real arithmetic algorithm is presented for reliable controller design and vibration suppression of flexible structures. Using first order sensitivity of the system eigenvalues with respect to the design parameters in conjunction with a continuation procedure, the method converts a nonlinear optimization problem into a maximization problem with linear inequality constraints. The method of linear programming is then applied to solve the converted linear optimization problem. The general efficiency of the linear programming approach allows the method to handle structural optimization problems with a large number of inequality constraints on the design vector. The method is demonstrated using a truss beam finite element model for the optimal sizing and placement of active/passive-structural members for damping augmentation. Results using both the sequential linear optimization approach and nonlinear optimization are presented and compared. The insensitivity to initial conditions of the linear optimization approach is also demonstrated.

Horta, L. G.↗

Passive stabilization for large space systems

The optimal tuning of multiple tuned-mass dampers for the transient vibration damping of large space structures is investigated. A multidisciplinary approach is used. Structural dynamic techniques are applied to gain physical insight into absorber/structure interaction and to optimize specific cases. Modern control theory and parameter optimization techniques are applied to the general optimization problem. A design procedure for multi-absorber multi-DOF vibration damping problems is presented. Classical dynamic models are extended to investigate the effects of absorber placement, existing structural damping, and absorber cross-coupling on the optimal design synthesis. The control design process for the general optimization problem is formulated as a linear output feedback control problem via the development of a feedback control canonical form. The techniques are applied to sample micro-g and pointing problems on the NASA dual keel space station.

Sesak, J. R.↗

Probabilistic boundary element method

The purpose of the Probabilistic Structural Analysis Method (PSAM) project is to develop structural analysis capabilities for the design analysis of advanced space propulsion system hardware. The boundary element method (BEM) is used as the basis of the Probabilistic Advanced Analysis Methods (PADAM) which is discussed. The probabilistic BEM code (PBEM) is used to obtain the structural response and sensitivity results to a set of random variables. As such, PBEM performs analogous to other structural analysis codes such as finite elements in the PSAM system. For linear problems, unlike the finite element method (FEM), the BEM governing equations are written at the boundary of the body only, thus, the method eliminates the need to model the volume of the body. However, for general body force problems, a direct condensation of the governing equations to the boundary of the body is not possible and therefore volume modeling is generally required.

Cruse, T. A.↗

Progress in integrated analysis with adaptive unstructured meshing

Design of lightweight structures and thermal protection systems for hypersonic vehicles depend on accurate prediction of aerothermal loads, structural temperatures and their gradients, and structural deformations and stresses. Concentration is on an alternative meshing technique which generates an entirely new adaptive unstructured mesh based on the solution obtained from the earlier mesh. The technique combined with the finite element method has been shown to significantly improve the efficiency and accuracy of the fluid, thermal, and structural analyses. Current capability of the adaptive unstructured meshing technique for the integrated fluid-thermal-structural analysis is described first. The technique was extended to transient thermal analysis of structures with time-dependent adaptive meshing to capture the detailed temperature response with a minimum number of unknowns and computational cost. Both linear and higher-order finite elements are implemented to demonstrate the generality of the technique and to investigate their solution accuracy. Currently, the adaptive meshing technique is being developed for plane structures that can be modeled with membrane elements and built-up structures modeled with membrane and bending elements. The capability of the technique to these different disciplinary problems is demonstrated by several examples.

Dechaumphai, Pramote↗

Pressure Distributions and Wave Drag Due to Two-Dimensional Fabrication-Type Surface Roughness on an Ogive Cylinder at Mach Numbers of 1.61 and 2.01

An investigation has been made at Mach numbers of 1.61 and 2.01 and over a range of free-stream Reynolds number per foot from about 1.2 x 10(exp 6) to 8.3 x 10(exp 6) to determine the pressure distributions and wave drags due to two-dimensional fabrication-type surface roughness. Ten types of surface roughness, including step, wave, crease, and swept configurations were investigated. The tests were made on an ogive cylinder of fineness ratio 12.2, the roughness elements covering the cylindrical portion of the model. The results indicate that wave drag is the major component of the drag due to roughness at supersonic speeds. The pressure distributions over the roughness elements were generally found to be in good agreement with linearized two-dimensional theory except for regions of the elements affected by boundary-layer separation and shock detachment. There was little or no effect of Reynolds number except on the pressures within the regions influenced by separation or shock detachment. Inasmuch as most of the roughness configurations were affected by flow separation and shock detachment, there was generally an effect of Reynolds number on the roughness wave drag. This wave drag decreased as the free-stream Reynolds number was decreased.

Czarnecki, K. R.↗

Analysis and prediction of reaction kinetics using the degree of rate control

“Degree of rate control” (DRC) analysis provides a quantitative approach for analysing the kinetics of multi-step reaction mechanisms that has been widely applied to both heterogeneous and homogeneous catalysis research, as well as electrocatalysis. The DRC of any given transition state or intermediate is defined as a partial derivative such that it approximately equals the fractional increase in net rate to the product of interest per differential decrease in its standard-state free energy for that species (÷RT), holding constant the standard-state free energies of all other transition states and intermediates. Even very complex mechanisms usually have only a few species with non-zero DRCs and are thus the species whose interactions with the catalyst most strongly affect the net rate. These key DRC values thus offer a simple and intuitive route to optimize catalyst materials, especially with the assistance of computational methods like density functional theory (DFT). These high-DRC species are also the species whose energetics must be most accurately measured or calculated to achieve an accurate kinetic model for any reaction mechanism. In simple cases with a single “rate-determining step”, the DRC for its transition state (TS) is + 1. Catalyst-bound intermediates, on the other hand, often have negative DRCs equal to a small integer times their fractional occupancy of catalyst sites. The apparent activation energy equals a weighted average of the standard-state enthalpies (relative to reactants) of all the species (intermediates, transition states and products) in the reaction mechanism, each weighted by its DRC (+RT). It has been shown that the apparent transfer coefficient in electrocatalysis, an inverted form of the Tafel slope, is a weighted average of the number of electrons transferred to generate each intermediate or product species in the mechanism, weighted again by the DRC. Quantitative analysis of kinetic isotope effects (KIEs, or the ratio of net rates for different reactant isotopes) in complex mechanisms has shown that the logarithm of the KIE equals the weighted average over all species in the mechanism of the difference between the two isotopes in their standard-state free energies (÷RT), again weighted by the DRC. The reaction orders with respect to fluid-phase concentrations of reactants, products and intermediates have also been proven to be directly related to DRCs. Thus, there are numerous experimental observables which equate to short linear combinations of DRCs, so that combinations of experimental measurements might provide access to DRC values. Since its invention for transition states in 1994 and its generalization to include intermediates in 2009, the DRC has thus far mainly been calculated for microkinetic models based either entirely on DFT or on DFT with the key energies (i.e., those for high-DRC species) being fine-tuned to match experiments. The relationships summarized in this work provide new opportunities for using experiments earlier in the development and optimization of microkinetic models that require input from computational catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗