Search NASA⌕ Search

SEARCH · Search NASA

Results for “small parameters”

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 163 records · Page 9

Results of including geometric nonlinearities in an aeroelastic model of an F/A-18

An integrated, nonlinear simulation model suitable for aeroelastic modeling of fixed-wing aircraft has been developed. While the author realizes that the subject of modeling rotating, elastic structures is not closed, it is believed that the equations of motion developed and applied herein are correct to second order and are suitable for use with typical aircraft structures. The equations are not suitable for large elastic deformation. In addition, the modeling framework generalizes both the methods and terminology of non-linear rigid-body airplane simulation and traditional linear aeroelastic modeling. Concerning the importance of angular/elastic inertial coupling in the dynamic analysis of fixed-wing aircraft, the following may be said. The rigorous inclusion of said coupling is not without peril and must be approached with care. In keeping with the same engineering judgment that guided the development of the traditional aeroelastic equations, the effect of non-linear inertial effects for most airplane applications is expected to be small. A parameter does not tell the whole story, however, and modes flagged by the parameter as significant also need to be checked to see if the coupling is not a one-way path, i.e., the inertially affected modes can influence other modes.

Buttrill, Carey S.↗

A Newton/upwind method and numerical study of shock wave/boundary layer interactions

The objective of the paper is two-fold. First, an upwind/central differencing method for solving the steady Navier-Stokes equations is described. The symmetric line relation method is used to solve the resulting algebraic system to achieve high computational efficiency. The grid spacings used in the calculations are determined from the triple-deck theory, in terms of Mach and Reynolds numbers and other flow parameters. Thus the accuracy of the numerical solutions is improved by comparing them with experimental, analytical, and other computational results. Secondly, the shock wave/boundary layer interactions are studied numerically, with special attention given to the flow separation. The concept of free interaction is confirmed. Although the separated region varies with Mach and Reynolds numbers, it is found that the transverse velocity component behind the incident shock, which has not been identified heretofore, is also an important parameter. A small change of this quantity is sufficient to eliminate the flow separation entirely.

Liou, Meng-Sing↗

Self-calibration of robot-sensor system

The process of finding the coordinate transformation between a robot and an external sensor system has been addressed. This calibration is equivalent to solving a nonlinear optimization problem for the parameters that characterize the transformation. A two-step procedure is herein proposed for solving the problem. The first step involves finding a nominal solution that is a good approximation of the final solution. A varational problem is then generated to replace the original problem in the next step. With the assumption that the variational parameters are small compared to unity, the problem that can be more readily solved with relatively small computation effort.

Yeh, Pen-Shu↗

A variational approach to dynamics of flexible multibody systems

This paper presents a variational formulation of constrained dynamics of flexible multibody systems, using a vector-variational calculus approach. Body reference frames are used to define global position and orientation of individual bodies in the system, located and oriented by position of its origin and Euler parameters, respectively. Small strain linear elastic deformation of individual components, relative to their body references frames, is defined by linear combinations of deformation modes that are induced by constraint reaction forces and normal modes of vibration. A library of kinematic couplings between flexible and/or rigid bodies is defined and analyzed. Variational equations of motion for multibody systems are obtained and reduced to mixed differential-algebraic equations of motion. A space structure that must deform during deployment is analyzed, to illustrate use of the methods developed.

Wu, Shih-Chin↗

The effect of suspended particles on Jean's criterion for gravitational instability

The effect that the proper inclusion of suspended particles has on Jeans' criterion for the self-gravitational instability of an unbounded nonrotating adiabatic gas cloud is examined by formulating the appropriate model system, introducing particular physically plausible equations of state and constitutive relations, performing a linear stability analysis of a uniformly expanding exact solution to these governing equations, and exploiting the fact that there exists a natural small material parameter for this problem given by N sub 1/n sub 1, the ratio of the initial number density for the particles to that for the gas. The main result of this investigation is the derivation of an altered criterion which can substantially reduce Jeans' original critical wavelength for instability. It is then shown that the existing discrepancy between Jeans' theoretical prediction using and actual observational data relevant to the Andromeda nebula M31 can be accounted for by this new criterion of assuming suspended particles of a reasonable grain size and distribution to be present.

Wollkind, David J.↗

Techniques for developing approximate optimal advanced launch system guidance

An extension to the authors' previous technique used to develop a real-time guidance scheme for the Advanced Launch System is presented. The approach is to construct an optimal guidance law based upon an asymptotic expansion associated with small physical parameters, epsilon. The trajectory of a rocket modeled as a point mass is considered with the flight restricted to an equatorial plane while reaching an orbital altitude at orbital injection speeds. The dynamics of this problem can be separated into primary effects due to thrust and gravitational forces, and perturbation effects which include the aerodynamic forces and the remaining inertial forces. An analytic solution to the reduced-order problem represented by the primary dynamics is possible. The Hamilton-Jacobi-Bellman or dynamic programming equation is expanded in an asymptotic series where the zeroth-order term (epsilon = 0) can be obtained in closed form.

Feeley, Timothy S.↗

A transfer matrix approach to vibration localization in mistuned blade assemblies

A study of mode localization in mistuned bladed disks is performed using transfer matrices. The transfer matrix approach yields the free response of a general, mono-coupled, perfectly cyclic assembly in closed form. A mistuned structure is represented by random transfer matrices, and the expansion of these matrices in terms of the small mistuning parameter leads to the definition of a measure of sensitivity to mistuning. An approximation of the localization factor, the spatially averaged rate of exponential attenuation per blade-disk sector, is obtained through perturbation techniques in the limits of high and low sensitivity. The methodology is applied to a common model of a bladed disk and the results verified by Monte Carlo simulations. The easily calculated sensitivity measure may prove to be a valuable design tool due to its system-independent quantification of mistuning effects such as mode localization.

Ottarson, Gisli↗

Density and mixture fraction measurements in a GO2/GH2 uni-element rocket chamber

In recent years, there has been a renewed interest in gas/gas injectors for rocket combustion. Specifically, the proposed new concept of full-flow oxygen rich preburner systems calls for the injection of both oxygen and hydrogen into the main chamber as gaseous propellants. The technology base for gas/gas injection must mature before actual booster class systems can be designed and fabricated. Since the data base for gas/gas injection is limited to studies focusing on the global parameters of small reaction engines, there is a critical need for experiment programs that emphasize studying the mixing and combustion characteristics of GO2 and GH2 propellants from a uni-element injector point of view. The experimental study of the combusting GO2/GH2 propellant combination in a uni-element rocket chamber also provides a simplified environment, in terms of both geometry and chemistry, that can be used to verify and validate computational fluid dynamic (CFD) models.

Moser, M. D.↗

Parametric instabilities of parallel-propagating Alfven waves: Some analytical results

We consider the stability of a circularly polarized Alfven wave (the pump wave) which propagates parallel to the ambient magnetic field. Only parallel-propagating perturbations are considered, and we ignore dispersive effects due to the ion cyclotron frequency. The dissipationless MHD equations are used throughout; thus possibibly important effects arising from Landau and transit time damping are omitted. We derive a series of analytical approximations to the dispersion relation using A = (Delta B/B(sub O))(exp 2) as a small expansion parameter; Delta B is the pump amplitude, and B(sub O) is the ambient magnetic field strength. We find that the plasma beta (the square of the ratio of the sound speed to the Alfven speed) plays a crucial role in determining the behavior of the parametric instabilities of the pump. If 0 less than beta less than 1 we find the familiar result that the pump decays into a forward propagating sound wave and a backward propagating Alfven wave with maximum growth rate gamma(sub max) varies A(sup 1/2), but beta cannot be too close to 0 or to 1. If beta approx. 1, we find gamma(sub max) varies A(sup 3/4), if beta greater than 1, we find gamma(sub max) varies A(sup 3/2), while if beta approx. 0, we obtain gamma(sub max) varies A(sup 1/3); moreover, if beta approx. 0 there is a nearly purely growing instability. In constrast to the familiar decay instability, for which the backward propagating Alfven wave has lower frequency and wavenumber than the pump, we find that if beta greater than or approx. equal to 1 the instability is really a beat instability which is dominated by a transverse wave which is forward propagating and has frequency and wavenumber which are nearly twice the pump values. Only the decay instability for 0 less than beta less than 1 can be regarded as producing two recognizable normal modes, namely, a sound wave and an Alfven wave. We discuss how the different characteristics of the instabilities may affect the evolution of Alfven waves in the solar wind. However, for a solar wind in which beta approx. 1 the growth times of the instabilities are probably too long for these instabilities to have an appreciable effect inside 1 AU.

Jayanti, V.↗

A Numerical Comparison of Barrier and Modified Barrier Methods for Large-Scale Bound-Constrained Optimization

When a classical barrier method is applied to the solution of a nonlinear programming problem with inequality constraints, the Hessian matrix of the barrier function becomes increasingly ill-conditioned as the solution is approached. As a result, it may be desirable to consider alternative numerical algorithms. We compare the performance of two methods motivated by barrier functions. The first is a stabilized form of the classical barrier method, where a numerically stable approximation to the Newton direction is used when the barrier parameter is small. The second is a modified barrier method where a barrier function is applied to a shifted form of the problem, and the resulting barrier terms are scaled by estimates of the optimal Lagrange multipliers. The condition number of the Hessian matrix of the resulting modified barrier function remains bounded as the solution to the constrained optimization problem is approached. Both of these techniques can be used in the context of a truncated-Newton method, and hence can be applied to large problems, as well as on parallel computers. In this paper, both techniques are applied to problems with bound constraints and we compare their practical behavior.

Nash, Stephen G.↗

STS-93 SSME Nozzle Tube Rupture Investigation

STS-93 was launched on July 23, 1999. There was an anomaly at the end of the launch in that the main engines shut down 0.16 second early because sensors detected a low level of oxidizer in the LOX tank (actually the duct from the tank to the vehicle). This resulted in a cutoff velocity for the vehicle that was 16 ft/sec low. It should have been 25872 ft/sec. The OMS engines were subsequently used to achieve the proper orbit. An investigation was immediately initiated into the cause of this LOX tank low level cutoff. It was noticed during the launch that the turbine temperatures for Main Engine 3 (E2019) were approximately 100 F higher than the preflight prediction. Linear Engine Model matching of the data indicated that a nozzle leak best fit the data. Post launch review of the data showed, that at approximately five seconds into the start, numerous parameters indicated small anomalous shifts. These shifts were all consistent with a rupture of nozzle tubes. Post launch review of the films showed that just after Space Shuttle Main Engine (SSME) ignition and just prior to liftoff a streak is seen in the exhaust plume of E2019. Just after liftoff the streak can be seen emanating from the nozzle wall. This photo confirmed that a leak was coming from the nozzle tubes. Based on the photo, the axial location of the leak was estimated to be 28 in. from the aft end of the nozzle and in line with nozzle coolant feed line #1. The streak continued to be visible during the launch.

Romine, W. Dennis↗

Inherent and Apparent Scattering Properties of Coated or Uncoated Spheres Embedded in an Absorbing Host Medium

The conventional Lorenz-Mie formalism is extended to the scattering process associated with a coated sphere embedded in an absorbing medium. It is shown that apparent and inherent scattering cross sections of a scattering particle, which are identical in the case of transparent host medium, are different if the host medium is absorptive. Here the inherent single-scattering properties are derived from the near-field information whereas the corresponding apparent counterparts are derived from the far-field asymptotic form of the scattered wave with scaling of host absorption that is assumed to be in an exponential form. The formality extinction and scattering efficiencies defined in the same manner as in the conventional sense can be unbounded. For a nonabsorptive particle embedded in an absorbing medium, the effect of host absorption on the phase matrix elements associated with polarization is significant. This effect, however, is largely reduced for strongly absorptive particles such as soot. For soot particles coated with water, the impurity can substantially reduce the single-scattering albedo of the particle if the size parameter is small. For water-coating soot and hollow ice spheres, it is shown that the phase matrix elements -P(sub 12)/P(sub 11) and P(sub 33)/P(sub 11) are unique if the shell is thin, as compared with the case for thick shell. Furthermore, the radiative transfer equation regarding a multidisperse particle system in an absorbing medium is discussed. It is illustrated that the conventional computation algorithms can be applied to solve the multiple scattering process if the scaled apparent single-scattering properties are applied.

Yang, P.↗

Electromagnetic Dissociation Cross Sections for High LET Fragments

Nuclear interaction cross sections are used in space radiation transport codes to calculate the probability of fragment emission in high energy nucleus-nucleus collisions. Strong interactions usually dominate in these collisions, but electromagnetic (EM) interactions can also sometimes be important. Strong interactions typically occur when the projectile nucleus hits a target nucleus, with a small impact parameter. For impact parameters larger than the sum of the nuclear radii, EM reactions dominate and the process is called electromagnetic dissociation (EMD) if one of the nuclei undergo fragmentation. Previous models of EMD have been used to calculate single proton (p) production, single neutron (n) production or light ion production, where a light ion is defined as an isotope of hydrogen (H) or helium (He), such as a deuteron (2H), a triton (3H), a helion (3He) or an alpha particle (4He). A new model is described which can also account for multiple nucleon production, such as 2p, 2n, 1p1n, 2p1n, 2p2n, etc. in addition to light ion production. Such processes are important to include for the following reasons. Consider, for example, the EMD reaction 56Fe + Al --> 52Cr + X + Al, for a 56Fe projectile impacting Al, which produces the high linear energy transfer (LET) fragment 52Cr. In this reaction, the most probable particles representing X are either 2p2n or 4He. Therefore, production of the high LET fragment 52Cr, must include the multiple nucleon production of 2p2n in addition to the light ion production of 4He. Previous models, such as the NUCFRG3 model, could only account for the 4He production process in this reaction and could not account for 2p2n. The new EMD model presented in this work accounts for both the light ion and multiple nucleon processes, and is therefore able to correctly account for the production of high LET products such as 52Cr. The model will be described and calculations will be presented that show the importance of light ion and multiple nucleon production. The work will also show that EMD reactions contribute most to those fragments with the highest LET.

Norbury, John↗

An information-matching approach to optimal experimental design and active learning

The efficacy of mathematical models heavily depends on the quality of the training data, yet collecting sufficient data is often expensive and challenging. Many modeling applications require inferring parameters only as a means to predict other quantities of interest (QoI). Because models often contain many unidentifiable (sloppy) parameters, QoIs often depend on a relatively small number of parameter combinations. Therefore, we introduce an information-matching criterion based on the Fisher information matrix to select the most informative training data from a candidate pool. This method ensures that the selected data contain sufficient information to learn only those parameters that are needed to constrain downstream QoIs. It is formulated as a convex optimization problem, making it scalable to large models and datasets. Here, we demonstrate the effectiveness of this approach across various modeling problems in diverse scientific fields, including power systems and underwater acoustics. Finally, we use information-matching as a query function within an active learning (AL) loop for materials science applications. In all these applications, we find that a relatively small set of optimal training data can provide the necessary information for achieving precise predictions. These results are encouraging for diverse future applications, particularly AL in large machine-learning models.

Materials science↗

Environment parameters and basic functions for floating-point computation

A language-independent proposal for environment parameters and basic functions for floating-point computation is presented. Basic functions are proposed to analyze, synthesize, and scale floating-point numbers. The model provides a small set of parameters and a small set of axioms along with sharp measures of roundoff error. The parameters and functions can be used to write portable and robust codes that deal intimately with the floating-point representation. Subject to underflow and overflow constraints, a number can be scaled by a power of the floating-point radix inexpensively and without loss of precision. A specific representation for FORTRAN is included.

Brown, W. S.↗

Resonant frequency calculations using a hybrid perturbation-Galerkin technique

A two-step hybrid perturbation Galerkin technique is applied to the problem of determining the resonant frequencies of one or several degree of freedom nonlinear systems involving a parameter. In one step, the Lindstedt-Poincare method is used to determine perturbation solutions which are formally valid about one or more special values of the parameter (e.g., for large or small values of the parameter). In step two, a subset of the perturbation coordinate functions determined in step one is used in Galerkin type approximation. The technique is illustrated for several one degree of freedom systems, including the Duffing and van der Pol oscillators, as well as for the compound pendulum. For all of the examples considered, it is shown that the frequencies obtained by the hybrid technique using only a few terms from the perturbation solutions are significantly more accurate than the perturbation results on which they are based, and they compare very well with frequencies obtained by purely numerical methods.

Geer, James F.↗

Resonant frequency calculations using a hybrid perturbation-Galerkin technique

A two-step hybrid perturbation Galerkin technique is applied to the problem of determining the resonant frequencies of one or several degrees of freedom nonlinear systems involving a parameter. In one step, the Lindstedt-Poincare method is used to determine perturbation solutions which are formally valid about one or more special values of the parameter (e.g., for large or small values of the parameter). In step two, a subset of the perturbation coordinate functions determined in step one is used in Galerkin type approximation. The technique is illustrated for several one degree of freedom systems, including the Duffing and van der Pol oscillators, as well as for the compound pendulum. For all of the examples considered, it is shown that the frequencies obtained by the hybrid technique using only a few terms from the perturbation solutions are significantly more accurate than the perturbation results on which they are based, and they compare very well with frequencies obtained by purely numerical methods.

Geer, James F.↗

One‐at‐a‐Time Parameter Perturbation Ensemble of the Community Land Model, Version 5.1

Comprehensive land models are subject to significant parametric uncertainty, which can be hard to quantify due to the large number of parameters and high model computational costs. We constructed a large parameter perturbation ensemble (PPE) for the Community Land Model version 5.1 with biogeochemistry configuration (CLM5.1-BGC). We performed more than 2,000 simulations perturbing 211 parameters across six forcing scenarios. This provides an expansive data set, which can be used to identify the most influential parameters on a wide range of output variables globally, by biome, or by plant functional type. We found that parameter effects can exceed scenario effects and that a small number of parameters explains a large fraction of variance across our ensemble. The most important parameters can differ regionally and also based on the forcing scenario. The software infrastructure developed for this experiment has greatly reduced the human and computer time needed for CLM PPEs, which can facilitate routine investigation of parameter sensitivity and uncertainty, as well as automated calibration.

Kennedy, Daniel [NSF National Center for Atmospher↗