Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “functional modeling”

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,153 records · Page 64

PolyFit: A C++ code for Polynomial Curve Fit with Calculation of Error Bars

In radiobiology, many dose-response results are modeled using the so-called linear-quadratic(LQ)model, which means that results are modeled as a function of dose Das 𝑅(𝐷)=𝛽0+𝛽1𝐷+𝛽2𝐷2. The coefficients 𝛽0, 𝛽1and 𝛽2are obtained from fitting a series of data points (𝑥𝑖,𝑦𝑖), which is usually done using a least-square method. The LQ and more generally the polynomial fit capability is implemented in many software that analyzes data. However, it is often convenient to do the fitting programmatically, especially when a large number of datasets should be analyzed. Furthermore, depending on the software used, some features may not be implemented. In this mini-review, I discuss the basis of polynomial fitting, including the calculation of errors on the coefficients and results, use of weighting and fixing the intercept value (the coefficient 𝛽0). A simple C++ code to perform the polynomial curve fitting is also provided. This code should be useful not only in radiobiology but in other fields of science as well.

Ianik Plante↗

Accuracy Assessment of Global Internal-Tide Models Using Satellite Altimetry Models Using Satellite Altimetry

Altimeter measurements are corrected for several geophysical parameters in order to access ocean signals of interest, like mesoscale or sub-mesoscale variability. The ocean tide is one of the most critical corrections due to the amplitude of the tidal elevations and to the aliasing phenomena of high-frequency signals into the lower-frequency band, but the internal-tide signatures at the ocean surface are not yet corrected globally. Internal tides can have a signature of several centimeters at the surface with wavelengths of about 50–250 km for the first mode and even smaller scales for higher-order modes. The goals of the upcoming Surface Water Ocean Topography (SWOT) mission and other high-resolution ocean measurements make the correction of these small-scale signals a challenge, as the correction of all tidal variability becomes mandatory to access accurate measurements of other oceanic signals. In this context, several scientific teams are working on the development of new internal-tide models, taking advantage of the very long altimeter time series now available, which represent an unprecedented and valuable global ocean database. The internal-tide models presented here focus on the coherent internal-tide signal and they are of three types: empirical models based upon analysis of existing altimeter missions, an assimilative model and a three-dimensional hydrodynamic model. A detailed comparison and validation of these internal-tide models is proposed using existing satellite altimeter databases. The analysis focuses on the four main tidal constituents: M2, K1, O1 and S2. The validation process is based on a statistical analysis of multi-mission altimetry including Jason-2 and Cryosphere Satellite-2 data. The results show a significant altimeter variance reduction when using internal-tide corrections in all ocean regions where internal tides are generating or propagating. A complementary spectral analysis also gives some estimation of the performance of each model as a function of wavelength and some insight into the residual non-stationary part of internal tides in the different regions of interest. This work led to the implementation of a new internal-tide correction (ZARON'one) in the next geophysical data records version-F (GDR-F) standards.

Loren Carrere↗

Initial Findings in Full-Shell Active Optics

The realization of active full-shell x-ray optics will provide ground-breaking capability for future missions. Future X-ray missions, such as Lynx, require high angular-resolution, large effective area, and as wide a field of view(FOV) as possible. It is currently not possible to perform high-resolution imaging of both wide and narrow FOVs with a single X-ray telescope. In this paper, we discuss the use of actuators to switch the optical surface of full-shell x-ray optics between prescriptions optimized for narrow-field and wide-field viewing, as well as correct for low-spatial-frequency errors in the optics. Previously, a pathfinder was developed using finite-element modeling(FEM). Using a combination of the model’s influence functions, the capacity of active full-shell optics was shown to switch between prescriptions with 92% accuracy. Another benefit of the analytical pathfinder was its ability to correct low-spatial frequency errors, which comprise a large percentage of the MSFC full-shell x-ray optics root-mean-square (RMS) slope error. In order to verify the pathfinder, a surface-parallel actuator was bonded to a nickel-replicated coupon. After the actuator was bonded to the coupon, the slope of the surface was measured while the actuator was activated. The proof-of-concept demonstrator development, the influence function data, and the resulting implications on the analytical pathfinder are reported and discussed.

X-ray optics↗

Numerical Investigation of the Effect of Trailing Edge Thickness of Simulated CMC Blades on Loss Profiles

To understand such an aerodynamic characteristic of the CMC blades, we numerically simulated three different CMC blades with varying trailing edge thicknesses. We utilized both Large Eddy simulation (LES) with the dynamic subgrid closure as well a recently implemented Reynolds-averaged Navier-Stokes (RANS) coupled with an intermittency function-based transition model. We considered two sets of experimental test points with different Reynolds numbers at a given high-free stream turbulence range (Tu=~10-13%). The LES was previously validated by comparing with experimental data. We validated the RANS transition model against the ERCOFTAC experimental zero-pressure gradient smooth flat plate boundary layer at low- and high freestream turbulence offered by the experimental data of the Stripf’s for a turbine blade. We confirmed that the predicted skin friction and the heat flux agree well with the data. Finally, we computed the predicted pressure loading profiles, total loss distributions, and integrated loss of the three different CMC blades and made comparisons against the experimental data acquired at NASA Glenn Transonic Turbine Blade Cascade Rig. Our paper will present the aforementioned results and discuss the ability of LES and RANS model of our choice to predict salient aerodynamic characteristics of the blades with the various trailing edge thickness effects.

LES↗

Numerical Investigation of the Effect of Trailing Edge Thickness of Simulated CMC Blades on Loss Profiles

To understand such an aerodynamic characteristic of the CMC blades, we numerically simulated three different CMC blades with varying trailing edge thicknesses. We utilized both Large Eddy simulation (LES) with the dynamic subgrid closure as well a recently implemented Reynolds-averaged Navier-Stokes (RANS) coupled with an intermittency function-based transition model. We considered two sets of experimental test points with different Reynolds numbers at a given high-free stream turbulence range (Tu=~10-13%). The LES was previously validated by comparing with experimental data. We validated the RANS transition model against the ERCOFTAC experimental zero-pressure gradient smooth flat plate boundary layer at low- and high freestream turbulence offered by the experimental data of the Stripf’s for a turbine blade. We confirmed that the predicted skin friction and the heat flux agree well with the data. Finally, we computed the predicted pressure loading profiles, total loss distributions, and integrated loss of the three different CMC blades and made comparisons against the experimental data acquired at NASA Glenn Transonic Turbine Blade Cascade Rig. Our paper will present the aforementioned results and discuss the ability of LES and RANS model of our choice to predict salient aerodynamic characteristics of the blades with the various trailing edge thickness effects.

LES↗

A Hybrid CFD/Engineering Model for Predicting Plume Induced Erosion and Cratering

The Descent Interpolated Gas-Granular Erosion Model (DIGGEM) is a hybrid computational fluid dynamics (CFD) / engineering model calibrated with flight data from Apollo. It is used for the prediction of plume induced viscous erosion and cratering during the descent of landing vehicles. The model assumes a functional relationship between local shear stress and eroded mass flux. It is used along with a suite of other tools to predict complex plume induced environments for Human Landing System (HLS) and Commercial Lunar Payload Services (CLPSS) vehicles.

PSI↗

Nitrogen in Ancient Mud: A Biosignature?

Nitrogen is an essential nutrient for all life on Earth and possibly elsewhere. Burial of nitrogen bound to organic matter constitutes the major flux of nitrogen into sediments today, which has led to the inference that nitrogen enrichments in sedimentary rocks may be a biosignature. However, abiotic processes such as lightning or volcanism can fix atmospheric N2 and contribute to sedimentary nitrogen burial in the absence of life. It is therefore uncertain whether observed nitrogen enrichments of up to 430 ppm in Paleoarchean metasedimentary biotite grains are indeed biogenic. This study seeks to address that problem with a numerical model. The NH4(+) concentration of an abiotic ocean is modeled as a function of source fluxes, pH-dependent NH3 volatilization, and equilibrated adsorption of NH4(+) onto clay particles. The results suggest that the observed nitrogen concentrations in Paleoarchean biotite can only be reconciled with purely abiotic processes if the ocean was more acidic (pH <6) and/or if the source fluxes from lightning and volcanism were at least an order of magnitude higher (≥10(exp 12) mol/yr) than previously thought. The bulk of the nitrogen is thus most likely of biological origin. While this does not necessitate a particular metabolism such as biological N2 fixation, the data provide evidence of nitrogen utilization back to 3.8 Gyr. Nitrogen abundances could thus provide useful information in extraterrestrial missions.

Early Earth↗

The Effect of Trailing Edge Grid Resolution in Large-Eddy Simulations

Ceramic Matrix Composite (CMC) with environmental barrier coating can tolerate significantly higher temperatures compared with the conventional metal alloy blades and thus can help achieve higher thermal efficiency in gas turbine engines. Due to the complex fabrication process of CMC blades, larger blade trailing edge thicknesses and larger leading-edge radii of curvature are expected. These features significantly alter aerodynamic performance especially near the trailing edge and thus the loss profile. In previous studies, the authors used the large eddy simulation (LES) as well as the Reynolds Averaged Navier-Stokes (RANS) coupled with an intermittency function-based transition model to simulate an aerodynamics performance of three different CMC blades and compared the predictions against the data acquired at NASA Glenn Transonic Turbine Blade Cascade Rig. Although we were able to successfully predict different characteristics of aerodynamic performance among three different blades (e.g., pressure loading and integrated losses), LES consistently predicted a small recirculation bubble at the end of the suction side. This was not inconsistent with the physics of the flow, but the pressure loading, locally, near the trailing edge, and as a result, the measured wake profiles did not match the experimental data. To delve further into these issues, a new LES grid with a much finer mesh near the trailing edge was generated where the spanwise extent of the grid was increased by 50%. An LES simulation for the CMC blade with the largest trailing edge thicknesses at the high-freestream turbulence (Tu=13%) condition showed that using such a mesh is necessary for LES to accurately capture the transition phenomena on the suction side. The resulting total pressure coefficient profile and the shift of the wake profile were greatly improved.

CMC↗

Assurance Equations: A Cost and Criticality Model for Optimizing Quality Assurance Surveillance

The cost of quality vs cost of failure correction has been a long-running topic of discussion within the Aerospace community. It leads directly to concepts of “risk tolerance”, and risk-based decision-making. It would be valuable if there was a way to compute the optimal investment in customer-executed quality assurance activities using defect significance with respect to performance objectives, the activities’ defect detection effectiveness, and the cost-penalty for late discovery of impactful defects. This optimization is particularly of interest to projects whose budget constraints significantly limit their risk management options.The cost to fix defects (i.e., failure correction) escalates as the project matures. There have been studies attempting to determine the relative cost of fixing defects discovered during various phases of a project life cycle with important implications, all of which suggest growth factors are large. The commonly referred to 1:10:100 rule represents a cost multiplier for repair/rework across the Design to Fab to Test hardware development phases. Cost premiums for QA activities also accumulate when they are treated as mandatory (due to schedule drag) or are performed later than their assigned phase.This paper describes the modeling of development phase -dependencies in the conduct of typical customer-executed quality assurance activities. Our initial modeling encompasses:• Distinct phases of the production lifecycle• Multiple kinds of Defects, each with some a-priori likelihood of being present• Each defect’s impact on performance Objectives for a type of hardware• The cost and efficacy of assurance techniques at detecting such Defects• The costs of fixing those Defects detected in a given phase of the production lifecycleThe model captures assurance activities’ abilities to Detect defects. Upon detection it is assumed that the Defect is immediately fixed. Defects that “escape” detection by some activity may thereafter be detected by a later activity, but by then the cost of fixing the Defect may have escalated. Defects are related to the performance Objectives they would detract from, were those Defects to remain present in the operating system.We have constructed and are exploring, a model that relates the importance of hardware system elements to mission objectives, the impact of types of Defects on those hardware types, the cost of customer-executed assurance activities (i.e., supplier controls) and their effectiveness towards reducing an impactful quality escape, and the cost of Defect correction across production phase. We describe the approach taken to select the key model aspects, why they are relevant to our NASA mission, and our efforts to populate it with relevant and contemporary data. We use a notional example to illustrate model design and function.

Plante, Jeannette↗

The Effect of Trailing Edge Grid Resolution in Large-Eddy Simulations

Ceramic Matrix Composite (CMC) with environmental barrier coating can tolerate significantly higher temperatures compared with the conventional metal alloy blades and thus can help achieve higher thermal efficiency in gas turbine engines. Due to the complex fabrication process of CMC blades, larger blade trailing edge thicknesses and larger leading-edge radii of curvature are expected. These features significantly alter aerodynamic performance especially near the trailing edge and thus the loss profile. In previous studies, the authors used the large eddy simulation (LES) as well as the Reynolds Averaged Navier-Stokes (RANS) coupled with an intermittency function-based transition model to simulate an aerodynamics performance of three different CMC blades and compared the predictions against the data acquired at NASA Glenn Transonic Turbine Blade Cascade Rig. Although we were able to successfully predict different characteristics of aerodynamic performance among three different blades (e.g., pressure loading and integrated losses), LES consistently predicted a small recirculation bubble at the end of the suction side. This was not inconsistent with the physics of the flow, but the pressure loading, locally, near the trailing edge, and as a result, the measured wake profiles did not match the experimental data. To delve further into these issues, a new LES grid with a much finer mesh near the trailing edge was generated where the spanwise extent of the grid was increased by 50%. An LES simulation for the CMC blade with the largest trailing edge thicknesses at the high-freestream turbulence (Tu=13%) condition showed that using such a mesh is necessary for LES to accurately capture the transition phenomena on the suction side. The resulting total pressure coefficient profile and the shift of the wake profile were greatly improved.

LES↗

Applicability of Fly Neighborly Noise Recommendations to UAM Quadrotors Undergoing Steady Maneuvers

This initial study examines whether the Helicopter Association International (HAI) Fly Neighborly operational recommendations that are based on single main rotor/tail rotor configurations will hold for non-conventional UAM rotorcraft with multiple rotors. The 6-occupant quadrotor concept vehicle designed under the NASA Revolutionary Vertical Lift Technology (RVLT) Project is studied. The tip speed is 550 ft/sec, with three blades per rotor (“550/3”). Predictions are made for three steady maneuvers: level turns, descending turns, and climbing turns. The RVLT Toolchain is exercised using CAMRAD II, pyaaron/AARON/ANOPP2 and a beta version of AMAT (ANOPP2 Mission Analysis Tool). AMAT provides functionality to acoustically model the curved flight paths associated with maneuvers. In addition to quadrotor trim and performance, this study includes analysis of azimuthal variations of the vertical blade loading and its derivative in the form of contour plots on a rotor plane and line plots at one radial location (r/R=0.765). Quadrotor noise trends are analyzed using maximum Overall Sound Pressure Level (OASPL) and Effective Perceived Noise Level (EPNL). The Fly Neighborly guideline that addresses descents (“Level turns are quieter than descending turns”) is predicted to hold for the RVLT Quadrotor as well.

UAM↗

Implementation of a Goal-Based Systems Engineering Process Using the Systems Modeling Language (SysML)

This paper describes the core framework used to implement a Goal-Function Tree (GFT) based systems engineering process using the Systems Modeling Language. It defines a set of principles built upon by the theoretical approach described in the InfoTech 2013 ISHM paper titled "Goal-Function Tree Modeling for Systems Engineering and Fault Management" presented by Dr. Stephen B. Johnson. Using the SysML language, the principles in this paper describe the expansion of the SysML language as a baseline in order to: hierarchically describe a system, describe that system functionally within success space, and allocate detection mechanisms to success functions for system protection.

Breckenridge, Jonathan T.↗

Discrete event simulation tool for analysis of qualitative models of continuous processing systems

An artificial intelligence design and qualitative modeling tool is disclosed for creating computer models and simulating continuous activities, functions, and/or behavior using developed discrete event techniques. Conveniently, the tool is organized in four modules: library design module, model construction module, simulation module, and experimentation and analysis. The library design module supports the building of library knowledge including component classes and elements pertinent to a particular domain of continuous activities, functions, and behavior being modeled. The continuous behavior is defined discretely with respect to invocation statements, effect statements, and time delays. The functionality of the components is defined in terms of variable cluster instances, independent processes, and modes, further defined in terms of mode transition processes and mode dependent processes. Model construction utilizes the hierarchy of libraries and connects them with appropriate relations. The simulation executes a specialized initialization routine and executes events in a manner that includes selective inherency of characteristics through a time and event schema until the event queue in the simulator is emptied. The experimentation and analysis module supports analysis through the generation of appropriate log files and graphics developments and includes the ability of log file comparisons.

Malin, Jane T.↗

Formulation of consumables management models: Mission planning processor payload interface definition

Consumables models required for the mission planning and scheduling function are formulated. The relation of the models to prelaunch, onboard, ground support, and postmission functions for the space transportation systems is established. Analytical models consisting of an orbiter planning processor with consumables data base is developed. A method of recognizing potential constraint violations in both the planning and flight operations functions, and a flight data file storage/retrieval of information over an extended period which interfaces with a flight operations processor for monitoring of the actual flights is presented.

Torian, J. G.↗

A Residuals Approach to Filtering, Smoothing and Identification for Static Distributed Systems

An approach for state estimation and identification of spatially distributed parameters embedded in static distributed (elliptic) system models is advanced. The method of maximum likelihood is used to find parameter values that maximize a likelihood functional for the system model, or equivalently, that minimize the negative logarithm of this functional. To find the minimum, a Newton-Raphson search is conducted that from an initial estimate generates a convergent sequence of parameter estimates. For simplicity, a Gauss-Markov approach is used to approximate the Hessian in terms of products of first derivatives. The gradient and approximate Hessian are computed by first arranging the negative log likelihood functional into a form based on the square root factorization of the predicted covariance of the measurement process. The resulting data processing approach, referred to here by the new term of predicted data covariance square root filtering, makes the gradient and approximate Hessian calculations very simple. A closely related set of state estimates is also produced by the maximum likelihood method: smoothed estimates that are optimal in a conditional mean sense and filtered estimates that emerge from the predicted data covariance square root filter.

Rodriguez, G.↗

Development of a Hybrid RANS/LES Method for Turbulent Mixing Layers

Significant research has been underway for several years in NASA Glenn Research Center's nozzle branch to develop advanced computational methods for simulating turbulent flows in exhaust nozzles. The primary efforts of this research have concentrated on improving our ability to calculate the turbulent mixing layers that dominate flows both in the exhaust systems of modern-day aircraft and in those of hypersonic vehicles under development. As part of these efforts, a hybrid numerical method was recently developed to simulate such turbulent mixing layers. The method developed here is intended for configurations in which a dominant structural feature provides an unsteady mechanism to drive the turbulent development in the mixing layer. Interest in Large Eddy Simulation (LES) methods have increased in recent years, but applying an LES method to calculate the wide range of turbulent scales from small eddies in the wall-bounded regions to large eddies in the mixing region is not yet possible with current computers. As a result, the hybrid method developed here uses a Reynolds-averaged Navier-Stokes (RANS) procedure to calculate wall-bounded regions entering a mixing section and uses a LES procedure to calculate the mixing-dominated regions. A numerical technique was developed to enable the use of the hybrid RANS-LES method on stretched, non-Cartesian grids. With this technique, closure for the RANS equations is obtained by using the Cebeci-Smith algebraic turbulence model in conjunction with the wall-function approach of Ota and Goldberg. The LES equations are closed using the Smagorinsky subgrid scale model. Although the function of the Cebeci-Smith model to replace all of the turbulent stresses is quite different from that of the Smagorinsky subgrid model, which only replaces the small subgrid turbulent stresses, both are eddy viscosity models and both are derived at least in part from mixing-length theory. The similar formulation of these two models enables the RANS and LES equations to be solved with a single solution scheme and computational grid. The hybrid RANS-LES method has been applied to a benchmark compressible mixing layer experiment in which two isolated supersonic streams, separated by a splitter plate, provide the flows to a constant-area mixing section. Although the configuration is largely two dimensional in nature, three-dimensional calculations were found to be necessary to enable disturbances to develop in three spatial directions and to transition to turbulence. The flow in the initial part of the mixing section consists of a periodic vortex shedding downstream of the splitter plate trailing edge. This organized vortex shedding then rapidly transitions to a turbulent structure, which is very similar to the flow development observed in the experiments. Although the qualitative nature of the large-scale turbulent development in the entire mixing section is captured well by the LES part of the current hybrid method, further efforts are planned to directly calculate a greater portion of the turbulence spectrum and to limit the subgrid scale modeling to only the very small scales. This will be accomplished by the use of higher accuracy solution schemes and more powerful computers, measured both in speed and memory capabilities.

Georgiadis, Nicholas J.↗

Deep Space Radiation Affects Neurovascular Functions in Human Organ-on-a-Chip Models

A major health risk for human deep space exploration is central nervous system (CNS) damage by galactic cosmic ray radiation. Simulated galactic cosmic rays or their components, especially the high- linear energy transfer (LET) particles such as 56 Fe ions, cause CNS damage, neuroinflammation and cognitive dysfunction in rodent models, but their effects on human CNS remain to be investigated. CNS damage from any insult, including ionizing radiation, is partially mediated by the blood-brain barrier (BBB), which regulates the interactions between CNS and the rest of the body. The main cellular regulators of BBB permeability are astrocytes, which also modulate neuronal health and neuroinflammation. However, there have been few studies on BBB and astrocyte functions in regulating CNS responses, especially in human tissue/organ analogs. Therefore, we utilized a high-throughput human 3D organ-on-a-chip system, seeded with induced pluripotent stem cell-derived endothelial cells, astrocytes and neurons, to study human neurovascular responses to simulated deep space radiation. We investigated BBB permeability, oxidative stress, cellular and tissue damage, and secreted factors over the time period of 24 hours-1 week after irradiation with 0.25-0.5 Gy 5-ion simplified simulated galactic cosmic rays and 0.3-0.8 Gy high-LET 600MeV/n 56 Fe particles, and compared the outcomes to low-LET irradiation with 0.1-1 Gy doses of X-rays and gamma rays. Both high and low-LET radiation increased neurovascular permeability, caused oxidative stress, damaged endothelial cells and tight junctions, and altered expression of inflammatory cytokines. Ionizing radiation- induced neurovascular permeability and oxidative stress peaked at 3 days after irradiation and were further exacerbated by the presence of astrocytes. Furthermore, in response to particle irradiation, astrocytes stimulated interleukin-1 signaling by inhibiting the expression of interleukin-1 receptor antagonist. Thus, we also evaluated interleukin-1 receptor antagonist as a potential countermeasure against particle radiation. Ultimately, our results may help develop countermeasures to mitigate human CNS damage in deep space exploration.

Sonali D Verma↗