SEARCH · Search NASA

Results for “nonlinear model reduction”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 253 records · Page 14

All sources NASA NTRS DOE OSTI

NASA NTRS1997

Controls/CFD Interdisciplinary Research Software Generates Low-Order Linear Models for Control Design From Steady-State CFD Results

The NASA Lewis Research Center is developing analytical methods and software tools to create a bridge between the controls and computational fluid dynamics (CFD) disciplines. Traditionally, control design engineers have used coarse nonlinear simulations to generate information for the design of new propulsion system controls. However, such traditional methods are not adequate for modeling the propulsion systems of complex, high-speed vehicles like the High Speed Civil Transport. To properly model the relevant flow physics of high-speed propulsion systems, one must use simulations based on CFD methods. Such CFD simulations have become useful tools for engineers that are designing propulsion system components. The analysis techniques and software being developed as part of this effort are an attempt to evolve CFD into a useful tool for control design as well. One major aspect of this research is the generation of linear models from steady-state CFD results. CFD simulations, often used during the design of high-speed inlets, yield high resolution operating point data. Under a NASA grant, the University of Akron has developed analytical techniques and software tools that use these data to generate linear models for control design. The resulting linear models have the same number of states as the original CFD simulation, so they are still very large and computationally cumbersome. Model reduction techniques have been successfully applied to reduce these large linear models by several orders of magnitude without significantly changing the dynamic response. The result is an accurate, easy to use, low-order linear model that takes less time to generate than those generated by traditional means. The development of methods for generating low-order linear models from steady-state CFD is most complete at the one-dimensional level, where software is available to generate models with different kinds of input and output variables. One-dimensional methods have been extended somewhat so that linear models can also be generated from two- and three-dimensional steady-state results. Standard techniques are adequate for reducing the order of one-dimensional CFD-based linear models. However, reduction of linear models based on two- and three-dimensional CFD results is complicated by very sparse, ill-conditioned matrices. Some novel approaches are being investigated to solve this problem.

Melcher, Kevin J.↗

NASA NTRS2002

In-Situ Measurements of HCN and CH3CN In the Pacific Troposphere: Sources, Sinks, and Comparisons with Satellite Observations

A new capillary gas chromatographic method using a Reduction Gas Detector was developed to measure HCN and CH3CN in the remote troposphere. This instrumental configuration was deployed for the very first time in the Trace-P field mission performed during the spring of 2001. The NASA DC-8 aircraft afforded an opportunity to measure HCN and CH3CN in polluted and pristine environments over the Pacific to a maximum altitude of 12 km. These are some of the first in situ measurements of the distribution of HCN and CH3CN over the Pacific. Large background concentrations of both nitriles were found to be present and significant variability was observed. The abundance of HCN and CH3CN was strongly impacted by outflow of pollution from Asia. In general there appeared to be a direct but nonlinear relationship between the mixing ratios of HCN and CH3CN. The vertical structure of these chemicals shows direct evidence of the presence of a significant oceanic sink. These observations will be compared with the column content HCN data from satellites and other available measurements. A large body of data have been collected and are being analyzed, both statistically and with the help of models, to better understand the sources and sinks of these nitriles. These results will be presented.

Singh, Hanwant B.↗

NASA NTRS1993

Moist wind relationships

Equations describing the temporal and spatial behavior of the kinematic moisture and heat flux are introduced. In these nonlinear equations, the contribution by diabatic processes to the large-scale flux is composed of two parts. One part is associated with a Rayleigh damping term, while the other arises from temporal and spatial changes in the pressure gradient term. The influence of diabatic processes on the large-scale moisture fluxes depends greatly on the degree of balance between forcing and damping terms in the governing equations. The existence of a near balance requires a reduction in the large-scale horizontal geostrophic wind speed. From a scale analysis of the moisture flux equations it is argued that reductions in the large-scale horizontal wind speed, observed within major cumulus cloud systems, help conserve large-scale moisture fluxes. The deviation of the wind from geostrophic conditions is easily estimated. This wind modification induces secondary vertical circulations that contribute to the convergence, creating or supporting long-lived mesoscale flows. In the tropics the wind modification has an antitriptic relationship. These diagnostic findings suggest possible modifications to the wind field in the application of cumulus parameterization, and may be important in diabatic initialization of numerical weather prediction models.

Raymond, William H.↗

NASA NTRS1999

A fluid--structure interaction finite element analysis of pulsatile blood flow through a compliant stenotic artery

A new model is used to analyze the fully coupled problem of pulsatile blood flow through a compliant, axisymmetric stenotic artery using the finite element method. The model uses large displacement and large strain theory for the solid, and the full Navier-Stokes equations for the fluid. The effect of increasing area reduction on fluid dynamic and structural stresses is presented. Results show that pressure drop, peak wall shear stress, and maximum principal stress in the lesion all increase dramatically as the area reduction in the stenosis is increased from 51 to 89 percent. Further reductions in stenosis cross-sectional area, however, produce relatively little additional change in these parameters due to a concomitant reduction in flow rate caused by the losses in the constriction. Inner wall hoop stretch amplitude just distal to the stenosis also increases with increasing stenosis severity, as downstream pressures are reduced to a physiological minimum. The contraction of the artery distal to the stenosis generates a significant compressive stress on the downstream shoulder of the lesion. Dynamic narrowing of the stenosis is also seen, further augmenting area constriction at times of peak flow. Pressure drop results are found to compare well to an experimentally based theoretical curve, despite the assumption of laminar flow.

Non-NASA Center↗

NASA NTRS1993

Space Shuttle Main Engine performance analysis

For a number of years, NASA has relied primarily upon periodically updated versions of Rocketdyne's power balance model (PBM) to provide space shuttle main engine (SSME) steady-state performance prediction. A recent computational study indicated that PBM predictions do not satisfy fundamental energy conservation principles. More recently, SSME test results provided by the Technology Test Bed (TTB) program have indicated significant discrepancies between PBM flow and temperature predictions and TTB observations. Results of these investigations have diminished confidence in the predictions provided by PBM, and motivated the development of new computational tools for supporting SSME performance analysis. A multivariate least squares regression algorithm was developed and implemented during this effort in order to efficiently characterize TTB data. This procedure, called the 'gains model,' was used to approximate the variation of SSME performance parameters such as flow rate, pressure, temperature, speed, and assorted hardware characteristics in terms of six assumed independent influences. These six influences were engine power level, mixture ratio, fuel inlet pressure and temperature, and oxidizer inlet pressure and temperature. A BFGS optimization algorithm provided the base procedure for determining regression coefficients for both linear and full quadratic approximations of parameter variation. Statistical information relative to data deviation from regression derived relations was also computed. A new strategy for integrating test data with theoretical performance prediction was also investigated. The current integration procedure employed by PBM treats test data as pristine and adjusts hardware characteristics in a heuristic manner to achieve engine balance. Within PBM, this integration procedure is called 'data reduction.' By contrast, the new data integration procedure, termed 'reconciliation,' uses mathematical optimization techniques, and requires both measurement and balance uncertainty estimates. The reconciler attempts to select operational parameters that minimize the difference between theoretical prediction and observation. Selected values are further constrained to fall within measurement uncertainty limits and to satisfy fundamental physical relations (mass conservation, energy conservation, pressure drop relations, etc.) within uncertainty estimates for all SSME subsystems. The parameter selection problem described above is a traditional nonlinear programming problem. The reconciler employs a mixed penalty method to determine optimum values of SSME operating parameters associated with this problem formulation.

Santi, L. Michael↗

NASA NTRS1990

Moist wind relationships

The equations describing the temporal and spatial behavior of the kinematic moisture and heat flux are described. In these nonlinear equations, the contribution by diabatic processes to the large-scale flux is composed of two parts. One part is associated with a Rayleigh damping term while the other arises from temporal and spatial changes in the pressure gradient term. It was found that the influence of diabatic processes on large-scale moisture fluxes depends greatly on the degree of balance between forcing and damping terms in the governing kinematic flux equations. The existence of a near balance requires a reduction in the large-scale horizontal geostrophic wind speed. Based on an examination of the moisture flux equations, it is argued that reductions in the large-scale horizontal wind speed observed within major cumulus cloud systems help conserve large-scale moisture fluxes. The deviation of the wind from geostrophic conditions is easily estimated for a near balanced state. This wind modification induces secondary vertical circulations that contribute to convergence, creating or supporting long-lived mesoscale flows. We believe this process to be a major supporter of the mesoscale circulations observed in severe storms and squall lines. In the tropics the wind modification has an antitriptic relationship. These diagnostic findings suggest possible modifications to the wind field in the application of a cumulus parameterization, and may be important in diabatic initialization of numerical weather prediction models.

Raymond, William H.↗

NASA NTRS1993

An investigation of fighter aircraft agility

This report attempts to unify in a single document the results of a series of studies on fighter aircraft agility funded by the NASA Ames Research Center, Dryden Flight Research Facility and conducted at the University of Kansas Flight Research Laboratory during the period January 1989 through December 1993. New metrics proposed by pilots and the research community to assess fighter aircraft agility are collected and analyzed. The report develops a framework for understanding the context into which the various proposed fighter agility metrics fit in terms of application and testing. Since new metrics continue to be proposed, this report does not claim to contain every proposed fighter agility metric. Flight test procedures, test constraints, and related criteria are developed. Instrumentation required to quantify agility via flight test is considered, as is the sensitivity of the candidate metrics to deviations from nominal pilot command inputs, which is studied in detail. Instead of supplying specific, detailed conclusions about the relevance or utility of one candidate metric versus another, the authors have attempted to provide sufficient data and analyses for readers to formulate their own conclusions. Readers are therefore ultimately responsible for judging exactly which metrics are 'best' for their particular needs. Additionally, it is not the intent of the authors to suggest combat tactics or other actual operational uses of the results and data in this report. This has been left up to the user community. Twenty of the candidate agility metrics were selected for evaluation with high fidelity, nonlinear, non real-time flight simulation computer programs of the F-5A Freedom Fighter, F-16A Fighting Falcon, F-18A Hornet, and X-29A. The information and data presented on the 20 candidate metrics which were evaluated will assist interested readers in conducting their own extensive investigations. The report provides a definition and analysis of each metric; details of how to test and measure the metric, including any special data reduction requirements; typical values for the metric obtained using one or more aircraft types; and a sensitivity analysis if applicable. The report is organized as follows. The first chapter in the report presents a historical review of air combat trends which demonstrate the need for agility metrics in assessing the combat performance of fighter aircraft in a modern, all-aspect missile environment. The second chapter presents a framework for classifying each candidate metric according to time scale (transient, functional, instantaneous), further subdivided by axis (pitch, lateral, axial). The report is then broadly divided into two parts, with the transient agility metrics (pitch lateral, axial) covered in chapters three, four, and five, and the functional agility metrics covered in chapter six. Conclusions, recommendations, and an extensive reference list and biography are also included. Five appendices contain a comprehensive list of the definitions of all the candidate metrics; a description of the aircraft models and flight simulation programs used for testing the metrics; several relations and concepts which are fundamental to the study of lateral agility; an in-depth analysis of the axial agility metrics; and a derivation of the relations for the instantaneous agility and their approximations.

Valasek, John↗

NASA NTRS2006

Diurnal Variability of Vertical Structure from a TRMM Passive Microwave "Virtual Radar" Retrieval

Robust description of the diurnal cycle from TRMM observations is complicated by the limitations of Low Earth Orbit (LEO) sampling; from a 'climatological' perspective, sufficient sampling must exist to control for both spatial and seasonal variability, before tackling an additional diurnal component (e.g., with 8 additional 3-hourly or 24 1-hourly bins). For documentation of vertical structure, the narrow sample swath of the TRMM Precipitation Radar limits the resolution of any of these components. A neural-network based 'virtual radar" retrieval has been trained and internally validated, using multifrequency / multipolarization passive microwave(TM1) brightness temperatures and textures parameters and lightning (LIS) observations, as inputs, and PR volumetric reflectivity as targets (outputs). By training the algorithms (essentially highly multivariate, nonlinear regressions) on a very large sample of high-quality co-located data from the center of the TRMM swath, 3D radar reflectivity and derived parameters (VIL, IWC, Echo Tops, etc.) can be retrieved across the entire TMI swath, good to 8-9% over the dynamic range of parameters. As a step in the retrieval (and as an output of the process), each TMI multifrequency pixel (at 85 GHz resolution) is classified into one of the 25 archetypal radar profile vertical structure "types", previously identified using cluster analysis. The dynamic range of retrieved vertical structure appears to have higher fidelity than the current (Version 6) experimental GPROF hydrometeor vertical structure retrievals. This is attributable to correct representation of the prior probabilities of vertical structure variability in the neural network training data, unlike the GPROF cloud-resolving model training dataset used in the V6 algorithms. The LIS lightning inputs are supplementary inputs, and a separate offline neural network has been trained to impute (predict) LIS lightning from passive-microwave-only data. The virtual radar retrieval is thus, in principle, extensible to Aqua/AMSR-E and NPOESS/CMIS passive microwave instruments. The virtual radar approach yields a threefold increase in effective sampling from the mission, albeit of lower-quality "retrieved" data, reducing the variance of local estimates by one third (or the standard deviation by-0.57). In this talk, the variance reduction is leveraged to more finely resolve global diurnal variability in both space and time (local hour).

Boccippio, Dennis J.↗

DOE OSTI2025

Investigating the Role of Accident Tolerant Cladding on Source Term Reduction for High-Burnup PWRs Using MELCOR

The use of accident tolerant fuel (ATF) cladding can increase coping times during and beyond design basis accidents. While such gains may be incremental, they provide a margin that can potentially be recovered to enable high-burnup (HBU) operation. Realizing such a margin requires demonstrating that the combination of HBU and ATF has not led to an overall increase in source term. This study investigates the influence of cladding technology (Zr-based, Cr-coated Zr, and FeCrAl) and fuel cycle length (18 and 24 months) on radiological dose at the boundary of the exclusion zone for a four-loop pressurized water reactor to investigate whether ATF claddings can provide such benefits. We analyze a recovered large break loss-of-coolant accident scenario to investigate the impact of transient timescale on the benefits of such coping time increases. The simulations have been performed using the MELCOR and MELCOR Accident Consequence Code System codes. For the cases analyzed, increased fuel cycle length did not necessarily increase radionuclide release and hydrogen generation, as these were found to be sensitive to the core power distribution. Similarly, off-site dose consequence is dominated by short-lived radionuclides that tend to saturate earlier in the burnup, so higher burnup operation did not necessarily increase the source term for the phenomena and transients analyzed here. Delays in recovery of the lowpressure safety injection system increase hydrogen production and radionuclide release, especially between 780 s and 1620 s, due to the nonlinear oxidation and core degradation behavior. Results show that Cr-coated Zr enhances safety by delaying heatup and gap release. Here, when uncertainty propagation on oxidation properties is considered, FeCrAl exhibits the lowest overall radionuclide release and off-site dose throughout the spectrum. However, while the considered “base model” performance is superior under delayed injection scenarios, upper-bound cases display hydrogen generation risk comparable to the Zr-based cladding.

Accident Tolerant Fuel↗

NASA NTRS2006

Nonlinear Spring Finite Elements for Predicting Mode I-Dominated Delamination Growth in Laminated Structure with Through-Thickness reinforcement

One particular concern of polymer matrix composite laminates is the relatively low resistance to delamination cracking, in particular when the dominant type of failure is mode I opening. One method proposed for alleviating this problem involves the insertion pultruded carbon pins through the laminate thickness. The pins, known as z-pins, are inserted into the prepreg laminate using an ultrasonic hammer prior to the curing process, resulting in a field of pins embedded normal to the laminate plane as illustrated in Figure. 1. Pin diameters range between 0.28-mm to 0.5-mm and standard areal densities range from 0.5% to 4%. The z-pins are provided by the manufacturer, Aztex(Registered TradeMark) , in a low-density foam preform, which acts to stabilize orientation of the pins during the insertion process [1-3]. Typical pin materials include boron and carbon fibers embedded in a polymer matrix. A number of methods have been developed for predicting delamination growth in laminates reinforced with z-pins. During a study on the effect of z-pin reinforcement on mode I delamination resistance, finite element analyses of z-pin reinforced double cantilever beam (DCB) specimens were performed by Cartie and Partridge [4]. The z-pin bridging stresses were modeled by applying equivalent forces at the pin locations. Single z-pin pull-out tests were performed to characterize the traction law of the pins under mode I loading conditions. Analytical solutions for delamination growth in z-pin reinforced DCB specimens were independently derived by Robinson and Das [5] and Ratcliffe and O'Brien [6]. In the former case, pin bridging stresses were modeled using a distributed load and in the latter example the bridging stresses were discretely modeled by way of grounded springs. Additionally, Robinson and Das developed a data reduction strategy for calculating mode I fracture toughness, G(sub Ic), from a z-pin reinforced DCB specimen test [5]. In both cases a traction law similar to that adopted by Cartie and Partridge was used to represent z-pin failure under mode I loading conditions. In the current work spring elements available in most commercial finite element codes were used to model z-pins. The traction law used in previous analyses [4-6] was employed to represent z-pin damage. This method is intended for and is limited to simulating z-pins in composite laminate structure containing mode I-dominated delamination cracking. The current technique differs from previous analyses in that spring finite elements (available in commercial codes) are employed for simulating zpins, reducing the complexity of the analysis construction process. Furthermore, the analysis method can be applied to general structure that experiences mode I-dominated delamination cracking, in contrast to existing analytical solutions that are only applicable to coupon DCB specimens.

Ratcliffe, James G.↗

NASA NTRS1999

Hyperspectral Remote Sensing and Ecological Modeling Research and Education at Mid America Remote Sensing Center (MARC): Field and Laboratory Enhancement

The purpose of this project was to establish a new hyperspectral remote sensing laboratory at the Mid-America Remote sensing Center (MARC), dedicated to in situ and laboratory measurements of environmental samples and to the manipulation, analysis, and storage of remotely sensed data for environmental monitoring and research in ecological modeling using hyperspectral remote sensing at MARC, one of three research facilities of the Center of Reservoir Research at Murray State University (MSU), a Kentucky Commonwealth Center of Excellence. The equipment purchased, a FieldSpec FR portable spectroradiometer and peripherals, and ENVI hyperspectral data processing software, allowed MARC to provide hands-on experience, education, and training for the students of the Department of Geosciences in quantitative remote sensing using hyperspectral data, Geographic Information System (GIS), digital image processing (DIP), computer, geological and geophysical mapping; to provide field support to the researchers and students collecting in situ and laboratory measurements of environmental data; to create a spectral library of the cover types and to establish a World Wide Web server to provide the spectral library to other academic, state and Federal institutions. Much of the research will soon be published in scientific journals. A World Wide Web page has been created at the web site of MARC. Results of this project are grouped in two categories, education and research accomplishments. The Principal Investigator (PI) modified remote sensing and DIP courses to introduce students to ii situ field spectra and laboratory remote sensing studies for environmental monitoring in the region by using the new equipment in the courses. The PI collected in situ measurements using the spectroradiometer for the ER-2 mission to Puerto Rico project for the Moderate Resolution Imaging Spectrometer (MODIS) Airborne Simulator (MAS). Currently MARC is mapping water quality in Kentucky Lake and vegetation in the Land-Between-the Lakes (LBL) using Landsat-TM data. A Landsat-TM scene of the same day was obtained to relate ground measurements to the satellite data. A spectral library has been created for overstory species in LBL. Some of the methods, such as NPDF and IDFD techniques for spectral unmixing and reduction of effects of shadows in classifications- comparison of hyperspectral classification techniques, and spectral nonlinear and linear unmixing techniques, are being tested using the laboratory.

Cetin, Haluk↗

DOE OSTI2025

Measuring the growth of Alfvén wave parametric decay instability using counter-propagating waves: Theory and simulations

The parametric decay instability (PDI) of Alfvén waves—where a pump Alfvén wave decays into a backward-propagating child Alfvén wave and a forward ion acoustic wave—is a fundamental nonlinear wave-wave interaction and holds significant implications for space and laboratory plasmas. However, to date there has been no direct experimental measurement of PDI. Here, we propose a novel and experimentally viable scheme to quantify the growth of Alfvén wave PDI on a linear device using a large pump Alfvén wave and a small counter-propagating seed Alfvén wave, with the seed-wave frequency tuned to match the backward Alfvén wave generated by standard PDI. Using hybrid simulations, we show that energy transfer from the pump to the seed reduces the latter's spatial damping. By comparing seed-wave amplitudes with and without the pump wave, this damping reduction can be used as a direct and reliable proxy for PDI growth. The method is validated in our simulations across a range of plasma and wave parameters and agrees well with theoretical predictions. Notably, the scheme exhibits no threshold for PDI excitation and is, in principle, readily implementable under current laboratory conditions. Finally, this scheme is a critical step toward solving the challenge of experimentally accessing Alfvén wave PDI and provides an elegant method that may be used to validate fundamental theories of parametric instabilities in controlled laboratory settings.

heliosphere↗

NASA NTRS1998

Bulk-Flow Analysis of Hybrid Thrust Bearings for Advanced Cryogenic Turbopumps

A bulk-flow analysis and computer program for prediction of the static load performance and dynamic force coefficients of angled injection, orifice-compensated hydrostatic/hydrodynamic thrust bearings have been completed. The product of the research is an efficient computational tool for the design of high-speed thrust bearings for cryogenic fluid turbopumps. The study addresses the needs of a growing technology that requires of reliable fluid film bearings to provide the maximum operating life with optimum controllable rotordynamic characteristics at the lowest cost. The motion of a cryogenic fluid on the thin film lands of a thrust bearing is governed by a set of bulk-flow mass and momentum conservation and energy transport equations. Mass flow conservation and a simple model for momentum transport within the hydrostatic bearing recesses are also accounted for. The bulk-flow model includes flow turbulence with fluid inertia advection, Coriolis and centrifugal acceleration effects on the bearing recesses and film lands. The cryogenic fluid properties are obtained from realistic thermophysical equations of state. Turbulent bulk-flow shear parameters are based on Hirs' model with Moody's friction factor equations allowing a simple simulation for machined bearing surface roughness. A perturbation analysis leads to zeroth-order nonlinear equations governing the fluid flow for the thrust bearing operating at a static equilibrium position, and first-order linear equations describing the perturbed fluid flow for small amplitude shaft motions in the axial direction. Numerical solution to the zeroth-order flow field equations renders the bearing flow rate, thrust load, drag torque and power dissipation. Solution to the first-order equations determines the axial stiffness, damping and inertia force coefficients. The computational method uses well established algorithms and generic subprograms available from prior developments. The Fortran9O computer program hydrothrust runs on a Windows 95/NT personal computer. The program, help files and examples are licensed by Texas A&M University Technology License Office. The study of the static and dynamic performance of two hydrostatic/hydrodynamic bearings demonstrates the importance of centrifugal and advection fluid inertia effects for operation at high rotational speeds. The first example considers a conceptual hydrostatic thrust bearing for an advanced liquid hydrogen turbopump operating at 170,000 rpm. The large axial stiffness and damping coefficients of the bearing should provide accurate control and axial positioning of the turbopump and also allow for unshrouded impellers, therefore increasing the overall pump efficiency. The second bearing uses a refrigerant R134a, and its application in oil-free air conditioning compressors is of great technological importance and commercial value. The computed predictions reveal that the LH2 bearing load capacity and flow rate increase with the recess pressure (i.e. increasing orifice diameters). The bearing axial stiffness has a maximum for a recess pressure rati of approx. 0.55. while the axial damping coefficient decreases as the recess pressure ratio increases. The computer results from three flow models are compared. These models are a) inertialess, b) fluid inertia at recess edges only, and c) full fluid inertia at both recess edges and film lands. The full inertia model shows the lowest flow rates, axial load capacity and stiffness coefficient but on the other hand renders the largest damping coefficients and inertia coefficients. The most important findings are related to the reduction of the outflow through the inner radius and the appearance of subambient pressures. The performance of the refrigerant hybrid thrust bearing is evaluated at two operating speeds and pressure drops. The computed results are presented in dimensionless form to evidence consistent trends in the bearing performance characteristics. As the applied axial load increases, the bearing film thickness and flow rate decrease while the recess pressure increases. The axial stiffness coefficient shows a maximum for a certain intermediate load while the damping coefficient steadily increases. The computed results evidence the paramount of centrifugal fluid inertia at low recess pressures (i.e. low loads), and where there is actually an inflow through the bearing inner diameter, accompanied by subambient pressures just downstream of the bearing recess edge. These results are solely due to centrifugal fluid inertia and advection transport effects. Recommendations include the extension of the computer program to handle flexure pivot tilting pad hybrid bearings and the ability to calculate moment coefficients for shaft angular misalignments.

SanAndres, Luis↗

DOE OSTI2024

Performance Improvements of the Griffin Solvers in FY24

The Griffin code is a MOOSE-based reactor physics application jointly developed by Idaho National Laboratory and Argonne National Laboratory under the Department of Energy Office of Nuclear Energy Nuclear Energy Advanced Modeling and Simulation Program. This fiscal year, we have made significant efforts to improve the performance of transport solver options and cross-section generation for the efficient use of Griffin in advanced reactor applications. For the HFEM-PN solver, the residual evaluations of HFEM kernels were optimized by utilizing the pre- computed averaged cross sections for individual elements. Numerical integration involving the evaluation of basis functions at quadrature points was bypassed by facilitating precomputed element mass matrices for response matrices. Red-black iterations were improved by introducing a new generalized minimum residual based solver. The memory usage of response matrix storage was significantly reduced by applying basis function rotations on interfaces and calculating volumetric odd-parity moments on the fly. Additionally, the adjoint flux and transient calculation capabilities of the HFEM-PN solver were successfully implemented and verified using the TWIGL benchmark problem. For the DFEM-SN solver, memory footprint and computation time were significantly reduced by not treating angular flux vectors as the MOOSE nonlinear system vectors. Specifically for IQS, scalar adjoint weighting was introduced to further eliminate angular adjoint flux storage in the MOOSE auxiliary system. It was demonstrated through the three-dimensional Advanced Burner Test Reactor core problem that the memory usage for transient calculations with the IQS method was reduced by over 7.5× compared to before the optimizations. For the self-shielding application programming interface, a new double-heterogeneity treatment method, named the Bell Function-Based Analytic Two-Region Slowing Down Method, was developed to efficiently flux-volume homogenize TRISO particles with the matrix. Additionally, optimizations were made to hyper- fine group (HFG) slowing down calculations by pretabulating collision probability coefficients and grouping isotopes, significantly reducing the computational time for calculating scattering sources per HFG. Lastly, the pin power reconstruction module was extended to account for temporal behavior in a microreactor analysis problem, specifically for a control drum transient. Verification tests for each of these improvements demonstrated significant performance enhancements and memory reduction.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS2020

Feasibility Study of SDAS Instrumentation's Ability to Identify Mobile Launcher (ML)/Crawler-Transporter (CT) Modes During Rollout Operations

The Space Launch System (SLS) and its Mobile Launcher (ML) will be transported to the launch pad via the Crawler-Transporter (CT) system. Rollout (i.e., transportation) loads produce structural loads on the integrated SLS/Orion Multi-Purpose Crew Vehicle (MPCV) launch vehicle which are of a concern with respect to fatigue. As part of the risk reduction process and in addition to the modal building block test approach that has been adopted by the SLS Program, acceleration data will be obtained during rollout for use in modal parameter estimation. There are several occurrences where the ML/CT will be transported either into the Vertical Assembly Building (VAB) or to the launch pad and back without the SLS stack as part of the Kennedy Space Center (KSC) Exploration Ground Systems (EGS) Integrated Test and Checkout (ITCO). NASA KSC EGS has instrumentation installed on both the ML and CT to record data during rollout, at the launch pad, and during liftoff. The EGS instrumentation on the ML, which includes accelerometers, is referred to as the Sensor Data Acquisition System (SDAS). The EGS instrumentation on the CT, which also includes accelerometers, is referred to as the CT Data Acquisition System (CTDAS). The forces and accelerations applied to the ML and CT during a rollout event will be higher than any of the planned building block modal tests. This can be very beneficial in helping identify nonlinear behavior in the structure. Developing modal parameters from the same test hardware in multiple boundary conditions and under multiple levels of excitation is a key step in developing a well correlated FEM. The purpose of this study was three fold. First, determine the target modes of the ML/CT in its rollout configuration. Second, determine if the test degrees of freedom (DOF) corresponding to the layout of the SDAS/CTDAS accelerometers (i.e. position and orientation) is sufficient to identify the target modes. Third, determine if the Generic Rollout Forcing Functions (GRFF's) is sufficient for identifying the ML/CT target modes accounting for variations in CT speed, modal damping, and sensor/ambient background noise levels. The finding from the first part of this study identified 28 target modes of the ML/CT rollout configuration based upon Modal Effective Mass Fractions (MEFF) and engineering judgement. The finding from the second part of this study showed that the SDAS/CTDAS accelerometers (i.e. position and orientation) are able to identify a sufficient number of the target modes to support model correlation of the ML/CT FEM. The finding from the third part of this study confirms the GRFFs sufficiently excite the ML/CT such that varying quantities of the defined target modes should be able to be extracted when utilizing an Experimental Modal Analysis (EMA) Multi-Input Multi-Output (MIMO) analysis approach. An EMA analysis approach was used because Operational Modal Analysis (OMA) tools were not available and the GRFFs were sufficiently uncorrelated. Two key findings from this third part of the study are that the CT speed does not show a significant impact on the ability to extract the modal parameters and that keeping the ambient background noise observed at each accelerometer location at or below 30 µgrms is essential to the success of this approach.

Winkel, James P.↗

NASA NTRS2005

Two-Photon Fluorescence Microscope for Microgravity Research

A two-photon fluorescence microscope has been developed for the study of biophysical phenomena. Two-photon microscopy is a novel form of laser-based scanning microscopy that enables three-dimensional imaging without many of the problems inherent in confocal microscopy. Unlike one-photon optical microscopy, two-photon microscopy utilizes the simultaneous nonlinear absorption of two near-infrared photons. However, the efficiency of two-photon absorption is much lower than that of one-photon absorption, so an ultra-fast pulsed laser source is typically employed. On the other hand, the critical energy threshold for two-photon absorption leads to fluorophore excitation that is intrinsically localized to the focal volume. Consequently, two-photon microscopy enables optical sectioning and confocal performance without the need for a signal-limiting pinhole. In addition, there is a reduction (relative to one-photon optical microscopy) in photon-induced damage because of the longer excitation wavelength. This reduction is especially advantageous for in vivo studies. Relative to confocal microscopy, there is also a reduction in background fluorescence, and, because of a reduction in Rayleigh scattering, there is a 4 increase of penetration depth. The prohibitive cost of a commercial two-photon fluorescence-microscope system, as well as a need for modularity, has led to the construction of a custom-built system (see Figure 1). This system includes a coherent mode-locked titanium: sapphire laser emitting 120-fs-duration pulses at a repetition rate of 80 MHz. The pulsed laser has an average output power of 800 mW and a wavelength tuning range of 700 to 980 nm, enabling the excitation of a variety of targeted fluorophores. The output from the laser is attenuated, spatially filtered, and then directed into a confocal scanning head that has been modified to provide for side entry of the laser beam. The laser output coupler has been replaced with a dichroic filter that reflects the longer-wavelength excitation light and passes the shorter-wavelength fluorescence light. Also, the confocal pinhole has been removed to increase the signal strength. The laser beam is scanned by a twoperpendicular- axis pair of galvanometer mirrors through a pupil transfer lens into the side port of an inverted microscope. Finally, the beam is focused by a 63-magnification, 1.3-numerical- aperture oil-immersion objective lens onto a specimen. The pupil transfer lens serves to match the intermediate image planes of the scanning head and the microscope, and its location is critical. In order to maximize the quality of the image, (that is, the point spread function of the objective lens for all scan positions), the entire system was modeled in optical-design software, and the various free design parameters (the parameters of the spatial-filter components as well as the separations of all of the system components) were determined through an iterative optimization process. A modular design was chosen to facilitate access to the optical train for future fluorescence correlation spectroscopy and fluorescence-lifetime experiments.

Fischer, David G.↗

NASA NTRS1999

Fractal analysis of heart rate dynamics as a predictor of mortality in patients with depressed left ventricular function after acute myocardial infarction. TRACE Investigators. TRAndolapril Cardiac Evaluation

A number of new methods have been recently developed to quantify complex heart rate (HR) dynamics based on nonlinear and fractal analysis, but their value in risk stratification has not been evaluated. This study was designed to determine whether selected new dynamic analysis methods of HR variability predict mortality in patients with depressed left ventricular (LV) function after acute myocardial infarction (AMI). Traditional time- and frequency-domain HR variability indexes along with short-term fractal-like correlation properties of RR intervals (exponent alpha) and power-law scaling (exponent beta) were studied in 159 patients with depressed LV function (ejection fraction <35%) after an AMI. By the end of 4-year follow-up, 72 patients (45%) had died and 87 (55%) were still alive. Short-term scaling exponent alpha (1.07 +/- 0.26 vs 0.90 +/- 0.26, p <0.001) and power-law slope beta (-1.35 +/- 0.23 vs -1.44 +/- 0.25, p <0.05) differed between survivors and those who died, but none of the traditional HR variability measures differed between these groups. Among all analyzed variables, reduced scaling exponent alpha (<0.85) was the best univariable predictor of mortality (relative risk 3.17, 95% confidence interval 1.96 to 5.15, p <0.0001), with positive and negative predictive accuracies of 65% and 86%, respectively. In the multivariable Cox proportional hazards analysis, mortality was independently predicted by the reduced exponent alpha (p <0.001) after adjustment for several clinical variables and LV function. A short-term fractal-like scaling exponent was the most powerful HR variability index in predicting mortality in patients with depressed LV function. Reduction in fractal correlation properties implies more random short-term HR dynamics in patients with increased risk of death after AMI.

Non-NASA Center↗

NASA NTRS2000

Optimization of Car Body under Constraints of Noise, Vibration, and Harshness (NVH), and Crash

To be competitive on the today's market, cars have to be as light as possible while meeting the Noise, Vibration, and Harshness (NVH) requirements and conforming to Government-man dated crash survival regulations. The latter are difficult to meet because they involve very compute-intensive, nonlinear analysis, e.g., the code RADIOSS capable of simulation of the dynamics, and the geometrical and material nonlinearities of a thin-walled car structure in crash, would require over 12 days of elapsed time for a single design of a 390K elastic degrees of freedom model, if executed on a single processor of the state-of-the-art SGI Origin2000 computer. Of course, in optimization that crash analysis would have to be invoked many times. Needless to say, that has rendered such optimization intractable until now. The car finite element model is shown. The advent of computers that comprise large numbers of concurrently operating processors has created a new environment wherein the above optimization, and other engineering problems heretofore regarded as intractable may be solved. The procedure, shown, is a piecewise approximation based method and involves using a sensitivity based Taylor series approximation model for NVH and a polynomial response surface model for Crash. In that method the NVH constraints are evaluated using a finite element code (MSC/NASTRAN) that yields the constraint values and their derivatives with respect to design variables. The crash constraints are evaluated using the explicit code RADIOSS on the Origin 2000 operating on 256 processors simultaneously to generate data for a polynomial response surface in the design variable domain. The NVH constraints and their derivatives combined with the response surface for the crash constraints form an approximation to the system analysis (surrogate analysis) that enables a cycle of multidisciplinary optimization within move limits. In the inner loop, the NVH sensitivities are recomputed to update the NVH approximation model while keeping the Crash response surface constant. In every outer loop, the Crash response surface approximation is updated, including a gradual increase in the order of the response surface and the response surface extension in the direction of the search. In this optimization task, the NVH discipline has 30 design variables while the crash discipline has 20 design variables. A subset of these design variables (10) are common to both the NVH and crash disciplines. In order to construct a linear response surface for the Crash discipline constraints, a minimum of 21 design points would have to be analyzed using the RADIOSS code. On a single processor in Origin 2000 that amount of computing would require over 9 months! In this work, these runs were carried out concurrently on the Origin 2000 using multiple processors, ranging from 8 to 16, for each crash (RADIOSS) analysis. Another figure shows the wall time required for a single RADIOSS analysis using varying number of processors, as well as provides a comparison of 2 different common data placement procedures within the allotted memories for each analysis. The initial design is an infeasible design with NVH discipline Static Torsion constraint violations of over 10%. The final optimized design is a feasible design with a weight reduction of 15 kg compared to the initial design. This work demonstrates how advanced methodology for optimization combined with the technology of concurrent processing enables applications that until now were out of reach because of very long time-to-solution.

Kodiyalam, Srinivas↗