Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reduced basis method”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

Reduced Stress Tensor and Dissipation and the Transport of Lamb Vector

We develop a methodology to ensure that the stress tensor, regardless of its number of independent components, can be reduced to an exactly equivalent one which has the same number of independent components as the surface force. It is applicable to the momentum balance if the shear viscosity is constant. A direct application of this method to the energy balance also leads to a reduction of the dissipation rate of kinetic energy. Following this procedure, significant saving in analysis and computation may be achieved. For turbulent flows, this strategy immediately implies that a given Reynolds stress model can always be replaced by a reduced one before putting it into computation. Furthermore, we show how the modeling of Reynolds stress tensor can be reduced to that of the mean turbulent Lamb vector alone, which is much simpler. As a first step of this alternative modeling development, we derive the governing equations for the Lamb vector and its square. These equations form a basis of new second-order closure schemes and, we believe, should be favorably compared to that of traditional Reynolds stress transport equation.

Wu, Jie-Zhi↗
NASA NTRS2012

A New Approach in Force-Limited Vibration Testing of Flight Hardware

The force-limited vibration test approaches discussed in NASA-7004C were developed to reduce overtesting associated with base shake vibration tests of aerospace hardware where the interface responses are excited coherently. This handbook outlines several different methods of specifying the force limits. The rationale for force limiting is based on the disparity between the impedances of typical aerospace mounting structures and the large impedances of vibration test shakers when the interfaces in general are coherently excited. Among these approaches, the semi-empirical method is presently the most widely used method to derive the force limits. The inclusion of the incoherent excitation of the aerospace structures at mounting interfaces has not been accounted for in the past and provides the basis for more realistic force limits for qualifying the hardware using shaker testing. In this paper current methods for defining the force limiting specifications discussed in the NASA handbook are reviewed using data from a series of acoustic and vibration tests. A new approach based on considering the incoherent excitation of the structural mounting interfaces using acoustic test data is also discussed. It is believed that the new approach provides much more realistic force limits that may further remove conservatism inherent in shaker vibration testing not accounted for by methods discussed in the NASA handbook. A discussion on using FEM/BEM analysis to obtain realistic force limits for flight hardware is provided.

force limiting responses↗
DOE OSTI2022

Modeling the Spectral Diversity of Quasars in the Sixteenth Data Release from the Sloan Digital Sky Survey

Abstract We present a new approach to capturing the broad diversity of emission-line and continuum properties in quasar spectra. We identify populations of spectrally similar quasars through pixel-level clustering on 12,968 high signal-to-noise ratio (S/N) spectra from the Sloan Digital Sky Survey (SDSS) in the redshift range of 1.57 < z < 2.4. Our clustering analysis finds 396 quasar spectra that are not assigned to any population, 15 misclassified spectra, and 6 quasars with incorrect redshifts. We compress the quasar populations into a library of 684 high-S/N composite spectra, anchored in redshift space by the Mg ii emission line. Principal component analysis on the library results in an eigenspectrum basis spanning 1067–4007 Å. We model independent samples of SDSS quasar spectra with the eigenbasis, allowing for a free redshift parameter. Our models achieve a median reduced χ 2 on non–broad absorption line quasar spectra that is reduced by 8.5% relative to models using the eigenspectra from the SDSS spectroscopic pipeline. A significant contribution to the relative improvement is from the ability to reconstruct the range of emission-line variation. The redshift estimates from our model are consistent with the Mg ii emission-line redshift with an average offset that displays 51.4% less redshift-dependent variation relative to the SDSS eigenspectra. Our method for developing quasar spectra models can improve automated classification and predict the intrinsic spectrum in regions affected by intervening absorbers such as Ly α , C iv , and Mg ii , thus benefiting studies of large-scale structure.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

CMB-S4: Foreground-cleaning Pipeline Comparison for Measuring Primordial Gravitational Waves

We compare multiple foreground-cleaning pipelines for estimating the tensor-to-scalar ratio, r, using simulated maps of the planned CMB-S4 experiment within the context of the South Pole Deep Patch. To evaluate robustness, we analyze bias and uncertainty on r across various foreground suites using map-based simulations. The foreground-cleaning methods include: a parametric maximum likelihood approach applied to auto- and cross-power spectra between frequency maps; a map-based parametric maximum-likelihood method; and a harmonic-space internal linear combination using frequency maps. We summarize the conceptual basis of each method to highlight their similarities and differences. To better probe the impact of foreground residuals, we implement an iterative internal delensing step, leveraging a map-based pipeline to generate a lensing B-mode template from the large aperture telescope frequency maps. Our results show that the performance of the three approaches is comparable for simple and intermediate-complexity foregrounds, with σ(r) ranging from 3–5 ×10 −4 . However, biases at the 1σ–2σ level appear when analyzing more complex forms of foreground emission. By extending the baseline pipelines to marginalize over foreground residuals, we demonstrate that contamination can be reduced to within statistical uncertainties, albeit with a pipeline-dependent impact on σ(r), which translates to a detection significance between 2σ and 4σ for an input value of r = 0.003. These findings suggest varying levels of maturity among the tested pipelines, with the auto- and cross-spectra-based approach demonstrating the best stability and overall performance. Moreover, given the extremely low noise levels, mutual validation of independent foreground-cleaning pipelines is essential to ensure the robustness of any potential detection.

astronomy data analysis↗
NASA NTRS1989

Aeroservoelastic modeling and applications using minimum-state approximations of the unsteady aerodynamics

The theoretical basis of flexible-aircraft modeling techniques encompassing aerodynamic, control, and elastic-structure effects is investigated analytically, with a focus on methods which employ minimum-state approximations for the unsteady aerodynamics. Rational-function approximations to generalized aerodynamic forces are reviewed; constraints and lag-coefficient optimization are explained; the problem of physical weighting in the minimum-state equations of motion is examined; and results of typical analyses from the NASA Active Flexible Wing project (Perry et al., 1988) are presented in extensive tables and graphs and discussed in detail. The minimum-state approach is shown to produce accurate models at significantly reduced computation costs.

Tiffany, Sherwood H.↗
NASA NTRS1992

A fundamental study of suction for Laminar Flow Control (LFC)

This report covers the period forming the first year of the project. The aim is to experimentally investigate the effects of suction as a technique for Laminar Flow Control. Experiments are to be performed which require substantial modifications to be made to the experimental facility. Considerable effort has been spent developing new high performance constant temperature hot-wire anemometers for general purpose use in the Fluid Mechanics Laboratory. Twenty instruments have been delivered. An important feature of the facility is that it is totally automated under computer control. Unprecedently large quantities of data can be acquired and the results examined using the visualization tools developed specifically for studying the results of numerical simulations on graphics works stations. The experiment must be run for periods of up to a month at a time since the data is collected on a point-by-point basis. Several techniques were implemented to reduce the experimental run-time by a significant factor. Extra probes have been constructed and modifications have been made to the traverse hardware and to the real-time experimental code to enable multiple probes to be used. This will reduce the experimental run-time by the appropriate factor. Hot-wire calibration drift has been a frustrating problem owing to the large range of ambient temperatures experienced in the laboratory. The solution has been to repeat the calibrations at frequent intervals. However the calibration process has consumed up to 40 percent of the run-time. A new method of correcting the drift is very nearly finalized and when implemented it will also lead to a significant reduction in the experimental run-time.

Watmuff, Jonathan H.↗
DOE OSTI2023

Rapid AI-based Dissection of Ashes using Raman and XRF Spectroscopy (RADAR-X)

Waste-to-energy (WTE) facilities incinerate ~35 million tons of municipal solid waste annually in the United States. The incineration process reduces the mass and the volume of the waste fraction by over 75 and 95%, respectively. The fraction left after incineration remains as ash residues and is referred to as WTE ash, compromising of bottom and fly ash. In the United States, ~10 million tons of WTE ashes are generated annually and predominantly landfilled because of the lack of secondary end-use pathways. This incurs a significant financial burden (landfilling costs) on U.S. WTE facilities and also results in the loss of materials to the landfill. The primary objective of this research is to understand better the elemental and mineralogical composition of WTE ashes from diverse origins and find composition dependent upcycling pathways for diverting these ashes from landfills. This primary objective was addressed through three research tasks –(Task I) An AI-led Multi-Modal Approach for Compositional Analysis, (Task II) Developing a Dissolution-Based Test for Real Time Analysis, and (Task III) Establishing composition-dependent end uses. The chemical composition of WTE ash is dependent on two factors, i.e., the input waste composition and the operational parameters of a WTE facility (combustion conditions). Amongst these two factors, the input waste composition will likely show spatial and temporal variation. As a result, the chemical composition of WTE ash will also fluctuate. To understand the spatial and temporal variation in WTE ash composition, in Task I, we collected 128 ash samples (62 bottom ash and 66 fly ash samples) from 11 WTE facilities located in 11 U.S. states and characterized them via X-ray Fluorescence, powder X-ray Diffraction, and Raman Spectroscopy. The findings from this extensive characterization work indicated that the key elements in WTE fly ashes are Ca, Cl (greater than 10 wt. %), Si, S, K, Zn ( between 1 and 10 wt. %), Mg, Al, P, Ti, Fe, Cu, Br, and Pb (between 0.1 and 1 wt. %). Similarly, the key elements in WTE bottom ash fraction finer than 45μm are Ca (greater than 10 wt. %), Mg, Al, Si, S, Cl, K, Ti, Fe, Zn (between 1 and 10 wt. %), P, V, Cr, Mn, Cu, Br, and Pb (between 0.1 and 1 wt. %). Here, we note that the dominant fraction of WTE bottom ash is the coarse fraction. The coarse WTE bottom ash fraction (rich in silicon) was not characterized in this study because of excessive grinding requirements and their unsuitability as a supplementary cementitious material due to their coarse nature. The elements in WTE bottom ashes are present as calcite, anhydrite, vaterite, hydroxyapatite, quartz, bassanite, gehlenite, akermanite, hydrocalumite, and portlandite. Similarly, the mineralogical species present in WTE fly ashes are calcium chloride hydroxide, halite, calcite, anhydrite, sylvite, hydrocalumite, vaterite, hannebachite, quartz, and bassanite. Temporal variation in ash composition may also result in significant fluctuations in chemical compositions. Therefore, a WTE facility may need to monitor the ash composition (elemental and mineralogical composition) in real time. In Task I, we evaluated the possibility of using a portable X-ray fluorescence (XRF) spectrometer to monitor the elemental composition in real-time. Specifically, we collected XRF data on identical specimens via a portable XRF spectrometer (low-end) and a lab-based benchtop XRF spectrometer (high-end). The collected data was used to train an A.I. algorithm (portable XRF data as an input and benchtop XRF data as an output) to predict accurate elemental composition using portable XRF data. Finally, we developed a 2-minute photobleaching protocol to monitor the mineralogical characteristics of WTE ashes via Raman spectroscopy. Overall, the activities in Task I improved our understanding of ash composition and developed techniques to monitor elemental and mineralogical composition in near real-time. Based on the findings of Task I, we find that WTE ashes exhibit wide variability in mineralogy. For ICP-based elemental analysis, all the mineralogical species in WTE ashes must be brought into solution. This is traditionally accomplished with acid digestion using a combination of multiple acids. However, acid digestion with multiple acids is time-consuming and often fails to ensure complete digestion of the ash matrix. To address this limitation, in Task II, we developed an alkali-fusion-based digestion protocol using a combination of lithium tetraborate, lithium metaborate, and their combinations as possible alkali fluxes for digesting WTE ashes entirely and rapidly. The validity of the developed method was evaluated on two standard ash specimens, i.e., SRM 1633c coal fly ash and BCR-176R incineration fly ash specimen. The findings suggest that the developed protocol can ensure complete digestion of elements such as Al, Ba, Ca, Cr, Cu, Mg, Mn, P, Sr, V, Zn, Be, K, and rare earth elements. The recent changes in the energy market towards renewables and increased metal recycling have resulted in reduced supplies of supplementary cementitious materials (coal fly ash and slag). Therefore, in Task III, we evaluated the possibility of employing WTE ashes as SCMs. As the chemical composition of WTE ashes varies temporally (on an hourly basis), there was also a need to develop tests that can evaluate the suitability of material to act as supplementary cementitious material rapidly, i.e., in a few minutes. Therefore, in Task II, we also developed a rapid test to assess the suitability of a material to act as an SCM in 5 minutes. This represents a significant advance over the state-of-the-art R 3 test, which takes ~144 hours. This test was initially validated on amorphous aluminosilicates, such as calcined clays, and could be extended to evaluate other industrial by-products, such as WTE ashes. In Task III, we evaluated the possibility of employing WTE ashes for two applications, i.e., as an SCM and a lime substitute for clay stabilization. The findings from Task I indicated that WTE ashes are enriched in chlorine and, therefore, cannot be used directly as an SCM due to corrosion-related risks and altered hydration kinetics. Accordingly, we developed an ash treatment protocol to reduce the solubility of chlorine-containing species in WTE ashes. The developed treatment protocol also immobilized lead in certain mineral forms. As a result of the treatment, WTE ashes can be used as SCMs without any corrosion or heavy metal leaching concerns. The second application examined in this study was clay stabilization. WTE ashes are calcium-rich and can be an adequate lime replacement for clay stabilization. Our findings reveal that the sum of the concentrations of Ca(OH) 2 and CaClOH controls the clay stabilization capability of WTE ashes. In summary, in this work, we evaluated the elemental and mineralogical characteristics of U.S. WTE ashes from diverse origins and developed tests to evaluate the chemical characteristics of these ashes in real time through a portable XRF and a benchtop Raman spectrometer. Based on the chemical characteristics of these ashes, we developed an ash treatment process to enable the use of WTE ashes as an SCM and also evaluated the possibility of employing these ashes for clay stabilization. Overall, the findings from this work enables the diversion of WTE ashes from landfills for multiple end-uses, i.e., as an SCM or a lime substitute for clay stabilization.

36 MATERIALS SCIENCE↗
NASA NTRS1977

NASTRAN hydroelastic modal studies. Volume 3: NASTRAN 3-D hydroelastic analysis and MESHGEN user's manuals

The method used to formulate the fluid/structure equations is described in the updates to the NASTRAN Theoretical Manual. The basis for defining the fluid are three dimensional finite elements connected to fluid grid points defining the Eulerian pressure at a point fixed in space. The use of a pressure single degree of freedom at each point rather than three displacements allows a finer mesh of elements with a reasonable matrix order. In the formulation of the fluid/structure system the interior fluid degrees of freedom are transformed and removed from the solution matrices. The eigenvalues of the combination are extracted from small, fully dense, symmetric mass and stiffness matrices, efficiently processed with the Givens method. The solution matrices are defined by the free surface displacements and the reduced structure coordinates.

Source record↗
NASA NTRS1986

ECUT (Energy Conversion and Utilization Technologies Program). Biocatalysis Project

Presented are the FY 1985 accomplishments, activities, and planned research efforts of the Biocatalysis Project of the U.S. Department of Energy, Energy Conversion and Utilization Technologies (ECUT) Program. The Project's technical activities were organized as follows: In the Molecular Modeling and Applied Genetics work element, research focused on (1) modeling and simulation studies to establish the physiological basis of high temperature tolerance in a selected enzyme and the catalytic mechanisms of three species of another enzyme, and (2) determining the degree of plasmid amplification and stability of several DNA bacterial strains. In the Bioprocess Engineering work element, research focused on (1) studies of plasmid propagation and the generation of models, (2) developing methods for preparing immobilized biocatalyst beads, and (3) developing an enzyme encapsulation method. In the Process Design and Analysis work element, research focused on (1) further refinement of a test case simulation of the economics and energy efficiency of alternative biocatalyzed production processes, (2) developing a candidate bioprocess to determine the potential for reduced energy consumption and facility/operating costs, and (3) a techno-economic assessment of potential advancements in microbial ammonia production.

Source record↗
NASA NTRS1992

Compiler analysis for irregular problems in FORTRAN D

We developed a dataflow framework which provides a basis for rigorously defining strategies to make use of runtime preprocessing methods for distributed memory multiprocessors. In many programs, several loops access the same off-processor memory locations. Our runtime support gives us a mechanism for tracking and reusing copies of off-processor data. A key aspect of our compiler analysis strategy is to determine when it is safe to reuse copies of off-processor data. Another crucial function of the compiler analysis is to identify situations which allow runtime preprocessing overheads to be amortized. This dataflow analysis will make it possible to effectively use the results of interprocedural analysis in our efforts to reduce interprocessor communication and the need for runtime preprocessing.

Vonhanxleden, Reinhard↗
NASA NTRS1998

An Improved Magnetic Field Model for Jupiter's Inner Magnetosphere Using a Microsignature of Amalthea

Observation of a particle absorption microsignature of the Jovian satellite Amalthea during the Pioneer II close flyby of Jupiter on December 3, 1974, has been described by McKibben et al. The microsignature was also observed by the University of Iowa/Pioneer II instrument in the distribution of protons in the kinetic energy range 0.61 less than T(sub p) less than 3.41 MeV but has not been previously reported. The finer time resolution and superior data quality of the latter observations provide a fresh basis for assessing the accuracy with which various published magnetic field models of Jupiter's magnetic field models describe Jupiter's magnetic field for radial distances of the order of or less than several planetary radii. The expected time of occurrence of the minimum of Amalthea's microsignature was calculated for each model and compared to the time of the observed minimum. The discrepancy between these two times was significant for each of the published models, but it could be reduced to zero by interpolating between two of the closely related models. The spherical harmonic coefficients of this combination model are tabulated. Also included is a novel method for calculating the intensity-time profile of the microsignature. In a later paper, the improved model is utilized to analyze the energetic particle measurements in Jupiter's innermost magnetosphere by the Galileo entry probe in December 1995.

Randall, Bruce A.↗
DOE OSTI2025

Exploiting a Shortcoming of Coupled-Cluster Theory: The Extent of Non-Hermiticity as a Diagnostic Indicator of Computational Accuracy

The fundamental non-Hermitian nature of the forms of the coupled-cluster (CC) theory widely used in quantum chemistry has usually been viewed as a negative, but the present paper shows how this can be used to an advantage. Specifically, the non-symmetric nature of the reduced one-particle density matrix (in the molecular orbital basis) is advocated as a diagnostic indicator of computational quality. In the limit of the full coupled-cluster theory [which is equivalent to full configuration interaction (FCI)], the electronic wave function and correlation energy are exact within a given one-particle basis set, and the symmetric character of the exact density matrix is recovered. The extent of the density matrix asymmetry is shown to provide a measure of “how difficult the problem is” (like the well-known T 1 diagnostic), but its variation with the level of theory also gives information about “how well this particular method works”, irrespective of the difficulty of the problem at hand. The proposed diagnostic is described and applied to a select group of small molecules, and an example of its overall utility for the practicing quantum chemist is illustrated through its application to the beryllium dimer (Be 2 ). Future application of this idea to excited states, open-shell systems, and symmetry-breaking problems and an extension of the method to the two-particle density are then proposed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1988

Element-by-element and implicit-explicit finite element formulations for computational fluid dynamics

Preconditioner algorithms to reduce the computational effort in FEM analyses of large-scale fluid-dynamics problems are presented. A general model problem is constructed on the basis of the convection-diffusion equation and the two-dimensional vorticity/stream-function formulation of the Navier-Stokes equations; this problem is then analyzed using element-by-element, implicit-explicit, and adaptive implicit-explicit approximation schemes. Numerical results for the two-dimensional advection and rigid-body rotation of a cosine hill, flow past a circular cylinder, and driven cavity flow are presented in extensive graphs and shown to be in good agreement with those obtained using implicit methods.

Tezduyar, T. E.↗
NASA NTRS1993

Global/local methods for probabilistic structural analysis

A probabilistic global/local method is proposed to reduce the computational requirements of probabilistic structural analysis. A coarser global model is used for most of the computations with a local more refined model used only at key probabilistic conditions. The global model is used to establish the cumulative distribution function (cdf) and the Most Probable Point (MPP). The local model then uses the predicted MPP to adjust the cdf value. The global/local method is used within the advanced mean value probabilistic algorithm. The local model can be more refined with respect to the g1obal model in terms of finer mesh, smaller time step, tighter tolerances, etc. and can be used with linear or nonlinear models. The basis for this approach is described in terms of the correlation between the global and local models which can be estimated from the global and local MPPs. A numerical example is presented using the NESSUS probabilistic structural analysis program with the finite element method used for the structural modeling. The results clearly indicate a significant computer savings with minimal loss in accuracy.

Millwater, H. R.↗
DOE OSTI2023

Constraint energy minimizing generalized multiscale finite element method for multi-continuum Richards equations

In fluid flow simulation, the multi-continuum model is a useful strategy. When the heterogeneity and contrast of coefficients are high, the system becomes multiscale, and some kinds of reduced order methods are demanded. Combining these techniques with nonlinearity, we will consider in this paper a dual-continuum model which is generalized as a multi-continuum model for a coupled system of nonlinear Richards equations as unsaturated flows, in complex heterogeneous fractured porous media; and we will solve it by a novel multiscale approach utilizing the constraint energy minimizing generalized multiscale finite element method (CEM-GMsFEM). In particular, such a nonlinear system will be discretized in time and then linearized by Picard iteration (whose global convergence is proved theoretically). Subsequently, we tackle the resulting linearized equations by the CEM-GMsFEM and obtain proper offline multiscale basis functions to span the multiscale space (which contains the pressure solution). More specifically, we first introduce two new sources of samples, and the GMsFEM is used over each coarse block to build local auxiliary multiscale basis functions via solving local spectral problems, that are crucial for detecting high-contrast channels. Second, per oversampled coarse region, local multiscale basis functions are created through the CEM as constrainedly minimizing an energy functional. Various numerical tests for our approach reveal that the error converges with the coarse-grid size and that only few oversampling layers as well as basis functions are needed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Comparison of Equilibrium and Nonequilibrium Approaches for Relative Binding Free Energy Predictions

Alchemical relative binding free energy calculations have recently found important applications in drug optimization. A series of congeneric compounds are generated from a preidentified lead compound, and their relative binding affinities to a protein are assessed in order to optimize candidate drugs. While methods based on equilibrium thermodynamics have been extensively studied, an approach based on nonequilibrium methods has recently been reported together with claims of its superiority. However, these claims pay insufficient attention to the basis and reliability of both methods. Here we report a comparative study of the two approaches across a large data set, comprising more than 500 ligand transformations spanning in excess of 300 ligands binding to a set of 14 diverse protein targets. Ensemble methods are essential to quantify the uncertainty in these calculations, not only for the reasons already established in the equilibrium approach but also to ensure that the nonequilibrium calculations reside within their domain of validity. If and only if ensemble methods are applied, we find that the nonequilibrium method can achieve accuracy and precision comparable to those of the equilibrium approach. Compared to the equilibrium method, the nonequilibrium approach can reduce computational costs but introduces higher computational complexity and longer wall clock times. There are, however, cases where the standard length of a nonequilibrium transition is not sufficient, necessitating a complete rerun of the entire set of transitions. This significantly increases the computational cost and proves to be highly inconvenient during large-scale applications. Our findings provide a key set of recommendations that should be adopted for the reliable implementation of nonequilibrium approaches to relative binding free energy calculations in ligand-protein systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 NTRS1996

User Interface Developed for Controls/CFD Interdisciplinary Research

The NASA Lewis Research Center, in conjunction with the University of Akron, is developing analytical methods and software tools to create a cross-discipline "bridge" between controls and computational fluid dynamics (CFD) technologies. Traditionally, the controls analyst has used simulations based on large lumping techniques to generate low-order linear models convenient for designing propulsion system controls. For complex, high-speed vehicles such as the High Speed Civil Transport (HSCT), simulations based on CFD methods are required to capture the relevant flow physics. The use of CFD should also help reduce the development time and costs associated with experimentally tuning the control system. The initial application for this research is the High Speed Civil Transport inlet control problem. A major aspect of this research is the development of a controls/CFD interface for non-CFD experts, to facilitate the interactive operation of CFD simulations and the extraction of reduced-order, time-accurate models from CFD results. A distributed computing approach for implementing the interface is being explored. Software being developed as part of the Integrated CFD and Experiments (ICE) project provides the basis for the operating environment, including run-time displays and information (data base) management. Message-passing software is used to communicate between the ICE system and the CFD simulation, which can reside on distributed, parallel computing systems. Initially, the one-dimensional Large-Perturbation Inlet (LAPIN) code is being used to simulate a High Speed Civil Transport type inlet. LAPIN can model real supersonic inlet features, including bleeds, bypasses, and variable geometry, such as translating or variable-ramp-angle centerbodies. Work is in progress to use parallel versions of the multidimensional NPARC code.

Source record↗