Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interface problems”

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,333 records · Page 74

Interface conditions for domain decomposition with radical grid refinement

Interface conditions for coupling the domains in a physically motivated domain decomposition method are discussed. The domain decomposition is based on an asymptotic-induced method for the numerical solution of hyperbolic conservation laws with small viscosity. The method consists of multiple stages. The first stage is to obtain a first approximation using a first-order method, such as the Godunov scheme. Subsequent stages of the method involve solving internal-layer problem via a domain decomposition. The method is derived and justified via singular perturbation techniques.

Scroggs, Jeffrey S.↗

GVIPS Models and Software

Two reports discuss, respectively, (1) the generalized viscoplasticity with potential structure (GVIPS) class of mathematical models and (2) the Constitutive Material Parameter Estimator (COMPARE) computer program. GVIPS models are constructed within a thermodynamics- and potential-based theoretical framework, wherein one uses internal state variables and derives constitutive equations for both the reversible (elastic) and the irreversible (viscoplastic) behaviors of materials. Because of the underlying potential structure, GVIPS models not only capture a variety of material behaviors but also are very computationally efficient. COMPARE comprises (1) an analysis core and (2) a C++-language subprogram that implements a Windows-based graphical user interface (GUI) for controlling the core. The GUI relieves the user of the sometimes tedious task of preparing data for the analysis core, freeing the user to concentrate on the task of fitting experimental data and ultimately obtaining a set of material parameters. The analysis core consists of three modules: one for GVIPS material models, an analysis module containing a specialized finite-element solution algorithm, and an optimization module. COMPARE solves the problem of finding GVIPS material parameters in the manner of a design-optimization problem in which the parameters are the design variables.

Arnold, Steven M.↗

Progress in Computational Simulation of Earthquakes

GeoFEST(P) is a computer program written for use in the QuakeSim project, which is devoted to development and improvement of means of computational simulation of earthquakes. GeoFEST(P) models interacting earthquake fault systems from the fault-nucleation to the tectonic scale. The development of GeoFEST( P) has involved coupling of two programs: GeoFEST and the Pyramid Adaptive Mesh Refinement Library. GeoFEST is a message-passing-interface-parallel code that utilizes a finite-element technique to simulate evolution of stress, fault slip, and plastic/elastic deformation in realistic materials like those of faulted regions of the crust of the Earth. The products of such simulations are synthetic observable time-dependent surface deformations on time scales from days to decades. Pyramid Adaptive Mesh Refinement Library is a software library that facilitates the generation of computational meshes for solving physical problems. In an application of GeoFEST(P), a computational grid can be dynamically adapted as stress grows on a fault. Simulations on workstations using a few tens of thousands of stress and displacement finite elements can now be expanded to multiple millions of elements with greater than 98-percent scaled efficiency on over many hundreds of parallel processors (see figure).

Donnellan, Andrea↗

CEA2022: A Modernization of NASA Glenn’s Software CEA (Chemical Equilibrium with Applications)

The software program “Chemical Equilibrium with Applications” (CEA) is used to solve chemical equilibrium, and compute thermodynamic and transport properties of the resulting mixture, and also has special solvers dedicated to rocket, shock, and detonation problems. We have recently completed a full re-write of CEA with modernization and improvements, called “CEA2022”. In this paper, we will give an overview of CEA2022’s features, and discuss some of the fundamental equations used by CEA2022, as well as the fundamental assumptions, in order to provide users with a complete understanding of the software’s methodology. Several enhancements have been made to the software, including modern software development practices, interface improvements, and additional features. The feature enhancements include: running cases in parallel with thread safe solves, thermodynamic and transport database updates, and allowing for negative and inert reactants. In terms of interface improvements, we have made CEA a reusable library by adding APIs for multiple languages, including Python, Matlab, Excel, Fortran, and C. The subroutine interface allows for integration with other applications, including flow-solver integration (i.e. with CFD). We also compare results between CEA2022 and the previous version (CEA2) as a validation of the new software.

chemical equilibrium↗

Using Immersive Virtual Environments for Certification

Immersive virtual environments (VEs) technology has matured to the point where it can be utilized as a scientific and engineering problem solving tool. In particular, VEs are starting to be used to design and evaluate safety-critical systems that involve human operators, such as flight and driving simulators, complex machinery training, and emergency rescue strategies.

virtual reality software certification human-compu↗

Adhesion, friction, wear, and lubrication research by modern surface science techniques.

The field of surface science has undergone intense revitalization with the introduction of low-energy electron diffraction, Auger electron spectroscopy, ellipsometry, and other surface analytical techniques which have been sophisticated within the last decade. These developments have permitted submono- and monolayer structure analysis as well as chemical identification and quantitative analysis. The application of a number of these techniques to the solution of problems in the fields of friction, lubrication, and wear are examined in detail for the particular case of iron; and in general to illustrate how the accumulation of pure data will contribute toward the establishment of physiochemical concepts which are required to understand the mechanisms that are operational in friction systems. In the case of iron, LEED, Auger and microcontact studies have established that hydrogen and light-saturated organic vapors do not establish interfaces which prevent iron from welding, whereas oxygen and some oxygen and sulfur compounds do reduce welding as well as the coefficient of friction. Interpretation of these data suggests a mechanism of sulfur interaction in lubricating systems.

Keller, D. V., Jr.↗

Stress intensity factors of composite orthotropic plates containing periodic buffer strips

The fracture problem of laminated plates which consist of bonded orthotropic layers is studied. The fields equations for an elastic orthotropic body are transformed to give the displacement and stress expressions for each layer or strip. The unknown functions in these expressions are found by satisfying the remaining boundary and continuity conditions. A system of singular integral equations is obtained from the mixed boundary conditions. The singular behavior around the crack tip and at the bimaterial interface is studied. The stress intensity factors are computed for various material combinations and various crack geometries. The results are discussed and are compared with those for isotropic materials.

Delale, F.↗

Modeling of Thermal Performance of Multiphase Nuclear Fuel Cell Under Variable Gravity Conditions

A unique numerical method has been developed to model the dynamic processes of bulk evaporation and condensation processes, associated with internal heat generation and natural convection under different gravity levels. The internal energy formulation, for the bulk liquid-vapor phase change problems in an encapsulated container, was employed. The equations, governing the conservation of mass, momentum and energy for both phases involved in phase change, were solved. The thermal performance of a multiphase uranium tetra-fluoride fuel element under zero gravity, micro-gravity and normal gravity conditions has been investigated. The modeling yielded results including the evolution of the bulk liquid-vapor phase change process, the evolution of the liquid-vapor interface, the formation and development of the liquid film covering the side wall surface, the temperature distribution and the convection flow field in the fuel element. The strong dependence of the thermal performance of such multiphase nuclear fuel cell on the gravity condition has been revealed. Under all three gravity conditions, 0-g, 10(exp -3)-g, and 1-g, the liquid film is formed and covers the entire side wall. The liquid film covering the side wall is more isothermalized at the wall surface, which can prevent the side wall from being over-heated. As the gravity increases, the liquid film is thinner, the temperature gradient is larger across the liquid film and smaller across the vapor phase. This investigation provides valuable information about the thermal performance of multi-phase nuclear fuel element for the potential space and ground applications.

Ding, Z.↗

PLUM: Parallel Load Balancing for Unstructured Adaptive Meshes

Dynamic mesh adaption on unstructured grids is a powerful tool for computing large-scale problems that require grid modifications to efficiently resolve solution features. Unfortunately, an efficient parallel implementation is difficult to achieve, primarily due to the load imbalance created by the dynamically-changing nonuniform grid. To address this problem, we have developed PLUM, an automatic portable framework for performing adaptive large-scale numerical computations in a message-passing environment. First, we present an efficient parallel implementation of a tetrahedral mesh adaption scheme. Extremely promising parallel performance is achieved for various refinement and coarsening strategies on a realistic-sized domain. Next we describe PLUM, a novel method for dynamically balancing the processor workloads in adaptive grid computations. This research includes interfacing the parallel mesh adaption procedure based on actual flow solutions to a data remapping module, and incorporating an efficient parallel mesh repartitioner. A significant runtime improvement is achieved by observing that data movement for a refinement step should be performed after the edge-marking phase but before the actual subdivision. We also present optimal and heuristic remapping cost metrics that can accurately predict the total overhead for data redistribution. Several experiments are performed to verify the effectiveness of PLUM on sequences of dynamically adapted unstructured grids. Portability is demonstrated by presenting results on the two vastly different architectures of the SP2 and the Origin2OOO. Additionally, we evaluate the performance of five state-of-the-art partitioning algorithms that can be used within PLUM. It is shown that for certain classes of unsteady adaption, globally repartitioning the computational mesh produces higher quality results than diffusive repartitioning schemes. We also demonstrate that a coarse starting mesh produces high quality load balancing, at a fraction of the cost required a fine initial mesh. Results indicate that our parallel load balancing strategy will remain viable on large numbers of processors.

Oliker, Leonid↗

Hydrodynamic Instability and Thermal Coupling in a Dynamic Model of Liquid-Propellant Combustion

For liquid-propellant combustion, the Landau/Levich hydrodynamic models have been combined and extended to account for a dynamic dependence of the burning rate on the local pressure and temperature fields. Analysis of these extended models is greatly facilitated by exploiting the realistic smallness of the gas-to-liquid density ratio rho. Neglecting thermal coupling effects, an asymptotic expression was then derived for the cellular stability boundary A(sub p)(k) where A(sub p) is the pressure sensitivity of the burning rate and k is the disturbance wavenumber. The results explicitly indicate the stabilizing effects of gravity on long-wave disturbances, and those of viscosity and surface tension on short-wave perturbations, and the instability associated with intermediate wavenumbers for critical negative values of A(sub p). In the limit of weak gravity, hydrodynamic instability in liquid-propellant combustion becomes a long-wave, instability phenomenon, whereas at normal gravity, this instability is first manifested through O(1) wavenumbers. In addition, surface tension and viscosity (both liquid and gas) each produce comparable effects in the large-wavenumber regime, thereby providing important modifications to the previous analyses in which one or more of these effects was neglected. For A(sub p)= O, the Landau/Levich results are recovered in appropriate limiting cases, although this typically corresponds to a hydrodynamically unstable parameter regime for p << 1. In addition to the classical cellular form of hydrodynamic stability, there exists a pulsating form corresponding to the loss of stability of steady, planar burning to time-dependent perturbations. This occurs for negative values of the parameter A(sub p), and is thus absent from the original Landau/Levich models. In the extended model, however, there exists a stable band of negative pressure sensitivities bounded above by the Landau type of instability, and below by this pulsating form of hydrodynamic instability. Indeed, nonsteady modes of combustion have been observed at low pressures in hydroxylammonium nitrate (HAN)-based liquid propellants, which often exhibit negative pressure sensitivities. While nonsteady combustion may correspond to secondary and higher-order bifurcations above the cellular boundary, it may also be a manifestation of this pulsating type of hydrodynamic instability. In the present work, a nonzero temperature sensitivity is incorporated into our previous asymptotic analyses. This entails a coupling of the energy equation to the previous purely hydrodynamic problem, and leads to a significant modification of the pulsating boundary such that, for sufficiently large values of the temperature-sensitivity parameter, liquid-propellant combustion can become intrinsically unstable to this alternative form of hydrodynamic instability. For simplicity, further attention is confined here to the inviscid version of the problem since, despite the fact that viscous and surface-tension effects are comparable, the qualitative nature of the cellular boundary remains preserved in the zero-viscosity limit, as does the existence of the pulsating boundary. The mathematical model adopts the classical assumption that there is no distributed reaction in either the liquid or gas phases, but now the reaction sheet, representing either a pyrolysis reaction or an exothermic decomposition at the liquid/gas interface, is assumed to depend on local conditions there.

Margolis, S. B.↗

Melt Motion Due to Peltier Marking During Bridgman Crystal Growth with an Axial Magnetic Field

This paper treats a liquid-metal flow inside an electrically insulating cylinder with electrically conducting solids above and below the liquid region. There is a uniform axial magnetic field, and there is an electric current through the liquid and both solids. Since the lower liquid-solid interface is concave into the solid and since the liquid is a better electrical conductor than the adjacent solid, the electric current is locally concentrated near the centerline. The return to a uniform current distribution involves a radial electric current which interacts with the axial magnetic field to drive an azimuthal flow. The axial variation of the centrifugal force due to the azimuthal velocity drives a meridional circulation with radial and axial velocities. This problem models the effects of Peltier marking during the vertical Bridgman growth of semiconductor crystals with an externally applied magnetic field, where the meridional circulation due to the Peltier Current may produce important mixing in the molten semiconductor.

Sellers, C. C.↗

A Self-Stabilizing Hybrid Fault-Tolerant Synchronization Protocol

This paper presents a strategy for solving the Byzantine general problem for self-stabilizing a fully connected network from an arbitrary state and in the presence of any number of faults with various severities including any number of arbitrary (Byzantine) faulty nodes. The strategy consists of two parts: first, converting Byzantine faults into symmetric faults, and second, using a proven symmetric-fault tolerant algorithm to solve the general case of the problem. A protocol (algorithm) is also present that tolerates symmetric faults, provided that there are more good nodes than faulty ones. The solution applies to realizable systems, while allowing for differences in the network elements, provided that the number of arbitrary faults is not more than a third of the network size. The only constraint on the behavior of a node is that the interactions with other nodes are restricted to defined links and interfaces. The solution does not rely on assumptions about the initial state of the system and no central clock nor centrally generated signal, pulse, or message is used. Nodes are anonymous, i.e., they do not have unique identities. A mechanical verification of a proposed protocol is also present. A bounded model of the protocol is verified using the Symbolic Model Verifier (SMV). The model checking effort is focused on verifying correctness of the bounded model of the protocol as well as confirming claims of determinism and linear convergence with respect to the self-stabilization period.

Malekpour, Mahyar R.↗

End-to-End Trade-Space Analysis for Designing Constellation

Multipoint measurement missions can provide a significant advancement in science return and this science interest coupled with as many recent technological advances are driving a growing trend in exploring distributed architectures for future NASA missions. Distributed Spacecraft Missions (DSMs) leverage multiple spacecraft to achieve one or more common goals. In particular, a constellation is the most general form of DSM with two or more spacecraft placed into specific orbit(s) for the purpose of serving a common objective (e.g., CYGNSS). Because a DSM architectural trade-space includes both monolithic and distributed design variables, DSM optimization is a large and complex problem with multiple conflicting objectives. Over the last two years, our team has been developing a Trade-space Analysis Tool for Constellations (TAT-C), implemented in common programming languages for pre-Phase A constellation mission analysis. By evaluating alternative mission architectures, TAT-C seeks to minimize cost and maximize performance for pre-defined science goals. This presentation will describe the overall architecture of TAT-C including: a User Interface (UI) at several levels of details and user expertise; Trade-space Search Requests that are created from the Science requirements gathered by the UI and validated by a Knowledge Base; a Knowledge Base to compare the current requests to prior mission concepts to potentially prune the trade-space; a Trade-space Search Iterator which, with inputs from the Knowledge Base, and, in collaboration with the Orbit & Coverage, Reduction & Metrics, and Cost& Risk modules, generates multiple potential architectures and their associated characteristics. TAT-C leverages the use of the Goddard Mission Analysis Tool (GMAT) to compute coverage and ancillary data, modeling orbits to balance accuracy and performance. The current version includes uniform and non-uniform Walker constellations as well as Ad-Hoc and precessing constellations, and its cost model represents an aggregate model consisting of Cost Estimating Relationships (CERs) from widely accepted models. The current GUI automatically generates graphics representing metrics such as average revisit time or coverage as a function of cost. The end-to-end system will be demonstrated as part of the presentation.

Mission Desig↗

Identification of structural interface characteristics using component mode synthesis

The inability to adequately model connections has limited the ability to predict overall system dynamic response. Connections between structural components are often mechanically complex and difficult to accurately model analytically. Improved analytical models for connections are needed to improve system dynamic predictions. This study explores combining Component Mode synthesis methods for coupling structural components with Parameter Identification procedures for improving the analytical modeling of the connections. Improvements in the connection properties are computed in terms of physical parameters so the physical characteristics of the connections can be better understood, in addition to providing improved input for the system model. Two sample problems, one utilizing simulated data, the other using experimental data from a rotor dynamic test rig are presented.

Huckelbridge, A. A.↗

Identification of structural interface characteristics using component mode synthesis

The inability to adequately model connections has limited the ability to predict overall system dynamic response. Connections between structural components are often mechanically complex and difficult to accurately model analytically. Improved analytical models for connections are needed to improve system dynamic predictions. This study explores combining Component Mode synthesis methods for coupling structural components with Parameter Identification procedures for improving the analytical modeling of the connections. Improvements in the connection properties are computed in terms of physical parameters so the physical characteristics of the connections can be better understood, in addition to providing improved input for the system model. Two sample problems, one utilizing simulated data, the other using experimental data from a rotor dynamic test rig, are presented.

Huckelbridge, A. A.↗

Identification of structural interface characteristics using component mode synthesis

The inability to adequately model connections has limited the ability to predict overall system dynamic response. Connections between structural components are often mechanically complex and difficult to accurately model analytically. Improved analytical models for connections are needed to improve system dynamic predictions. This study explores combining Component Mode synthesis methods for coupling structural components with Parameter Identification procedures for improving the analytical modeling of the connections. Improvements in the connection properties are computed in terms of physical parameters so the physical characteristics of the connections can be better understood, in addition to providing improved input for the system model. Two sample problems, one utilizing simulated data, the other using experimental data from a rotor dynamic test rig, are presented.

Huckelbridge, A. A.↗

Digital adaptive flight control design using single stage model following indices

Simple mechanical linkages are often unable to cope with many control problems associated with high-performance aircraft. This has led to the development of digital fly-by-wire control systems and in particular digital adaptive controllers that can be efficiently adjusted during system operation. To this effect, a control law has been derived based upon the minimization of a single-stage weighted combination of control energy and the squared error between the states of a linear plant and model. This control logic is interfaced with an on-line weighted least-squares estimator and a Kalman state filter. The utility of the resultant control system is illustrated by its application to the linearized dynamics of a typical fighter aircraft.

Alag, G.↗

Object impedance control for cooperative manipulation - Theory and experimental results

The dynamic control module being developed in the Dynamic and Strategic Control of Cooperative Manipulators (DASCCOM) project at the Stanford University Aerospace Robotics Laboratory is described. First, the cooperative manipulation problem is analyzed from a systems perspective, and the desirable features of a control system for cooperative manipulation are discussed. Next, a control policy is developed that enforces a controlled impedance not of the individual arm endpoints, but of the manipulated object itself. A parallel implementation for a multiprocessor system is presented. The controller fully compensates for the system dynamics and directly controls the object internal forces. Most importantly, it presents a simple, powerful, intuitive interface to the strategic controller. Experimental results for a dual two-link arm robotic system are presented to verify the controllers performance, for both free-motion slews and environmental contact.

Schneider, Stanley A.↗