Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “state dependent variable”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

The Winds of B Supergiants

We present the most suitable data sets available in the International Ultraviolet Explorer (IUE) archive for the study of time-dependent stellar winds in early B supergiants. The UV line profile variability in 11 B0 to B3 stars is analyzed, compared and discussed, based on 16 separate data sets comprising over 600 homogeneously reduced high-resolution spectrograms. The targets include 'normal' stars with moderate rotation rates and examples of rapid rotators. A gallery of grey-scale images (dynamic spectra) is presented, which demonstrates the richness and range of wind variability and highlights different structures in the winds of these stars. This work emphasizes the suitability of B supergiants for wind studies, under-pinned by the fact that they exhibit unsaturated wind lines for a wide range of ionization. The wind activity of B supergiants is substantial and has highly varied characteristics. The variability evident in individual stars is classified and described in terms of discrete absorption components, spontaneous absorption, bowed structures, recurrence, and ionization variability and stratification. Similar structures can occur in stars of different fundamental parameters but also different structures may occur in the same star at a given epoch. We discuss the physical phenomena that may be associated with the spectral signatures, and highlight the challenges that these phenomena present to theoretical studies of time-dependent outflows in massive stars. In addition, SEI line-synthesis modelling of the UV wind lines is used to provide further information about the state of the winds in our program stars. Typically the range, implied by the line profile variability, in the product of mass-loss rate and ion fraction (M qi) is a factor of approximately 1.5, when integrated between 0.2 and 0.9 v infinity; it it can however be several times larger over localized velocity regions. At a given effective temperature the mean relative ion ratios can differ by a factor of 5. The general excess in predicted (forward-scattered) emission in the low velocity regime is discussed in turns of structured outflows. Mean ion fractions are estimated over the B0 to B1 spectral classes, and trends in the ionic ratios as a function of wind velocity are described. The low values obtained for the ion fractions of UV resonance lines may reflect the role of clumping in the wind.

Massa, Derck↗
NASA NTRS2002

The Winds of B Supergiants

We present the most suitable data sets available in the International Ultraviolet Explorer (IUE) archive for the study of time-dependent stellar winds in early B supergiants. The UV line profile variability in 11 B0 to B3 stars is analyzed, compared and discussed, based on 16 separate data sets comprising over 600 homogeneously reduced high-resolution spectrograms. The targets include 'normal' stars with moderate rotation rates and examples of rapid rotators. A gallery of grey-scale images (dynamic spectra) is presented, which demonstrates the richness and range of wind variability and highlights different structures in the winds of these stars. This work emphasises the suitability of B supergiants for wind studies, under-pinned by the fact that they exhibit unsaturated wind lines for a wide range of ionization. The wind activity of B supergiants is substantial and has highly varied characteristics. The variability evident in individual stars is classified and described in terms of discrete absorption components, spontaneous absorption, bowed structures, recurrence, and ionization variability and stratification. Similar structures can occur in stars of different fundamental parameters, but also different structures may occur in the same star at a given epoch. We discuss the physical phenomena that may be associated with the spectral signatures, and highlight the challenges that these phenomena present to theoretical studies of time-dependent outflows in massive stars. In addition, SEI line-synthesis modelling of the UV wind lines is used to provide further information about the state of the winds in our program stars. Typically the range, implied by the line profile variability, in the product of mass-loss rate and ion fraction (M (dot) q(sub i)) is a factor of approximately 1.5, when integrated between 0.2 and 0.9 v infinity; it can however be several times larger over localized velocity regions. At a given effective temperature the mean relative ion ratios can differ by a factor of 5. The general excess in predicted (forward-scattered) emission in the low velocity regime is discussed in terms of structured outflows. Mean ion fractions are estimated over the B0 to B1 spectral classes, and trends in the ionic ratios as a function of wind velocity are described. The low values obtained for the ion fractions of UV resonance lines may reflect the role of clumping in the wind.

Massa, D.↗
NASA NTRS1991

Leading-edge receptivity to a vortical freestream disturbance: A numerical analysis

The receptivity to freestream vorticity of the boundary layer over a flat plate with an elliptic leading edge is investigated numerically. The flow is simulated by solving the incompressible Navier-Stokes system in general curvilinear coordinates with the vorticity and stream function as dependent variables. A finite-difference scheme which is second-order accurate in both space and time is used. As a first step, the steady basic-state solution is computed. Then a small amplitude vortical disturbance is introduced at the upstream boundary and the governing equations are solved time-accurately to evaluate the spatial and temporal growth of the perturbations leading to instability waves (Tollmien-Schlichting waves) inside the boundary layer. Preliminary results for a symmetric, 2-D disturbance reveal the presence of Tollmien-Schlichting waves aft of the flat-plate/ellipse juncture.

Buter, Thomas A.↗
NASA NTRS1998

Experimental Investigation of Organic Synthesis in Hydrothermal Environments

Seafloor hydrothermal systems may be the most likely locations on the early Earth for the emergence of life. Because of the disequilibrium inherent in such dynamic, mixing environments, abundant chemical energy would have been available for formation of the building blocks of life. In addition, theoretical studies suggest that organic compounds in these conditions would reach metastable states, due to kinetic barriers to the formation of stable equilibrium products (CO2 and methane). The speciation of organic carbon in metastable states is highly dependent on the oxidation state, pH, temperature, pressure and bulk composition of the system. The goal of our research is to investigate the effects of a number external variables on the formation, transformation, and stability of organic compounds at hydrothermal conditions. We have begun experimental work to attempt to control the oxidation state of simulated hydrothermal systems by using buffers composed of mineral powders and gas mixtures. We are also beginning to test the stability of organic compounds under these conditions. The experiments are being performed using the hydrothermal bomb apparatus at the U.S. Geological Survey in Menlo Park, CA and the supercritical water oxidizer (SCWO) at NASA Ames Research Center in Moffet Field, CA. The amino acids decomposed rapidly. Even after the approximately 15 minutes between addition of the amino acids and the first sampling, no amino acids were detected in the PPM system by GC- MS, while in the FeFeO system the amino acids were present at a level of less than 50% of original. Carboxylic acids, ammonia, and CO2 were the main products, along with some unidentified compounds. The ratios of carboxylic acids and concentrations of other products seem to have remained stable during the experiments, consistent with observations of other metastable systems and theoretical predictions.

Shock, Everett L.↗
NASA NTRS1998

Experimental Investigation of Organic Synthesis in Hydrothermal Environments

The results of the investigation were presented at a Astrobiology Institute General Meeting. Seafloor hydrothermal systems may be the most likely locations on the early Earth for the emergence of life. Because of the disequilibrium inherent in such dynamic, mixing environments, abundant chemical energy would have been available for formation of the building blocks of life. In addition, theoretical studies suggest that organic compounds in these conditions would reach metastable states, due to kinetic barriers to the formation of stable equilibrium products (CO2 and methane). The speciation of organic carbon in metastable states is highly dependent on the oxidation state, pH, temperature, pressure and bulk composition of the system. The goal of our research is to investigate the effects of a number external variables on the formation, transformation, and stability of organic compounds at hydrothermal conditions. We have begun experimental work to attempt to control the oxidation state of simulated hydrothermal systems by using buffers composed of mineral powders and gas mixtures. We are also beginning to test the stability of organic compounds under these conditions.

Shock, Everett L.↗
NASA NTRS2008

In-Field Diffuse Ultraviolet Spectroscopy and Imaging of the Stardust Sample Return Capsule

In-field diffuse Ultraviolet (UV) spectroscopy and imaging systems were developed for the purposes of evaluating the surface chemical composition of spacecraft thermal control coatings and materials. The investigation of these systems and the compilation of an associated UV reflectance and luminescence database were conducted using the Stardust Sample Return Capsule (SRC), located at the Johnson Space Center. Spectral responses of the surfaces of the Stardust forebody and aftbody in both reflectance and fluorescence modes were examined post-flight. In this paper, we report on two primary findings of in-field diffuse UV spectroscopy and imaging: (1) deduction of the thermal history of thermal control coatings of the forebody and (2) bond line variations in the aftbody. In the forebody, the thermal history of thermal control coatings may be deduced from the presence of particular semiconducting defect states associated with ZnO, a common emissivity constituent in thermal control coatings. A spatial dependence of this history was mapped for these regions. In the aftbody, luminescing defect states, associated with Si and SiO2 color centers were found along regions of bond variability.

Pugel, D. Elizabeth↗
NASA NTRS1992

Volume 2: Explicit, multistage upwind schemes for Euler and Navier-Stokes equations

The objective of this study was to develop a high-resolution-explicit-multi-block numerical algorithm, suitable for efficient computation of the three-dimensional, time-dependent Euler and Navier-Stokes equations. The resulting algorithm has employed a finite volume approach, using monotonic upstream schemes for conservation laws (MUSCL)-type differencing to obtain state variables at cell interface. Variable interpolations were written in the k-scheme formulation. Inviscid fluxes were calculated via Roe's flux-difference splitting, and van Leer's flux-vector splitting techniques, which are considered state of the art. The viscous terms were discretized using a second-order, central-difference operator. Two classes of explicit time integration has been investigated for solving the compressible inviscid/viscous flow problems--two-state predictor-corrector schemes, and multistage time-stepping schemes. The coefficients of the multistage time-stepping schemes have been modified successfully to achieve better performance with upwind differencing. A technique was developed to optimize the coefficients for good high-frequency damping at relatively high CFL numbers. Local time-stepping, implicit residual smoothing, and multigrid procedure were added to the explicit time stepping scheme to accelerate convergence to steady-state. The developed algorithm was implemented successfully in a multi-block code, which provides complete topological and geometric flexibility. The only requirement is C degree continuity of the grid across the block interface. The algorithm has been validated on a diverse set of three-dimensional test cases of increasing complexity. The cases studied were: (1) supersonic corner flow; (2) supersonic plume flow; (3) laminar and turbulent flow over a flat plate; (4) transonic flow over an ONERA M6 wing; and (5) unsteady flow of a compressible jet impinging on a ground plane (with and without cross flow). The emphasis of the test cases was validation of code, and assessment of performance, as well as demonstration of flexibility.

Elmiligui, Alaa↗
NASA NTRS1998

A Second Law Based Unstructured Finite Volume Procedure for Generalized Flow Simulation

An unstructured finite volume procedure has been developed for steady and transient thermo-fluid dynamic analysis of fluid systems and components. The procedure is applicable for a flow network consisting of pipes and various fittings where flow is assumed to be one dimensional. It can also be used to simulate flow in a component by modeling a multi-dimensional flow using the same numerical scheme. The flow domain is discretized into a number of interconnected control volumes located arbitrarily in space. The conservation equations for each control volume account for the transport of mass, momentum and entropy from the neighboring control volumes. In addition, they also include the sources of each conserved variable and time dependent terms. The source term of entropy equation contains entropy generation due to heat transfer and fluid friction. Thermodynamic properties are computed from the equation of state of a real fluid. The system of equations is solved by a hybrid numerical method which is a combination of simultaneous Newton-Raphson and successive substitution schemes. The paper also describes the application and verification of the procedure by comparing its predictions with the analytical and numerical solution of several benchmark problems.

Majumdar, Alok↗
NASA NTRSDate not supplied

Understanding Resilience Optimization Architectures With an Optimization Problem Repository

Optimizing a system’s resilience can be challenging, especially when it involves considering both the inherent resilience of a robust design and the active resilience of a health management system to a set of computationally-expensive hazard simulations. While prior work has developed specialized architectures to effectively and efficiently solve combined design and resilience optimization problems, the comparison of these architectures has been limited to a single case study. To further study resilience optimization formulations, this work develops a problem repository which includes previously-developed resilience optimization problems and additional problems presented in this work: a notional system resilience model, a pandemic response model, and a cooling tank hazard prevention model. This work then uses models in the repository at large to understand the characteristics of resilience optimization problems and study the applicability of optimization architectures and decomposition strategies. Based on the comparisons in the repository, applying an optimization architecture effectively requires understanding the alignment and coupling relationships between the design and resilience models, as well as the efficiency characteristics of the algorithms. While alignment determines the necessity of a surrogate of resilience cost in the upper-level design problem, coupling determines the overall applicability of a sequential, alternating, or bilevel structure. Additionally, the application of decomposition strategies is dependent on there being limited interactions between variable sets, which often does not hold when a resilience policy is parameterized in terms of actions to take in hazardous model states rather than specific given scenarios.

Resilience↗
NASA NTRS1975

Pressure- and buoyancy-driven thermal convection in a rectangular enclosure

Results are presented for unsteady laminar thermal convection in compressible fluids at various reduced levels of gravity in a rectangular enclosure which is heated on one side and cooled on the opposite side. The results were obtained by solving numerically the equations of conservation for a viscous, compressible, heat-conducting, ideal gas in the presence of a gravitational body force. The formulation differs from the Boussinesq simplification in that the effects of variable density are completely retained. A conservative, explicit, time-dependent, finite-difference technique was used and good agreement was found for the limited cases where direct comparison with previous investigations was possible. The solutions show that the thermally induced motion is acoustic in nature at low levels of gravity and that the unsteady-state rate of heat transfer is thereby greatly enhanced relative to pure conduction. The nonlinear variable density profile skews the streamlines towards the cooler walls but is shown to have little effect on the steady-state isotherms.

Spradley, L. W.↗
NASA NTRS1993

Explicit 3D continuum fracture modeling with smooth particle hydrodynamics

Impact phenomena shaped our solar system. As usual for most solar system processes, the scales are far different than we can address directly in the laboratory. Impact velocities are often much higher than we can achieve, sizes are often vastly larger, and most impacts take place in an environment where the only gravitational force is the mutual pull of the impactors. The Smooth Particle Hydrodynamics (SPH) technique has been applied in the past to the simulations of giant impacts. In these simulations, the colliding objects were so massive (at least a sizeable fraction of the Earth's mass) that material strength was negligible compared to gravity. This assumption can no longer be made when the bodies are much smaller. To this end, we have developed a 3D SPH code that includes a strength model to which we have added a von Mises yielding relation for stresses beyond the Hugoniot Elastic Limit. At the lower stresses associated with brittle failure, we use a rate-dependent strength based on the nucleation of incipient flaws whose number density is given by a Weibull distribution. Following Grady and Kipp and Melosh et al., we introduce a state variable D ('damage'), 0 less than D less than 1, which expresses the local reduction in strength due to crack growth under tensile loading. Unfortunately for the hydrodynamics, Grady and Kipp's model predicts which fragments are the most probable ones and not the ones that are really formed. This means, for example, that if a given laboratory experiment is modeled, the fragment distribution obtained from the Grady-Kipp theory would be equivalent to a ensemble average over many realizations of the experiment. On the other hand, the hydrodynamics itself is explicit and evolves not an ensemble average but very specific fragments. Hence, there is a clear incompatibility with the deterministic nature of the hydrodynamics equations and the statistical approach of the Grady-Kipp dynamical fracture model. We remedy these shortcomings by making the incipient flaw distribution explicit, i.e., particles carry activation strains which are distributed at random with a probability of occurrence given by the Weibull distribution. If the local principal axis strain exceeds this limit, damage starts to grow. By growing explicit cracks together with statistical cracks (damage) at the sub-particle scale, we ensure that material strength and fragmentation is independent of model resolution. We tested our scheme by simulating laboratory impact experiments on basalt spheres.

Benz, W.↗
NASA NTRS2002

Incorporation of Mean Stress Effects into the Micromechanical Analysis of the High Strain Rate Response of Polymer Matrix Composites

The results presented here are part of an ongoing research program, to develop strain rate dependent deformation and failure models for the analysis of polymer matrix composites subject to high strain rate impact loads. A micromechanics approach is employed in this work, in which state variable constitutive equations originally developed for metals have been modified to model the deformation of the polymer matrix, and a strength of materials based micromechanics method is used to predict the effective response of the composite. In the analysis of the inelastic deformation of the polymer matrix, the definitions of the effective stress and effective inelastic strain have been modified in order to account for the effect of hydrostatic stresses, which are significant in polymers. Two representative polymers, a toughened epoxy and a brittle epoxy, are characterized through the use of data from tensile and shear tests across a variety of strain rates. Results computed by using the developed constitutive equations correlate well with data generated via experiments. The procedure used to incorporate the constitutive equations within a micromechanics method is presented, and sample calculations of the deformation response of a composite for various fiber orientations and strain rates are discussed.

Goldberg, Robert K.↗
NASA NTRS2005

Fully Associative, Nonisothermal, Potential-Based Unified Viscoplastic Model for Titanium-Based Matrices

A number of titanium matrix composite (TMC) systems are currently being investigated for high-temperature air frame and propulsion system applications. As a result, numerous computational methodologies for predicting both deformation and life for this class of materials are under development. An integral part of these methodologies is an accurate and computationally efficient constitutive model for the metallic matrix constituent. Furthermore, because these systems are designed to operate at elevated temperatures, the required constitutive models must account for both time-dependent and time-independent deformations. To accomplish this, the NASA Lewis Research Center is employing a recently developed, complete, potential-based framework. This framework, which utilizes internal state variables, was put forth for the derivation of reversible and irreversible constitutive equations. The framework, and consequently the resulting constitutive model, is termed complete because the existence of the total (integrated) form of the Gibbs complementary free energy and complementary dissipation potentials are assumed a priori. The specific forms selected here for both the Gibbs and complementary dissipation potentials result in a fully associative, multiaxial, nonisothermal, unified viscoplastic model with nonlinear kinematic hardening. This model constitutes one of many models in the Generalized Viscoplasticity with Potential Structure (GVIPS) class of inelastic constitutive equations.

Source record↗
NASA NTRS2023

Investigating the Response of Land-Atmosphere Interactions and Feedbacks to Spatial Representation of Irrigation in a Coupled Modeling Framework

The transport of water, heat, and momentum from the surface to the atmosphere is dependent in part on the 10 characteristics of the land surface. Together with the model physics, parameterization schemes, and parameters employed, land datasets determine the spatial variability in land surface states (i.e., soil moisture and temperature) and fluxes. Despite the importance of these datasets, they are often chosen out of convenience or regional limitations without due assessment of their impacts on model results. Irrigation is an anthropogenic form of land heterogeneity that has been shown to alter the land surface energy balance, ambient weather, and local circulations. As such, irrigation schemes are becoming more 15 prevalent in weather and climate models with rapid developments in dataset availability and parameterization scheme complexity. Thus, to address pragmatic issues related to modeling irrigation, this study uses a high-resolution, regional coupled modeling system to investigate the impacts of irrigation dataset selection on land-atmosphere (L-A) coupling using a case study from the Great Plains Irrigation Experiment (GRAINEX) field campaign. The simulations are assessed in the context of irrigated versus non-irrigated regions, subregions across the irrigation gradient, and sub-grid scale process 20 representation in coarser scale models. The results show that L-A coupling is sensitive to the choice of irrigation dataset and resolution and that the irrigation impact on surface fluxes and near surface meteorology can be dominant, conditioned on the details of the irrigation map (i.e., boundaries, heterogeneity, etc), or minimal. A consistent finding across several analyses was that even a low percentage of irrigation fraction (i.e., 4-16%) can have significant local and downstream atmospheric impacts (e.g., lower PBL height), suggesting that representation of boundaries and heterogeneous areas within irrigated 25 regions is particularly important for the modeling of irrigation impacts on the atmosphere in this model. When viewing the simulations presented here as a proxy for ‘ideal’ tiling in a Earth System Model scale gridbox, the results show that some ‘tiles’ will reach critical nonlinear moisture and planetary boundary layer (PBL) thresholds that could be important for clouds and convection, implying that heterogeneity resulting from irrigation should be taken into consideration in new sub-grid land-atmosphere exchange parameterizations.

Patricia Lawston-Parker↗
NASA NTRS2021

A variable magnetic disc wind in the black hole X-ray binary GRS 1915+105?

Context. GRS 1915+105 being one of the brightest transient black hole binaries (BHBs) in the X-rays offers a unique testbed for the study of the connection between accretion and ejection mechanisms in BHBs. In particular, this source can be used to study the accretion disc wind and its dependence on the state changes in BHBs. Aims. Our aim is to investigate the origin and geometry of the accretion disc wind in GRS 1915+105. This study will provide a basis for planning future observations with the X-ray Imaging Spectroscopy Mission (XRISM), and may also provide important parameters for estimating the polarimetric signal with the upcoming Imaging X-ray Polarimetry Explorer (IXPE). Methods. We analysed the spectra of GRS 1915+105 in the soft ϕ and hard χ classes using the high-resolution spectroscopy offered by Chandra HETGS. In the soft state, we find a series of wind absorption lines that follow a non-linear dependence of velocity width, velocity shift, and equivalent width with respect to ionisation, indicating a multiple component or stratified outflow. In the hard state we find only a faint Fe XXVI absorption line. We model the absorption lines in both the states using a dedicated magneto-hydrodynamic (MHD) wind model to investigate a magnetic origin of the wind and to probe the cause of variability in the observed line flux between the two states. Conclusions. The MHD disc wind model provides a good fit for both states, indicating the possibility of a magnetic origin of the wind. The multiple ionisation components of the wind are well characterised as a stratification of the same magnetic outflow. We find that the observed variability in the line flux between soft and hard states cannot be explained by photo-ionisation alone but is most likely due to a large (three orders of magnitude) increase in the wind density. We find the mass outflow rate of the wind to be comparable to the accretion rate, suggesting an intimate link between accretion and ejection processes that lead to state changes in BHBs.

A. Ratheesh↗
NASA NTRS1973

Application of digital computer APU modeling techniques to control system design.

Study of the required controls for a H2-O2 auxiliary power unit (APU) technology program for the Space Shuttle. A steady-state system digital computer program was prepared and used to optimize initial system design. Analytical models of each system component were included. The program was used to solve a nineteen-dimensional problem, and then time-dependent differential equations were added to the computer program to simulate transient APU system and control. Some system parameters were considered quasi-steady-state, and others were treated as differential variables. The dynamic control analysis proceeded from initial ideal control modeling (which considered one control function and assumed the others to be ideal), stepwise through the system (adding control functions), until all of the control functions and their interactions were considered. In this way, the adequacy of the final control design over the required wide range of APU operating conditions was established.

Bailey, D. A.↗
NASA NTRS1994

Applications of automatic differentiation in computational fluid dynamics

Automatic differentiation (AD) is a powerful computational method that provides for computing exact sensitivity derivatives (SD) from existing computer programs for multidisciplinary design optimization (MDO) or in sensitivity analysis. A pre-compiler AD tool for FORTRAN programs called ADIFOR has been developed. The ADIFOR tool has been easily and quickly applied by NASA Langley researchers to assess the feasibility and computational impact of AD in MDO with several different FORTRAN programs. These include a state-of-the-art three dimensional multigrid Navier-Stokes flow solver for wings or aircraft configurations in transonic turbulent flow. With ADIFOR the user specifies sets of independent and dependent variables with an existing computer code. ADIFOR then traces the dependency path throughout the code, applies the chain rule to formulate derivative expressions, and generates new code to compute the required SD matrix. The resulting codes have been verified to compute exact non-geometric and geometric SD for a variety of cases. in less time than is required to compute the SD matrix using centered divided differences.

Green, Lawrence L.↗
NASA NTRS2004

Modeling the Nonlinear, Strain Rate Dependent Deformation of Woven Ceramic Matrix Composites With Hydrostatic Stress Effects Included

An analysis method based on a deformation (as opposed to damage) approach has been developed to model the strain rate dependent, nonlinear deformation of woven ceramic matrix composites with a plain weave fiber architecture. In the developed model, the differences in the tension and compression response have also been considered. State variable based viscoplastic equations originally developed for metals have been modified to analyze the ceramic matrix composites. To account for the tension/compression asymmetry in the material, the effective stress and effective inelastic strain definitions have been modified. The equations have also been modified to account for the fact that in an orthotropic composite the in-plane shear stiffness is independent of the stiffness in the normal directions. The developed equations have been implemented into a commercially available transient dynamic finite element code, LS-DYNA, through the use of user defined subroutines (UMATs). The tensile, compressive, and shear deformation of a representative plain weave woven ceramic matrix composite are computed and compared to experimental results. The computed values correlate well to the experimental data, demonstrating the ability of the model to accurately compute the deformation response of woven ceramic matrix composites.

Goldberg, Robert K.↗