Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reduced Order 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,693 records · Page 94
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

A staggered-grid finite-difference scheme optimized in the time–space domain for modeling scalar-wave propagation in geophysical problems

For modeling scalar-wave propagation in geophysical problems using finite-difference schemes, optimizing the coefficients of the finite-difference operators can reduce numerical dispersion. Most optimized finite-difference schemes for modeling seismic-wave propagation suppress only spatial but not temporal dispersion errors. We develop a novel optimized finite-difference scheme for numerical scalar-wave modeling to control dispersion errors not only in space but also in time. Our optimized scheme is based on a new stencil that contains a few more grid points than the standard stencil. We design an objective function for minimizing relative errors of phase velocities of waves propagating in all directions within a given range of wavenumbers. Dispersion analysis and numerical examples demonstrate that our optimized finite-difference scheme is computationally up to 2.5 times faster than the optimized schemes using the standard stencil to achieve the similar modeling accuracy for a given 2D or 3D problem. Compared with the high-order finite-difference scheme using the same new stencil, our optimized scheme reduces 50 percent of the computational cost to achieve the similar modeling accuracy. This new optimized finite-difference scheme is particularly useful for large-scale 3D scalar-wave modeling and inversion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1995

Intermittent turbulence in solar wind from the south polar hole

The magnetic fields measured by the Ulysses spacecraft are used to study solar wind turbulence in the fast solar wind from the south polar hole. The spacecraft was at about 46 deg south latitude and 3.9 AU. For a magnetic field with a Gaussian distribution the power spectrum (second-order structure function) is sufficient to completely characterize the turbulence. However, the actual distribution is non-Gaussian so that the effects of intermittency must be taken into account. The observed spectral exponents include effects of intermittency and cannot be directly compared with the standard second-order spectral theories such as the Kolmogorov and Kraichnan theories. To permit a better comparison of the observations with the theoretical models, we study the structure characteristics of the data. We find the exponents of the second-order structure functions (power spectra) and the higher-order normalized structure functions for the components of the magnetic fields. We show that these sets of exponents can be approximately described by two basic numbers: the spectral exponent and the intermittency exponent. The intermittency exponent characterizes correlation properties of the energy cascade from large to small scales. Before comparing the observations to the theoretically expected values, a reduction must be made to the observed spectral exponent. The amount of the reduction depends on both the intermittency exponent and the model of the energy cascade assumed in the turbulence theory. We reduce the measured spectral indices according to a simple model for Alfven turbulence that is described here. We then compare our reduced spectral indices with second-order spectral theory. The reduced spectral indices for the period range of 1 min to about a half hour are remarkably constant and in good agreement with the value of 3/2. Thus our treatment is self-consistent. Our tentative conclusion is that the high-frequency turbulence appears to agree with the model of random-phased Alfven waves. This tentative conclusion must be tested by further theoretical and observational work.

Ruzmaikin, A. A.↗
DOE OSTI2025

When to vaccinate for seasonal influenza: check the peak forecast

Background Seasonal influenza infects 5-20% of people every year in the United States, resulting in hospitalizations, deaths, and adverse economic impacts. To mitigate these impacts, influenza vaccines are developed and distributed annually; however, growing evidence suggests that vaccine effectiveness (VE) wanes over the course of a flu season. Delaying influenza vaccination for older adults has attracted attention as a potential public health strategy. However, given the uncertainties in seasonal peak, vaccine effectiveness, and waning rates, postponing vaccination could also lead to increased morbidity, motivating an evaluation of a range of potential scenarios. The aim of this study was to investigate favorable age group-specific vaccination schedules that could lead to the greatest disease burden reduction. Methods We systematically investigated a broad range of vaccination start times for five age groups under six combinations of initial effectiveness and waning rates, based on influenza cases and vaccine uptake data from 10 influenza seasons. We defined the most favorable vaccination schedule as the one that resulted in the greatest reduction in disease burden. Results In scenarios with fast waning, all age groups benefit from delaying vaccination regardless of initial VE and peak timing. In scenarios with slower waning, results are mixed. For the ≥65 group, high initial VE and slow waning suggests that in early-peaking seasons, early vaccination most effectively reduces disease burden, while in late-peaking seasons delaying vaccination is most effective. For the ≥65 group in medium and low initial VE, and slow waning scenarios, delaying vaccination appears to prevent the greatest number of cases, regardless of whether the season peaks early or late. Conclusion The most favorable vaccination schedule is sensitive to changes in initial VE, waning rate, and peak timing. Given estimates of these quantities from statistical and immunological models and observations, our methods can inform vaccination recommendations in order to most effectively reduce the annual disease burden caused by seasonal influenza. Specifically, accurate peak timing forecasts for the upcoming season have the potential to guide decisions on when to vaccinate.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2000

Adding In-Plane Flexibility to the Equations of Motion of a Single Rotor Helicopter

This report describes a way to add the effects of main rotor blade flexibility in the in- plane or lead-lag direction to a large set of non-linear equations of motion for a single rotor helicopter with rigid blades(l). Differences between the frequency of the regressing lag mode predicted by the equations of (1) and that measured in flight (2) for a UH-60 helicopter indicate that some element is missing from the analytical model of (1) which assumes rigid blades. A previous study (3) noted a similar discrepancy for the CH-53 helicopter. Using a relatively simple analytical model in (3), compared to (1), it was shown that a mechanical lag damper increases significantly the coupling between the rigid lag mode and the first flexible mode. This increased coupling due to a powerful lag damper produces an increase in the lowest lag frequency when viewed in a frame rotating with the blade. Flight test measurements normally indicate the frequency of this mode in a non-rotating or fixed frame. This report presents the additions necessary to the full equations of motion, to include main rotor blade lag flexibility. Since these additions are made to a very complex nonlinear dynamic model, in order to provide physical insight, a discussion of the results obtained from a simplified set of equations of motion is included. The reduced model illustrates the physics involved in the coupling and should indicate trends in the full model.

Curtiss, H. C., Jr.↗
NASA NTRS2000

RFI and SCRIMP Model Development and Verification

Vacuum-Assisted Resin Transfer Molding (VARTM) processes are becoming promising technologies in the manufacturing of primary composite structures in the aircraft industry as well as infrastructure. A great deal of work still needs to be done on efforts to reduce the costly trial-and-error methods of VARTM processing that are currently in practice today. A computer simulation model of the VARTM process would provide a cost-effective tool in the manufacturing of composites utilizing this technique. Therefore, the objective of this research was to modify an existing three-dimensional, Resin Film Infusion (RFI)/Resin Transfer Molding (RTM) model to include VARTM simulation capabilities and to verify this model with the fabrication of aircraft structural composites. An additional objective was to use the VARTM model as a process analysis tool, where this tool would enable the user to configure the best process for manufacturing quality composites. Experimental verification of the model was performed by processing several flat composite panels. The parameters verified included flow front patterns and infiltration times. The flow front patterns were determined to be qualitatively accurate, while the simulated infiltration times over predicted experimental times by 8 to 10%. Capillary and gravitational forces were incorporated into the existing RFI/RTM model in order to simulate VARTM processing physics more accurately. The theoretical capillary pressure showed the capability to reduce the simulated infiltration times by as great as 6%. The gravity, on the other hand, was found to be negligible for all cases. Finally, the VARTM model was used as a process analysis tool. This enabled the user to determine such important process constraints as the location and type of injection ports and the permeability and location of the high-permeable media. A process for a three-stiffener composite panel was proposed. This configuration evolved from the variation of the process constraints in the modeling of several different composite panels. The configuration was proposed by considering such factors as: infiltration time, the number of vacuum ports, and possible areas of void entrapment.

Loos, Alfred C.↗
NASA NTRS1973

Resonant oscillations of intermediate frequency in a stratified atmosphere.

A class of solutions to a model of forced oscillations in a rotating stratified atmospheric layer is derived and analyzed. The basic model is found to reduce to a boundary value problem with a second-order linear partial differential equation of the hyperbolic type for this range of forcing frequencies. The forced solutions are shown to exhibit resonances with the normal modes of oscillation of the layer. The characteristics of the resonant modes are analyzed in terms of mean tropospheric values of temperature, temperature lapse, wind speed, horizontal and vertical wind shears, latitude, and the frequency and horizontal wavelength of the forcing mechanism. These solutions are compared with solutions to the model for a different (subinertial) range of forcing frequencies. This comparison leads to an elliptic boundary value problem. The solutions in that case do not exhibit the same type of resonance and generally decay away from the region of forcing.

Thompson, O. E.↗
DOE OSTI2023

Mshock HED Council Request FY24 [Slides]

ICF and HEDP experiments experience complicated loadings including successive shocks, but modeling and analytic work has mostly engaged the simpler shock or shock-reshock case. This is because the co propagating case is difficult to achieve with conventional (non-HED) drivers. Successive shocks will be challenging to model in BHR or modal model interpretations: 1) Ex: Consider the case where long wavelength modes re-invert coherently, but short wavelength nonlinear modes are spun up turbulently. 2) Necessitates more advanced diagnostics and analysis than simple “mixing layers.” High resolution spectral information (concomitant with Ω-EP campaigns) and higher-order moment analyses will be necessary to stress and validate the reduced BHR-type models.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Scalar QED as a toy model for higher-order effects in classical gravitational scattering

Quantum Electrodynamics (QED) serves as a useful toy model for classical observables in gravitational two-body systems with reduced complexity due to the linearity of QED. We investigate scattering observables in scalar QED at the sixth order in the charges (two-loop order) in a classical regime analogous to the post-Minkowskian expansion in General Relativity. We employ modern scattering amplitude tools and extract classical observables by both eikonal methods and the formalism of Kosower, Maybee, and O’Connell (KMOC). In addition, we provide a simplified approach to extracting the radial action beyond the conservative sector.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1985

On-orbit systems identification of flexible spacecraft

The process of on-orbit systems identification of flexible spacecraft is examined in terms of the difficulties that are expected because of the potentially unmanageable number of unknown model parameters due to the very high order system models involved. A Jordon block canonical form and global model parameters are used to reduce the number of unknown parameters to manageable numbers. A Bayesian approach is discussed which enables the merging of theoretical models with ground or on-orbit test results by using Unconditional Maximum Likelihood Parameter Estimation. A Modified Newton-Raphson technique is proposed for determining the model parameter estimates and an expected fit error criterion is recommended for the determination of the model structure and order reduction.

Taylor, L. W., Jr.↗
DOE OSTI2020

Distribution System Segmented Model Simplification with Independent Dynamically Changing End-Use Loads

With more emerging technologies in power distribution systems, to support quasi-static time-series (QSTS) simulation over short time frames, accurately modeling the dynamic behavior of end-use load under various conditions is im- portant. However, time and cost can be a burden for the QSTS simulation of thousands of houses representing dynamic load behaviors on large distribution feeders. There are two main contributions of this paper in order to reduce the QSTS simulation complexity. First, the paper presents a method to stochasti- cally model the dynamic behavior of end-use load for QSTS simulations. Second, using the stochastically modeled load behavior as an input to a distribution sys- tem, and building on our previous QSTS simulation simplification/acceleration method that relies on uniform load behavior, this paper presents a segment substitution approach for the simplification of distribution system models with independent dynamically changing end-use loads. Our numerical results show that compared to GridLAB-D, a software tool that simulates end-use load be- havior through physical models, the proposed stochastic modeling approach achieved an 89.3% simulation time reduction for end-use loads, while accurately capturing the temporally changing total load behavior. Moreover, the QSTS simulation of a large test distribution system through OpenDSS showed that the proposed simplification method decreased the feeder simulation time by a factor of three at less than 0.3% voltage error. Compared to our previous method with independent dynamic loads, the voltage error is reduced by a factor of four.

Power distribution system modeling, stochastic pro↗
NASA NTRS2002

Experimental Feedback Control of Flow Induced Cavity Tones

An experimental study of the application of discrete-time, linear quadratic control design methods to the cavity tone problem is described. State space models of the dynamics from a synthetic jet actuator at the leading edge of the cavity to two pressure sensors in the cavity were computed from experimental data. Variations in model order, control order, control bandwidth, and properties of a Kalman state estimator were studied. Feedback control reduced the levels of multiple cavity tones at Mach 0.275, 0.35, and 0.45. Closed loop performance was often limited by excitation of sidebands of cavity tones, and creation of new tones in the spectrum. State space models were useful for explaining some of these limitations, but were not able to account for non-linear dynamics, such as interactions between tones at different frequencies.

Cabell, Randolph H.↗
DOE OSTI2021

Machine learning based rate optimization under geologic uncertainty

We propose a novel approach for rate optimization during a waterflood under geologic uncertainty in reservoir properties such as permeability and porosity. The traditional approach typically involves several runs of the forward simulator. This may not scale well when the optimization is to be performed at the full field-level and over multiple geologic realizations. A machine-learning (ML) based approach which is quick and scalable for rate optimization over multiple geologic realizations is proposed instead. The training data for the model is generated by running the forward simulator with randomly assigned well rates using multiple geologic realizations. A reduced order representation of the permeability heterogeneity in each of the realizations is derived using a grid connectivity transformation (GCT). This step involves finding basis functions corresponding to the different modal frequencies of the grid connectivity represented by the grid Laplacian. The projection of the heterogeneous property field along these basis functions gives the basis coefficients that form the reduced order representation. Subsequently, for each training datapoint, streamlines are traced and the minimum time of flight (TOF) representing the tracer breakthrough time at each producer is recorded. The basis coefficients and well rates are fed to a machine learning model as input and the minimum TOF at the producers forms the output of the model. This trained model can then be used along with an optimizer for computing the optimal injection rates to maximize the injection sweep efficiency. This corresponds to minimizing the variance in the minimum TOF within each well group. Different architectures of neural network are tested using 5-fold cross validation to decide the best ML model to compute the streamline time of flight. The trained model is used to perform well rate optimization over multiple realizations of geology by using a risk tolerance penalty. The optimal well rates thus obtained are compared with two cases: a) equal well rates assigned to all injectors and producers and b) well rates obtained by optimizing over a single realization without considering the uncertainty in geology. The optimal well rates are seen to offer better oil recovery and sweep efficiency than both cases.

02 PETROLEUM↗
DOE OSTI2022

Transient behavior of drift and ionization in atmospheric pressure nitrogen discharge

The fluid models are frequently used to describe a non-thermal plasma such as a streamer discharge. The required electron transport data and rate coefficients for the fluid model are parametrized using the local field approximation (LFA) in first order models and the local-mean-energy approximation (LMEA) in second order models. We performed Monte Carlo simulations in Nitrogen gas with step changes in the E/N (reduced electric field) to study the behavior of the transport properties in the transient phase. During the transient phase of the simulation, we extract the instantaneous electron mean energy, which is different from the steady state mean electron energy, and the corresponding transport parameters and rate coefficients. Our results indicate that the mean electron energy is not a suitable parameter for mobility/drift of electrons due to big difference in momentum relaxation and energy relaxation. Furthermore, the high energy threshold rates such as ionization show a strong correlation to mean electron energy. In second order models where the energy-balance equation is solved, we suggest that it would rather be appropriate to use the local electric field to find electron drift velocity in gases such as Nitrogen and the local mean electron energy to determine the ionization and excitation rates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1996

Modeling Orifice Pulse Tube Coolers

We have developed a calculational model that treats all the components of an orifice pulse tube cooler. We base our analysis on 1-dimensional thermodynamic equations for the regenerator and we assume that all mass flows, pressure oscillations and temperature oscillations are small and sinusoidal. Non-linear pressure drop effects are included in the regenerator to account for finite pressure amplitude effects. The resulting mass flows and pressures are matched at the boundaries with the other components of the cooler: compressor, aftercooler, cold heat exchanger, pulse tube, hot heat exchanger, orifice and reservoir. The results of the calculation are oscillating pressures, mass flows and enthalpy flows in the main components of the cooler. By comparing with the calculations of other available models, we show that our model is very similar to REGEN 3 from NIST and DeltaE from Los Alamos National Lab. Our model is much easier to use than other available models because of its simple graphical interface and the fact that no guesses are required for the operating pressures or mass flows. In addition, the model only requires a few minutes of running time allowing many parameters to be optimized in a reasonable time. A version of the model is available for use over the World Wide Web at http://irtek.arc.nasa.gov. Future enhancements include adding a bypass orifice and including second order terms in steady mass streaming and steady heat transfer. A two-dimensional anelastic approximation of the fluid equations will be used as the basis for the latter analysis. Preliminary results are given in dimensionless numbers appropriate for oscillating compressible flows. The model shows how transverse heat transfer reduces enthalpy flow, particularly for small pulse tubes. The model also clearly shows mass recirculation in the open tube on the order of the tube length. They result from the higher order Reynolds stresses. An interesting result of the linearized approach is that the steady mass streaming does not affect the enthalpy flow at second order. The major effect of recirculating mass streaming is to increase transverse temperature gradients, which leads to higher entropy production and reduced efficiency.

Roach, Kittel P.↗
NASA NTRS2007

Numerical Modeling of Propellant Boiloff in Cryogenic Storage Tank

This Technical Memorandum (TM) describes the thermal modeling effort undertaken at Marshall Space Flight Center to support the Cryogenic Test Laboratory at Kennedy Space Center (KSC) for a study of insulation materials for cryogenic tanks in order to reduce propellant boiloff during long-term storage. The Generalized Fluid System Simulation program has been used to model boiloff in 1,000-L demonstration tanks built for testing the thermal performance of glass bubbles and perlite insulation. Numerical predictions of boiloff rate and ullage temperature have been compared with the measured data from the testing of demonstration tanks. A satisfactory comparison between measured and predicted data has been observed for both liquid nitrogen and hydrogen tests. Based on the experience gained with the modeling of the demonstration tanks, a numerical model of the liquid hydrogen storage tank at launch complex 39 at KSC was built. The predicted boiloff rate of hydrogen has been found to be in good agreement with observed field data. This TM describes three different models that have been developed during this period of study (March 2005 to June 2006), comparisons with test data, and results of parametric studies.

Majumdar, A. K.↗
DOE OSTI2013

Reduction of discretization error for ray tracing of MOC through a correction on collision probabilities

A new correction model for ray tracing of the method of characteristics is proposed in order to reduce discretization error. As the ray tracing parameters such as azimuthal angle division, polar angle division and ray separation are considered in this study. In the method of characteristics, region average scalar fluxes can be implicitly expressed by collision probabilities, although these collision probabilities are not directly treated in the ordinary calculation scheme. From this viewpoint, difference between a coarse ray tracing condition and a detailed one can be interpreted as the difference in the estimation of collision probabilities. In other words, the discretization error for ray tracing can be recognized as a consequence of inaccurate collision probabilities caused by coarse ray tracing. This discussion suggests that accurate region average scalar flux can be obtained through an appropriate correction on collision probabilities. In this paper, a correction model on collision probabilities is theoretically derived based on the neutron balance equation, and its validity is confirmed through typical single assembly calculations. The effectiveness of the present correction method is also discussed in this paper. It is confirmed that discretization error for ray tracing can be significantly reduced by the present correction method in a multi-assembly calculation, though the correction factor is estimated in single assembly geometry. (authors)

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Expediting model-based optoacoustic reconstructions with tomographic symmetries

Purpose: Image quantification in optoacoustic tomography implies the use of accurate forward models of excitation, propagation, and detection of optoacoustic signals while inversions with high spatial resolution usually involve very large matrices, leading to unreasonably long computation times. The development of fast and memory efficient model-based approaches represents then an important challenge to advance on the quantitative and dynamic imaging capabilities of tomographic optoacoustic imaging. Methods: Herein, a method for simplification and acceleration of model-based inversions, relying on inherent symmetries present in common tomographic acquisition geometries, has been introduced. The method is showcased for the case of cylindrical symmetries by using polar image discretization of the time-domain optoacoustic forward model combined with efficient storage and inversion strategies. Results: The suggested methodology is shown to render fast and accurate model-based inversions in both numerical simulations andpost mortem small animal experiments. In case of a full-view detection scheme, the memory requirements are reduced by one order of magnitude while high-resolution reconstructions are achieved at video rate. Conclusions: By considering the rotational symmetry present in many tomographic optoacoustic imaging systems, the proposed methodology allows exploiting the advantages of model-based algorithms with feasible computational requirements and fast reconstruction times, so that its convenience and general applicability in optoacoustic imaging systems with tomographic symmetries is anticipated.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2013

A Diagnoser Algorithm for Anomaly Detection in DEDS under Partial Unreliable Observations: Characterization and Inclusion in Sensor Configuration Optimizaton

Complex engineering systems have to be carefully monitored to meet demanding performance requirements, including detecting anomalies in their operations. There are two major monitoring challenges for these systems. The first challenge is that information collected from the monitored system is often partial and/or unreliable, in the sense that some occurred events may not be reported and/or may be reported incorrectly (e.g., reported as another event). The second is that anomalies often consist of sequences of event patterns separated in space and time. This paper introduces and analyzes a diagnoser algorithm that meets these challenges for detecting and counting occurrences of anomalies in engineering systems. The proposed diagnoser algorithm assumes that models are available for characterizing plant operations (via stochastic automata) and sensors (via probabilistic mappings) used for reporting partial and unreliable information. Methods for analyzing the effects of model uncertainties on the diagnoser performance are also discussed. In order to select configurations that reduce sensor costs, while satisfying diagnoser performance requirements, a sensor configuration selection algorithm developed in previous work is then extended for the proposed diagnoser algorithm. The proposed algorithms and methods are then applied to a multi-unit-operation system, which is derived from an actual facility application. Results show that the proposed diagnoser algorithm is able to detect and count occurrences of anomalies accurately and that its performance is robust to model uncertainties. Furthermore, the sensor configuration selection algorithm is able to suggest optimal sensor configurations with significantly reduced costs, while still yielding acceptable performance for counting the occurrences of anomalies.

97 MATHEMATICS AND COMPUTING↗