Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “domain structure”

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,027 records · Page 57

A CFD study of complex missile and store configurations in relative motion

An investigation was conducted from May 16, 1990 to August 31, 1994 on the development of computational fluid dynamics (CFD) methodologies for complex missiles and the store separation problem. These flowfields involved multiple-component configurations, where at least one of the objects was engaged in relative motion. The two most important issues that had to be addressed were: (1) the unsteadiness of the flowfields (time-accurate and efficient CFD algorithms for the unsteady equations), and (2) the generation of grid systems which would permit multiple and moving bodies in the computational domain (dynamic domain decomposition). The study produced two competing and promising methodologies, and their proof-of-concept cases, which have been reported in the open literature: (1) Unsteady solutions on dynamic, overlapped grids, which may also be perceived as moving, locally-structured grids, and (2) Unsteady solutions on dynamic, unstructured grids.

Baysal, Oktay↗

TE/TM Simulations of Interferometric Measurements

Interferometric synthetic aperture radar (IFSAR) measurements at X-, C-, L-, and P-band are used to derive ground topography at meter level resolution. Interpretation of the derived topography requires attention due to the complex interaction of the radar signal with ground cover. The presence of penetrable surfaces such as vegetation, and tree canopies poses a challenge since the depth of penetration depends on a number of parameters such as the operating radar frequency, polarization, incident angle, as well as terrain structure. The dependence of the reconstructed topography on polarization may lead to the characterization of the ground cover. Simulation of interferometric measurements is useful for interpretation of the derived topography (B. Houshmand, Proceedings of URSI, 314, 1997). In this talk , time domain simulations for interferometric measurement for TE- and TM- polarization are presented. Time domain simulation includes the effects of the surface material property as well geometry comparable the radar signal wavelength (B. Houshmand, Proceedings of the URSI, 25, 1998). The IFSAR simulation is carried out in two steps. First, the forward scattering data is generated based on full wave analysis. Next, the electromagnetic information is inverted to generate surface topography. This inversion is based on the well known IFSAR processing technique which is composed of signal compression, and formation of an interferogram. The full wave forward scattering data is generated by the scattered-field formulation of the FDTD algorithm. The simulation is carried out by exciting the computational domain by a radar signal. The scattered field is then computed and translated to the receiving interferometric antennas using the time-domain Huygen's principle. The inversion process starts by compressing the time-domain data. The range compressed data from both receivers are then coregistered to form an interferogram. The resulting interferogram is then related to the ground topography using the radar imaging geometry. In this talk, the simulation results are compared with the C-band TM IFSAR derived topography, and the TE/TM SAR images at L-Band.

Houshmand, Bijan↗

A Framework for Extending the Science Traceability Matrix: Application to the Planned Europa Mission

One of the most critical functions of the systems engineering requirements process for a large multi-instrument science-driven space mission is to successfully communicate customer expectations into a comprehensive and traceable science requirements flowdown. These requirements are essential to communicating the constraints on the scope of the science investigations and clarifying how multiple instruments contribute to a given science goal. They also provide insight into how the science goals of the whole mission are affected by design choices. There is little specific guidance available on best practices for developing this science-driven flowdown. A unified Science Traceability Matrix (USTM) contains a significant amount of information that can be leveraged for that purpose, but the USTM was not designed to directly produce a complete science requirements flowdown. Thus, starting with the principles codified in a USTM, the authors propose a framework that directly maps into the requirements flowdown and supports broader systems engineering processes while retaining its meaning to the science team. This Science Traceability and Alignment Framework, or STAF, defines a set of common definitions and valid relationships to structure communication across the project. In addition, STAF populates a network of information that can be useful to support complex mission analysis activities such as fault protection. This work discusses the highest-level implementation of the STAF, the project-domain or P-STAF, which describes an approach to decomposing customer requirements into science requirements. The planned Europa Mission is used as a case study for the implementation of this framework and its potential benefits to a project.

Susca, Sara↗

Quantifying Structural Relationships of Metal-Binding Sites Suggests Origins of Biological Electron Transfer

Biological redox reactions drive planetary biogeochemical cycles. Using a novel, structure-guided sequence analysis of proteins, we explored the patterns of evolution of enzymes responsible for these reactions. Our analysis reveals that the folds that bind transition metal–containing ligands have similar structural geometry and amino acid sequences across the full diversity of proteins. Similarity across folds reflects the availability of key transition metals over geological time and strongly suggests that transition metal–ligand binding had a small number of common peptide origins. We observe that structures central to our similarity network come primarily from oxidoreductases, suggesting that ancestral peptides may have also facilitated electron transfer reactions. Last, our results reveal that the earliest biologically functional peptides were likely available before the assembly of fully functional protein domains over 3.8 billion years ago. Thus, life is a special, very complex form of motion of matter, but this form did not always exist, and it is not separated from inorganic nature by an impassable abyss; rather, it arose from inorganic nature as a new property in the process of evolution of the world. We must study the history of this evolution if we want to solve the problem of the origin of life.

Yana Bromberg↗

Active control for vibration suppression in a flexible beam using a modal domain optical fiber sensor

An account is given of the use of a modal-domain (MD) fiber-optic sensor as an active control system component for vibration suppression, whose output is proportional to the integral of the axial strain along the optical fiber. When an MD sensor is attached to, or embedded in, a flexible structure, it senses the strain in the structure along its gage length. On the basis of the present integration of the sensor model into a flexible-structure model, it becomes possible to design a control system with a dynamic compensator which adds damping to the low-order modes of the flexible structure. This modeling procedure has been experimentally validated.

Cox, D. E.↗

On a cost functional for H2/H(infinity) minimization

A cost functional is proposed and investigated which is motivated by minimizing the energy in a structure using only collocated feedback. Defined for an H(infinity)-norm bounded system, this cost functional also overbounds the H2 cost. Some properties of this cost functional are given, and preliminary results on the procedure for minimizing it are presented. The frequency domain cost functional is shown to have a time domain representation in terms of a Stackelberg non-zero sum differential game.

Macmartin, Douglas G.↗

Evaluation of piezoceramic actuators for control of aircraft interior noise

Results of an experiment to evaluate piezoceramic actuators as the control actuator for active control of interior noise in a large-scale fuselage model are presented. Control was demonstrated for tonal excitation using a time domain least mean squares algorithm. A maximum of four actuator channels and six error signals were used. The actuators were employed for control of noise at frequencies where interior cavity modes were the dominant response and for driven acoustic responses where a structure resonance was dominant. Global reductions of 9 to 12 dB were obtained for the cases examined. The most effective configuration of skin-mounted actuators was found to be a pure in-plane forcing function as opposed to a bending excitation. The frame-mounted actuators were found to be equally effective as the skin-mounted actuators. However, both configurations resulted in local regions of unacceptably high vibration response in the structure.

Silcox, Richard J.↗

Jupiter Thermospheric General Circulation Model (JTGCM): Global Structure and Dynamics Driven by Auroral and Joule Heating

A growing multispectral database plus recent Galileo descent measurements are being used to construct a self-consistent picture of the Jupiter thermosphere/ionosphere system. The proper characterization of Jupiter s upper atmosphere, embedded ionosphere, and auroral features requires the examination of underlying processes, including the feedbacks of energetics, neutral-ion dynamics, composition, and magnetospheric coupling. A fully 3-D Jupiter Thermospheric General Circulation Model (JTGCM) has been developed and exercised to address global temperatures, three-component neutral winds, and neutral-ion species distributions. The domain of this JTGCM extends from 20-microbar (capturing hydrocarbon cooling) to 1.0 x 10(exp -4) nbar (including aurora/Joule heating processes). The resulting JTGCM has been fully spun-up and integrated for greater than or equal to40 Jupiter rotations. Results from three JTGCM cases incorporating moderate auroral heating, ion drag, and moderate to strong Joule heating processes are presented. The neutral horizontal winds at ionospheric heights vary from 0.5 km/s to 1.2 km/s, atomic hydrogen is transported equatorward, and auroral exospheric temperatures range from approx.1200-1300 K to above 3000 K, depending on the magnitude of Joule heating. The equatorial temperature profiles from the JTGCM are compared with the measured temperature structure from the Galileo AS1 data set. The best fit to the Galileo data implies that the major energy source for maintaining the equatorial temperatures is due to dynamical heating induced by the low-latitude convergence of the high-latitude-driven thermospheric circulation. Overall, the Jupiter thermosphere/ionosphere system is highly variable and is shown to be strongly dependent on magnetospheric coupling which regulates Joule heating.

Bougher, S. W.↗

Global Weather States and Their Properties from Passive and Active Satellite Cloud Retrievals

In this study, the authors apply a clustering algorithm to International Satellite Cloud Climatology Project (ISCCP) cloud optical thickness-cloud top pressure histograms in order to derive weather states (WSs) for the global domain. The cloud property distribution within each WS is examined and the geographical variability of each WS is mapped. Once the global WSs are derived, a combination of CloudSat and Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) vertical cloud structure retrievals is used to derive the vertical distribution of the cloud field within each WS. Finally, the dynamic environment and the radiative signature of the WSs are derived and their variability is examined. The cluster analysis produces a comprehensive description of global atmospheric conditions through the derivation of 11 WSs, each representing a distinct cloud structure characterized by the horizontal distribution of cloud optical depth and cloud top pressure. Matching those distinct WSs with cloud vertical profiles derived from CloudSat and CALIPSO retrievals shows that the ISCCP WSs exhibit unique distributions of vertical layering that correspond well to the horizontal structure of cloud properties. Matching the derived WSs with vertical velocity measurements shows a normal progression in dynamic regime when moving from the most convective to the least convective WS. Time trend analysis of the WSs shows a sharp increase of the fair-weather WS in the 1990s and a flattening of that increase in the 2000s. The fact that the fair-weather WS is the one with the lowest cloud radiative cooling capability implies that this behavior has contributed excess radiative warming to the global radiative budget during the 1990s.

histograms↗

Verifying PLC Programs via Monitors: Extending the Integration of FRET and PLCverif

Verification of Programmable Logic Controller (PLC) programs requires reasoning about propositions qualified in terms of time. CERN’s PLCverif, an open-source tool for the analysis of safety-critical PLC systems, uses Linear Temporal Logic (LTL) for the specification of properties. Until now, PLCverif depended on third-party tools that accept LTL specifications to perform verification. However, our experience with industrial PLC programs shows that, to overcome analysis limitations, a wide range of techniques are needed to successfully verify complex properties. In this paper, we extend PLCverif to enable PLC program verification of pure-past LTL (PLTL) safety properties with assertion-based verification tools. To this end, we take an algorithm from the runtime-monitoring domain, apply it to bounded model checking of PLC programs, and implement it in PLCverif. We extend the integration of NASA’s Formal Requirements Elicitation Tool (FRET) into PLCverif to use PLTL properties generated with FRET. In addition, we leverage the program structure induced by the PLC scan-cycle for a state-space reduction. Finally, we expose the algorithm to a real-world case study of critical systems at CERN.

Formal verification↗

Structural stability augmentation system design using BODEDIRECT: A quick and accurate approach

A methodology is presented for a modal suppression control law design using flight test data instead of mathematical models to obtain the required gain and phase information about the flexible airplane. This approach is referred to as BODEDIRECT. The purpose of the BODEDIRECT program is to provide a method of analyzing the modal phase relationships measured directly from the airplane. These measurements can be achieved with a frequency sweep at the control surface input while measuring the outputs of interest. The measured Bode-models can be used directly for analysis in the frequency domain, and for control law design. Besides providing a more accurate representation for the system inputs and outputs of interest, this method is quick and relatively inexpensive. To date, the BODEDIRECT program has been tested and verified for computational integrity. Its capabilities include calculation of series, parallel and loop closure connections between Bode-model representations. System PSD, together with gain and phase margins of stability may be calculated for successive loop closures of multi-input/multi-output systems. Current plans include extensive flight testing to obtain a Bode-model representation of a commercial aircraft for design of a structural stability augmentation system.

Goslin, T. J.↗

Intergration of system identification and robust controller designs for flexible structures in space

An approach is developed using experimental data to identify a reduced-order model and its model error for a robust controller design. There are three steps involved in the approach. First, an approximately balanced model is identified using the Eigensystem Realization Algorithm, which is an identification algorithm. Second, the model error is calculated and described in frequency domain in terms of the H(infinity) norm. Third, a pole placement technique in combination with a H(infinity) control method is applied to design a controller for the considered system. A set experimental data from an existing setup, namely the Mini-Mast system, is used to illustrate and verify the approach.

Juang, Jer-Nan↗

Communications network design and costing model technical manual

This computer model provides the capability for analyzing long-haul trunking networks comprising a set of user-defined cities, traffic conditions, and tariff rates. Networks may consist of all terrestrial connectivity, all satellite connectivity, or a combination of terrestrial and satellite connectivity. Network solutions provide the least-cost routes between all cities, the least-cost network routing configuration, and terrestrial and satellite service cost totals. The CNDC model allows analyses involving three specific FCC-approved tariffs, which are uniquely structured and representative of most existing service connectivity and pricing philosophies. User-defined tariffs that can be variations of these three tariffs are accepted as input to the model and allow considerable flexibility in network problem specification. The resulting model extends the domain of network analysis from traditional fixed link cost (distance-sensitive) problems to more complex problems involving combinations of distance and traffic-sensitive tariffs.

Logan, K. P.↗

Linear baroclinic instability in the presence of large scale topography

The effect of a planetary-scale, wavenumber 2 topography on baroclinically active disturbances is investigated for a channel domain in a two-layer, quasi-geostrophic context. When the lower-layer zonal velocity is nonzero, the topography influences the disturbances by forcing a stationary wave, and the topography and the forced wave influence the growth rates and the spatial structures of the time-dependent solutions. The case of zero zonal velocity in the lower layer was also investigated, for which no forced wave exists. Asymptotic forms of the equations, valid when the topographic effect (governed by the ratio of the nondimensional topographic height to the rotational Froude number) is small, are used to obtain both the stationary and time-dependent solutions. The latter are also obtained using a numerical approach, in which is determined the eigenvalues and eigenfunctions of a matrix representing the dynamical equations. Agreement is good between the two approaches. Recent laboratory experiments with a baroclinic annulus in which there is a false bottom with wavenumber 2 topography, are used to select governing parameters. The simultaneous presence of a stationary forced wave of wavenumber 2 and a time-dependent baroclinic wave of wavenumber 4, which has wavenumber 2 and 6 sidebands due to the topography, yields a flow field that exhibits some principal features of the laboratory experiments. The position of the forced wave and the location of an excursion in latitude of the storm track show qualitative resemblance to those features observed in the atmosphere.

Reynolds, Nathaniel Dunton↗

A stable procedure to calculate the transient scattering by conducting surfaces of arbitrary shape

In this paper, a solution procedure is presented to obtain the transient scattering by arbitrarily shaped conducting objects directly in the time-domain using the marching-on-in-time method. The late-time oscillations are eliminated by a simple stabilization procedure which involves negligible amount of extra computation. Numerical results for surface current density and far-scattered fields are given for various structures and compared with other methods.

Vechinski, Douglas A.↗

Ishtar deformed belts: Evidence for deformation from below?

The mountain belts of Ishtar Terra are unique on Venus. Models for their formation include mantle upwelling, mantle downwelling, and horizontal convergence. The present forms of these models are too simple to predict surface strain, topography, or gravity. More detailed models will require specific constraints as imposed by geologic relations. In order to develop specific constraints for geodynamic models, we examine the geology of Ishtar Terra as viewed in Magellan SAR imagery in an attempt to interpret regional surface strain patterns. In this paper, we present geologic and structural relations that leads us to postulate that Ishtar deformed belts result from shear forces within the mantle acting on the lithosphere, and not by horizontal forces from colliding plates. We propose that the surface strains result from differential strain and displacement of domains within the upper mantle, and that further analysis of Ishtar deformation may allow us to identify individual domains within the mantle, and to constrain displacement trajectories between domains.

Hansen, V. L.↗

Interaction of a Boundary Layer with a Turbulent Wake

The objective of this grant was to study the transition mechanisms on a flat-plate boundary layer interacting with the wake of a bluff body. This is a simplified configuration presented and designed to exemplify the phenomena that occur in multi-element airfoils, in which the wake of an upstream element impinges on a downstream one. Some experimental data is available for this configuration at various Reynolds numbers. The first task carried out was the implementation and validation of the immersed-boundary method. This was achieved by performing calculations of the flow over a cylinder at low and moderate Reynolds numbers. The low-Reynolds number results are discussed, which is enclosed as Appendix A. The high-Reynolds number results are presented in a paper in preparation for the Journal of Fluid Mechanics. We performed calculations of the wake-boundary-layer interaction at two Reynolds numbers, Re approximately equal to 385 and 1155. The first case is discussed and a comparison of the two calculations is reported. The simulations indicate that at the lower Reynolds number the boundary layer is buffeted by the unsteady Karman vortex street shed by the cylinder. This is shown: long streaky structures appear in the boundary layer in correspondence of the three-dimensionalities in the rollers. The fluctuations, however, cannot be self-sustained due to the low Reynolds-number, and the flow does not reach a turbulent state within the computational domain. In contrast, in the higher Reynolds-number case, boundary-layer fluctuations persist after the wake has decayed (due, in part, to the higher values of the local Reynolds number Re achieved in this case); some evidence could be observed that a self-sustaining turbulence generation cycle was beginning to be established. A third simulation was subsequently carried out at a higher Reynolds number, Re=3900. This calculation gave results similar to those of the Re=l155 case. Turbulence was established at fairly low Reynolds number, as a consequence of the high level of the free-stream perturbation. An instantaneous flow visualization for that case is shown. A detailed examination of flow statistics in the transitional and turbulent regions, including the evolution of the turbulent kinetic energy (TKE) budget and frequency spectra showed the formation and evolution of turbulent spots characteristic of the bypass transition mechanism. It was also observed that the turbulent eddies achieved an equilibrium, fully developed turbulent states first, as evidenced by the early agreement achieved by the terms in the TKE budget with those observed in turbulent flows. Once a turbulent Reynolds stress profile had been established, the velocity profile began to resemble a turbulent one, first in the inner region and later in the outer region of the wall layer. An extensive comparison of the three cases, including budgets, mean velocity and Reynolds stress profiles and flow visualization, is included. The results obtained are also presented.

Piomelli, Ugo↗

Ares I-X In-Flight Modal Identification

Operational modal analysis is a procedure that allows the extraction of modal parameters of a structure in its operating environment. It is based on the idealized premise that input to the structure is white noise. In some cases, when free decay responses are corrupted by unmeasured random disturbances, the response data can be processed into cross-correlation functions that approximate free decay responses. Modal parameters can be computed from these functions by time domain identification methods such as the Eigenvalue Realization Algorithm (ERA). The extracted modal parameters have the same characteristics as impulse response functions of the original system. Operational modal analysis is performed on Ares I-X in-flight data. Since the dynamic system is not stationary due to propellant mass loss, modal identification is only possible by analyzing the system as a series of linearized models over short periods of time via a sliding time-window of short time intervals. A time-domain zooming technique was also employed to enhance the modal parameter extraction. Results of this study demonstrate that free-decay time domain modal identification methods can be successfully employed for in-flight launch vehicle modal extraction.

Bartkowicz, Theodore J.↗