Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Complexity theory”

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 1,045 records · Page 58

Physics-Based Approach to Predict the Solar Activity Cycles

Observations of the complex highly non-linear dynamics of global turbulent flows and magnetic fields are currently available only from Earth-side observations. Recent progress in helioseismology has provided us some additional information about the subsurface dynamics, but its relation to the magnetic field evolution is not yet understood. These limitations cause uncertainties that are difficult take into account, and perform proper calibration of dynamo models. The current dynamo models have also uncertainties due to the complicated turbulent physics of magnetic field generation, transport and dissipation. Because of the uncertainties in both observations and theory, the data assimilation approach is natural way for the solar cycle prediction and estimating uncertainties of this prediction. I will discuss the prediction results for the upcoming Solar Cycle 25 and their uncertainties and affect of Ensemble Kalman Filter parameters to resulting predictions.

Physics-Based↗

Production and Potential Detection of Functionalized Hexamethylene-Tetramine Compounds in Space

Laboratory studies have shown that exposure of mixed ices of astrophysical interest to ionizing radiation such as ultraviolet (UV) photons or energetic particles (electrons, protons) leads to the production of large numbers of new, more complex compounds. A significant portion of these new species appear to belong to a family of molecules that consist of hexamethylenetetramine (HMT; C6N4H12) and HMT to which different chemical side groups have been substituted for a peripheral H atom. This work presents the identification of HMT-methanol (HMT-CH2OH), one of these HMT variants, in organic residues produced from the UV irradiation of astrophysically relevant ice mixtures at < 20 K. We also present the infrared (IR) spectra of HMT, HMTCH2OH, and a number of other HMT variants computed using density functional theory (DFT) computations. These spectra can be compared with each other and show similarities that can be used to search for this family of compounds in space.

UV Irradiation↗

S-Wave Dispersion Relations: Exact Left Hand E-Plane Discontinuity from the Born Series

We show, for a superposition of Yukawa potentials, that the left hand cut discontinuity in the complex E plane of the (S-wave) scattering amplitude is given exactly, in an interval depending on n, by the discontinuity of the Born series stopped at order n. This also establishes an inverse and unexpected correspondence of the Born series at positive high energies and negative low energies. We can thus construct a viable dispersion relation (DR) for the partial (S-) wave amplitude. The high numerical precision achievable by the DR is demonstrated for the exponential potential at zero scattering energy. We also briefly discuss the extension of our results to Field Theory.

Bessis, D.↗

Subgrid Scale Modeling in Solar Convection Simulations using the ASH Code

The turbulent solar convection zone has remained one of the most challenging and important subjects in physics. Understanding the complex dynamics in the solar con- vection zone is crucial for gaining insight into the solar dynamo problem. Many solar observatories have generated revealing data with great details of large scale motions in the solar convection zone. For example, a strong di erential rotation is observed: the angular rotation is observed to be faster at the equator than near the poles not only near the solar surface, but also deep in the convection zone. On the other hand, due to the wide range of dynamical scales of turbulence in the solar convection zone, both theory and simulation have limited success. Thus, cutting edge solar models and numerical simulations of the solar convection zone have focused more narrowly on a few key features of the solar convection zone, such as the time-averaged di erential rotation. For example, Brun & Toomre (2002) report computational finding of differential rotation in an anelastic model for solar convection. A critical shortcoming in this model is that the viscous dissipation is based on application of mixing length theory to stellar dynamics with some ad hoc parameter tuning. The goal of our work is to implement the subgrid scale model developed at CTR into the solar simulation code and examine how the differential rotation will be a affected as a result. Specifically, we implement a Smagorinsky-Lilly subgrid scale model into the ASH (anelastic spherical harmonic) code developed over the years by various authors. This paper is organized as follows. In x2 we briefly formulate the anelastic system that describes the solar convection. In x3 we formulate the Smagorinsky-Lilly subgrid scale model for unstably stratifed convection. We then present some preliminary results in x4, where we also provide some conclusions and future directions.

Young, Y.-N.↗

Composite laminate optimization program suitable for microcomputers

The design of a composite panel requires some way of finding the minimum thickness laminate which will withstand the load requirements without failure. The mathematical complexity of this problem dictates the use of nonlinear optimization techniques. Specialized laminate optimization programs were developed which are compact and efficient enough to run on microcomputers. Only stresses at a point and inplane loads and deflections are considered. The programs are simple to use and require no knowledge of optimization. Techniques are developed which find minimum thickness laminates with either ply ratios or ply angles as design variables. A method is presented for finding the optimum orientation for the axis of symmetry of an orthotropic laminate. The orthotropic laminate program uses an approximate failure theory, which speed up computations dramatically.

Flanagan, G. V.↗

Uncertainty modelling for the control of an active structure

Two approaches to uncertainty modeling of a precision active truss structure for control design are presented. The first approach is based on the H(infinity) optimal control theory, which is used for robust control design with an additive uncertainty model. It is designed to provide vibration suppression to selected modes of the structure at the accelerometer locations. The uncertainty model can be used as an adjustable design parameter which will enable a study of the trade-offs between performance and robust stability. The second approach is based on mu-synthesis. Although the structured singular value can be used to analyze the robustness of any type of uncertainty, there exist difficulties in synthesis if the uncertainty block contains both real and complex variations. Mu-synthesis techniques do not allow for real-parameter variations in the synthesis procedure. The use of a complex variation description enclosing the real variation is suggested.

Blackwood, G. H.↗

A unified process for systems identification based on performance assessment

Many modern spacecraft are complex multi-body dynamic systems where the bodies are connected by several active control systems for pointing and isolation. Mission requirements indicate that many structural modes and possibly some nonlinear effects will require characterization. Thus, even characterization at the subsystem level will become more difficult than usual. System level characterization difficulties will be compounded by the fact that only limited ground testing will be possible on the full up system, and flight testing will be restricted by an extremely limited measurements set. The object of the present discussion is the application of matrix majorant theory to the problem of assessing dynamic system performance when knowledge of the system is uncertain. We show how majorants provide an effective tool to relate required performance output to system identification test quality in terms of residual uncertainty in input-output relations, parameter values, nonlinearities, and interactions. The underlying machinery consists of the block-norm matrix which is a nonnegative matrix each of whose elements is the norm of a block of a suitably partitioned matrix. A matrix which bounds the block-norm matrix in the sense of nonnegative matrices, i.e., element by element is known as a majorant.

Hyland, C.↗

Geomagnetic Cutoff Rigidity Computer Program: Theory, Software Description and Example

The access of charged particles to the earth from space through the geomagnetic field has been of interest since the discovery of the cosmic radiation. The early cosmic ray measurements found that cosmic ray intensity was ordered by the magnetic latitude and the concept of cutoff rigidity was developed. The pioneering work of Stoermer resulted in the theory of particle motion in the geomagnetic field, but the fundamental mathematical equations developed have 'no solution in closed form'. This difficulty has forced researchers to use the 'brute force' technique of numerical integration of individual trajectories to ascertain the behavior of trajectory families or groups. This requires that many of the trajectories must be traced in order to determine what energy (or rigidity) a charged particle must have to penetrate the magnetic field and arrive at a specified position. It turned out the cutoff rigidity was not a simple quantity but had many unanticipated complexities that required many hundreds if not thousands of individual trajectory calculations to solve. The accurate calculation of particle trajectories in the earth's magnetic field is a fundamental problem that limited the efficient utilization of cosmic ray measurements during the early years of cosmic ray research. As the power of computers has improved over the decades, the numerical integration procedure has grown more tractable, and magnetic field models of increasing accuracy and complexity have been utilized. This report is documentation of a general FORTRAN computer program to trace the trajectory of a charged particle of a specified rigidity from a specified position and direction through a model of the geomagnetic field.

Smart, D. F.↗

Error-Rate Bounds for Coded PPM on a Poisson Channel

Equations for computing tight bounds on error rates for coded pulse-position modulation (PPM) on a Poisson channel at high signal-to-noise ratio have been derived. These equations and elements of the underlying theory are expected to be especially useful in designing codes for PPM optical communication systems. The equations and the underlying theory apply, more specifically, to a case in which a) At the transmitter, a linear outer code is concatenated with an inner code that includes an accumulator and a bit-to-PPM-symbol mapping (see figure) [this concatenation is known in the art as "accumulate-PPM" (abbreviated "APPM")]; b) The transmitted signal propagates on a memoryless binary-input Poisson channel; and c) At the receiver, near-maximum-likelihood (ML) decoding is effected through an iterative process. Such a coding/modulation/decoding scheme is a variation on the concept of turbo codes, which have complex structures, such that an exact analytical expression for the performance of a particular code is intractable. However, techniques for accurately estimating the performances of turbo codes have been developed. The performance of a typical turbo code includes (1) a "waterfall" region consisting of a steep decrease of error rate with increasing signal-to-noise ratio (SNR) at low to moderate SNR, and (2) an "error floor" region with a less steep decrease of error rate with increasing SNR at moderate to high SNR. The techniques used heretofore for estimating performance in the waterfall region have differed from those used for estimating performance in the error-floor region. For coded PPM, prior to the present derivations, equations for accurate prediction of the performance of coded PPM at high SNR did not exist, so that it was necessary to resort to time-consuming simulations in order to make such predictions. The present derivation makes it unnecessary to perform such time-consuming simulations.

Moision, Bruce↗

Strategic Technology Investment Analysis: An Integrated System Approach

Complex technology investment decisions within NASA are increasingly difficult to make such that the end results are satisfying the technical objectives and all the organizational constraints. Due to a restricted science budget environment and numerous required technology developments, the investment decisions need to take into account not only the functional impact on the program goals, but also development uncertainties and cost variations along with maintaining a healthy workforce. This paper describes an approach for optimizing and qualifying technology investment portfolios from the perspective of an integrated system model. The methodology encompasses multi-attribute decision theory elements and sensitivity analysis. The evaluation of the degree of robustness of the recommended portfolio provides the decision-maker with an array of viable selection alternatives, which take into account input uncertainties and possibly satisfy nontechnical constraints. The methodology is presented in the context of assessing capability development portfolios for NASA technology programs.

decision support↗

Compass Concurrent Engineering Lessons Learned in Remote and Hybrid Environments

The Compass Team at NASA’s Glenn Research Center (GRC) is a concurrent engineering team which specializes in conceptual spacecraft mission designs. Detailed descriptions of the team, its history, and its operating model can be found in [1] and [2]. During the COVID-19 pandemic, the team was required to move to remote (virtual) operation from their in-person model for approximately 22 months. As the restrictions began to lift and team members were able to return in-person to the Compass Lab, the team moved into a hybrid mode of operation, with some participants still tying in remotely some or all of the time. This paper discusses many of the lessons learned from these experiences, highlighting improvements, outstanding challenges and the tools and methodologies used to address both. Throughout this discussion the terms “in-person”, “hybrid”, “virtual” and “remote” will be used. For the purposes of this paper, “in-person” will be understood to mean when team members are interacting simultaneously, physically within the Compass Lab. “Remote” or “virtual” will refer to when interactions are happening between people who are not co-located using only technology to interface. “Hybrid” will refer to when two or more participants are physically located in the Compass Lab and one or more participant(s) is participating remotely. Media richness is described as “a medium’s ability to communicate effectively based of four factors. They are the capacity for immediate feedback, the number of cues and channels it utilizes, the degree of personalization it affords, and its ability to communicate using natural language” [3]. This theory will be referenced and discussed in multiple of the following sections due to its relevance when choosing how to operate in remote and hybrid modes, as well as in making tool selections. The key to selecting the best mode of communication lays in how complicated the discussion is and the level of ambiguity involved. Not all conversations or interactions require media rich mediums. For example, providing information about which there is little to no ambiguity can easily be done in less rich methods - such as email. A conversation including high levels of ambiguity and/or complex information is better suited to a richer medium, such as in-person or a video call with shared screens.

Concurrent Engineering↗

Theory and Experiment of Multielement Airfoils: A Comparison

A detailed comparison of computed and measured pressure distributions, velocity profiles, transition onset, and Reynolds shear stresses for multi-element airfoils is presented. It is shown that the transitional k-zeta model, which is implemented into CFL3D, does a good job of predicting pressure distributions, transition onset, and velocity profiles with the exception of velocities in the slat wake region. Considering the fact that the hot wire used was not fine enough to resolve Reynolds stresses in the boundary layer, comparisons of turbulence stresses varied from good to fair. It is suggested that the effects of unsteadiness be thoroughly evaluated before more complicated transition/turbulence models are used. Further, it is concluded that the present work presents a viable and economical method for calculating laminar/transitional/turbuient flows over complex shapes without user interface.

Czerwiec, Ryan↗

Theory and Simulations of Solar System Plasmas

"Theory and simulations of solar system plasmas" aims to highlight results from microscopic to global scales, achieved by theoretical investigations and numerical simulations of the plasma dynamics in the solar system. The theoretical approach must allow evidencing the universality of the phenomena being considered, whatever the region is where their role is studied; at the Sun, in the solar corona, in the interplanetary space or in planetary magnetospheres. All possible theoretical issues concerning plasma dynamics are welcome, especially those using numerical models and simulations, since these tools are mandatory whenever analytical treatments fail, in particular when complex nonlinear phenomena are at work. Comparative studies for ongoing missions like Cassini, Cluster, Demeter, Stereo, Wind, SDO, Hinode, as well as those preparing future missions and proposals, like, e.g., MMS and Solar Orbiter, are especially encouraged.

Goldstein, Melvyn L.↗

Rapid Aero Modeling for Urban Air Mobility Aircraft in Computational Experiments

Rapid Aero Modeling (RAM) applied to computational testing, RAM-C, is an approach to efficiently and automatically obtain aerodynamic models during computational investigations. RAM-C is designed to estimate models appropriate for flight dynamics studies and simulations. The approach responds to a demand for experimental efficiency and model fidelity that has increased with growing aircraft complexity and aerodynamic nonlinearities associated with hybrid and electric vertical takeoff and landing (eVTOL) aircraft. In an Urban Air Mobility (UAM) transportation system, it is expected that aircraft will embrace many features from both airplanes and rotorcraft. These vehicles present many more factors than conventional aircraft which can lead to increased computational costs and missed key factor interactions when applying traditional testing and modeling methods. RAM-C provides feedback loops around computational codes to rapidly guide testing toward aerodynamic models meeting user-defined fidelity goals. It combines and extends concepts from design of experiment theory and aircraft system identification theory that allow the user the freedom to choose, in advance of the test, a specific level of fidelity in terms of prediction error. RAM-C only collects enough data required to meet the user-specified prediction error requirements thus saving computational time and resources. The overall achievable fidelity of the final model also depends on the accuracy of the test facility, or in this case, the computational modeling approach. Previous studies to support development of the RAM-T process were conducted in wind tunnel tests to assess potential metrics, algorithms, and procedures. This paper presents results from the next steps taken and tests conducted for the development of RAM-C technology and highlights some of the unique features of RAM applied eVTOL configurations in a computational study.

Aerodynamics↗

Characteristics of finite amplitude stationary gravity waves in the atmosphere of Venus

This paper extends the study of stationary gravity waves generated near the surface of Venus reported previously by Young et al. to include finite amplitude effects associated with large amplitude waves. Waves are forced near the surface of Venus by periodic forcing. The height-dependent profiles of static stability and mean wind in the Venus atmosphere play a very important role in the evolution of the nonlinear behavior of the waves, just as they do in the linear wave solutions. Certain wave properties are qualitatively consistent with linear wave theory, such as wave trapping, resonance, and wave evanescence for short horizontal wavelenghts. However, the finite amplitude solutions also exhibit many other interesting features. In particular, for forcing amplitudes representative of those that could be expected in mountainous regions such as Aphrodite Terra, waves generated near the surface can reach large amplitudes at and above cloud levels, with clear signatures in the circulation pattern. At still higher levels, the waves can reach large enough amplitude to break, unless damping rates above the clouds are sufficient to limit wave amplitude growth. Well below cloud levels the waves develop complex flow patterns as the result of finite amplitude wave-wave interactions, and waves are generated having considerably shorter horizontal wavelenghts than that associated with the forcing near the surface. Nonlinear interactions can excite waves that are resonant with the background wind and static stability fields even when the primary surface forcing does not, and these waves can dominate the wave spectrum near cloud levels. A global map of Venus topographic slopes derived from Magellan altimetry data shows that slopes of magnitude comparable to or exceeding that used to force the model are ubiquitous over the surface.

Young, Richard E.↗

Origin of the Galactic Disk 6.7 kev Line Emission

The goal of this program was to determine if the extended FeXXV 6.7 kev line emission might possibly be produced and confined by the hot wind-shocked bubbles to accompany UC HII regions. The main result of this study are: (1) FeXXV is detected in the W3 complex, but at a level that could only explain a small fraction of the galactic disk emission if all UC HII regions emit at about the same intensity as the W3 complex; (2) Two X-ray sources are detected in W3. W3-X 1 coincides with the radio image of this region, but W3-X2 has no radio, optical, or infrared counterpart; (3) There is no evidence for variability of W3-X1 during the period of observations (approx, 40,000 sec); (4) The X-ray spectrum of W3-X1 has no emission shortward of 1 kev, it peaks at approx. 2 kev and show significant emission out to approx. 6 kev. No individual lines are resolved. There is currently no generally accepted theory for extended hard X-ray emission in HII regions. Perhaps the most significant discovery of this program has been the detection of extended hard X-rays and the realization that some entirely new processes must be invoked to understand this; and (5)A minimum (chi)(sup 2) fit of the spectrum implies a H absorbing column of N(sub H) approx, equals to 2.1 x 10(exp 22)/ cm, a temperature of the emitting plasma of 7 x 10(exp 7) K, and a luminosity of approx. equal to 10(33)erg/s.

Churchwell, Ed↗

The Relaxation Matrix for Symmetric Tops with Inversion Symmetry. I. Effects of Line Coupling on Self-Broadened v (sub 1) and Pure Rotational Bands of NH3

The Robert-Bonamy formalism has been commonly used to calculate half-widths and shifts of spectral lines for decades. This formalism is based on several approximations. Among them, two have not been fully addressed: the isolated line approximation and the neglect of coupling between the translational and internal motions. Recently, we have shown that the isolated line approximation is not necessary in developing semi-classical line shape theories. Based on this progress, we have been able to develop a new formalism that enables not only to reduce uncertainties on calculated half-widths and shifts, but also to model line mixing effects on spectra starting from the knowledge of the intermolecular potential. In our previous studies, the new formalism had been applied to linear and asymmetric-top molecules. In the present study, the method has been extended to symmetric-top molecules with inversion symmetry. As expected, the inversion splitting induces a complete failure of the isolated line approximation. We have calculated the complex relaxation matrices of selfbroadened NH3. The half-widths and shifts in the ν1 and the pure rotational bands are reported in the present paper. When compared with measurements, the calculated half-widths match the experimental data very well, since the inapplicable isolated line approximation has been removed. With respect to the shifts, only qualitative results are obtained and discussed. Calculated off-diagonal elements of the relaxation matrix and a comparison with the observed line mixing effects are reported in the companion paper (Paper II).

Trajectory models↗

Composite propellant combustion modeling with a porous plate burner

A burner is designed to model on a large scale (millimeters) the complex vapor phase processes that take place on awkwardly small scales (about a hundred microns) in the combustion of practical AP/ composite propellants at conventional pressures. Binder vapor evolution is modeled with gaseous fuel (ethane in the experiments reported) flow through a porous plate and oxidizer vapor with the flow of a gaseous oxidizer (air and enriched air in the experiments reported) through discrete holes in the porous plate. Measured flame standoff distance and surface temperature variations are consistent with theoretical predictions at atmospheric pressure. Data obtained at several atmospheres are discussed in the light of the current theories of propellant burning that place varying emphasis on the roles of chemical kinetics and fluid dynamic diffusion/mixing in the vapor phase of a burning composite propellant. The potential and proposed future applications of the burner are indicated.

Kumar, R. N.↗