Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reduce 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,225 records · Page 68

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.↗

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↗

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↗

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.↗

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.↗

Clock model interpolation and symmetry breaking in O(2) models

The $q$-state clock model is a classical spin model that corresponds to the Ising model when $q=2$ and to the $XY$ model when $q\to\infty$. The integer-$q$ clock model has been studied extensively and has been shown to have a single phase transition when $q=2$,$3$,$4$ and two phase transitions when $q>4$.We define an extended $q$-state clock model that reduces to the ordinary $q$-state clock model when $q$ is an integer and otherwise is a continuous interpolation of the clock model to noninteger $q$. We investigate this class of clock models in 2D using Monte Carlo (MC) and tensor renormalization group (TRG) methods, and we find that the model with noninteger $q$ has a crossover and a second-order phase transition. We also define an extended-$O(2)$ model (with a parameter $\gamma$) that reduces to the $XY$ model when $\gamma=0$ and to the extended $q$-state clock model when $\gamma\to\infty$, and we begin to outline the phase diagram of this model. These models with noninteger $q$ serve as a testbed to study symmetry breaking in situations corresponding to quantum simulators where experimental parameters can be tuned continuously.

Hostetler, Leon↗

The MicroBooNE Single-Photon Low-Energy Excess Search (Public Note 1087)

MicroBooNE is a short baseline neutrino experiment at Fermilab designed to address the low energy excess observed by the MiniBooNE experiment. This note describes and presents preliminary results for the MicroBooNE analysis developed to address this excess as a single photon plus one or zero protons in the final state. The analysis assumes neutrino neutral current Δ resonance production followed by Δ radiative decay on argon (NC Δ → Nγ ) as the "signal model"; event reconstruction and selection have been developed and optimized in order to maximize efficiency and reduce cosmogenic and other beam-related backgrounds to the NC Δ → Nγ signal. We present the analysis methodology and validation checks performed on limited-statistics open data sets, corresponding to 5 x10 19 protons on target (POT), following a blind analysis, as well as the projected sensitivities for testing the Standard Model (SM) predicted rate for the NC Δ → Nγ process and for testing the interpretation of the previously observed MiniBooNE low energy excess as NC Δ → Nγ events, using the full anticipated MicroBooNE data set of 12.25 x 10 20 POT.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The MicroBooNE Single-Photon Low-Energy Excess Search

MicroBooNE is a short baseline neutrino experiment at Fermilab designed to address the low energy excess observed by the MiniBooNE experiment. This note describes and presents preliminary results for the MicroBooNE analysis developed to address this excess as a single photon plus one or zero protons in the final state. The analysis assumes neutrino neutral current Δ resonance production followed by Δ radiative decay on argon (NC Δ → $N_γ$) as the “signal model”; event reconstruction and selection have been developed and optimized in order to maximize efficiency and reduce cosmogenic and other beam-related backgrounds to the NC Δ → $N_γ$ signal. We present the analysis methodology and validation checks performed on limited-statistics open data sets, corresponding to 5×10 19 protons on target (POT), following a blind analysis, as well as the projected sensitivities for testing the Standard Model (SM) predicted rate for the NC Δ → $N_γ$ process and for testing the interpretation of the previously observed MiniBooNE low energy excess as NC Δ → $N_γ$ events, using the full anticipated MicroBooNE data set of 12.25×10 20 POT.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Explicit modeling and concurrent processing in the simulation of multibody dynamic systems

The objective is to present the activities at TRW in developing the capability to simulate the behavior of large flexible multibody space structures. The features of the simulation tools are: (1) to accommodate all rigid/flexible body degrees-of-freedom which incorporate the control system models and external forces, (2) to provide the flexibility to incorporate engineering-defined models and to retain parameters of significance to the engineer, (3) to reduce the computation cost by one order of magnitude (two orders of magnitude compared to a CRAY 1S), and (4) to keep it versatile so that radical variations in anticipated space structures can be accommodated. The current computer tools to simulate multibody systems appear not only to be very costly and time consuming, but also do not produce the desired fidelity of the mathematical models. In summary, a multibody simulation tool will be developed in the near future which will allow solution of the dynamics and controls of the deployment of the LDR backup structure, or the problem associated with the robotic assembly of the structure. The tools will allow the engineer to define the modeling technique and solve problems in less time and at reduced cost.

Gluck, R.↗

Reducing model uncertainty effects in flexible manipulators through the addition of passive damping

An important issue in the control of practical systems is the effect of model uncertainty on closed loop performance. This is of particular concern when flexible structures are to be controlled, due to the fact that states associated with higher frequency vibration modes are truncated in order to make the control problem tractable. Digital simulations of a single-link manipulator system are employed to demonstrate that passive damping added to the flexible member reduces adverse effects associated with model uncertainty. A controller was designed based on a model including only one flexible mode. This controller was applied to larger order systems to evaluate the effects of modal truncation. Simulations using a Linear Quadratic Regulator (LQR) design assuming full state feedback illustrate the effect of control spillover. Simulations of a system using output feedback illustrate the destabilizing effect of observation spillover. The simulations reveal that the system with passive damping is less susceptible to these effects than the untreated case.

Alberts, T. E.↗

Machine Learning-Based Upscaling of Finite-Size Molecular Dynamics Diffusion Simulations for Binary Fluids

Molecular diffusion coefficients calculated using molecular dynamics (MD) simulations suffer from finite-size (i.e., finite box size and finite particle number) effects. Results from finite-sized MD simulations can be upscaled to infinite simulation size by applying a correction factor. For self-diffusion of single-component fluids, this correction has been well-studied by many researchers including Yeh and Hummer (YH); for binary fluid mixtures, a modified YH correction was recently proposed for correcting MD-predicted Maxwell–Stephan (MS) diffusion rates. In this study we use both empirical and machine learning methods to identify improvements to the finite-size correction factors for both self-diffusion and MS diffusion of binary Lennard-Jones (LJ) fluid mixtures. Using artificial neural networks (ANNs), the error in the corrected LJ fluid diffusion is reduced by an order of magnitude versus existing YH corrections, and the ANN models perform well for mixtures with large dissimilarities in size and interaction energies where the YH correction proves insufficient.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Stacked-mosaic amplifier and diode delivery concept for kJ-class inertial fusion energy laser driver

We report on a laser amplifier architecture designed as a modular element for scaling to multi-megajoule laser facilities intended for inertial fusion energy (IFE) power plants. This kilojoule-class module features a stacked mosaic of gain media integrated with a diode delivery system for efficient optical pumping of the mosaic-structured gain medium, enabling high-repetition-rate operation at high wall plug efficiencies, necessary for an IFE driver. Details of a diode delivery system capable of pumping the mosaic architecture are presented. Comprehensive numerical modeling demonstrates that the stacked-mosaic approach reduces transverse gain by five orders of magnitude compared to a conventional full-aperture Yb:YAG based amplifier, substantially suppressing transverse amplified spontaneous emission (TASE) and enabling enhanced longitudinal energy extraction. Detailed analysis of thermal management and wavefront distortion in a 4 × 4 mosaic array indicates that temperature gradients and thermally induced aberrations are effectively controlled using gas cooling and commercially available phase-plate and deformable mirror technologies. We further discuss the applicability of a stacked-mosaic architecture to direct-drive IFE schemes, where broad spectral bandwidth is critical, and its compatibility with frequency conversion modules for up-conversion to blue wavelengths. Finally, an example point design for a 10 kJ, 10 Hz Yb:YAG module operating at 175 K with wall-plug efficiency exceeding 10 % is presented, underscoring the feasibility of this approach for next-generation high-energy lasers for IFE drivers. The results establish the stacked-mosaic amplifier as a scalable, robust platform not only for IFE but also for a broad range of advanced scientific and industrial laser applications.

Lasers↗

Altimeter error sources at the 10-cm performance level

Error sources affecting the calibration and operational use of a 10 cm altimeter are examined to determine the magnitudes of current errors and the investigations necessary to reduce them to acceptable bounds. Errors considered include those affecting operational data pre-processing, and those affecting altitude bias determination, with error budgets developed for both. The most significant error sources affecting pre-processing are bias calibration, propagation corrections for the ionosphere, and measurement noise. No ionospheric models are currently validated at the required 10-25% accuracy level. The optimum smoothing to reduce the effects of measurement noise is investigated and found to be on the order of one second, based on the TASC model of geoid undulations. The 10 cm calibrations are found to be feasible only through the use of altimeter passes that are very high elevation for a tracking station which tracks very close to the time of altimeter track, such as a high elevation pass across the island of Bermuda. By far the largest error source, based on the current state-of-the-art, is the location of the island tracking station relative to mean sea level in the surrounding ocean areas.

Martin, C. F.↗

Effects of moisture, residual thermal curing stresses and mechanical load on the damage development in quasi-isotropic laminates

Classical laminate theory and a finite element model were used to predict stress states prior to the first formation of damage in laminates fabricated from T/300/5208. Crack patterns characteristic of the laminate in a wet or dry condition were also predicted using a shear lag model. Development of edge damage was recorded and observed during the test by transferring an image of the damage from the edge surface on to a thin acetate sheet such that the damage imprinted could be immediately viewed on a microfiche card reader. Moisture was shown to significantly alter the interior and edge dry stress states due to swelling and a reduction of elastic properties and to reduce the transverse strength in 90 deg plies. A model was developed in order to predict changes in first ply failure laminate loads due to differences in stacking sequence together with a wet or dry environmental condition.

Kriz, R. D.↗

The delineation and interpretation of the Earth's gravity field

The observed changes in velocity with time are reduced relative to the well-determined low degree and order GEM field model and accelerations are found by analytical differentiation of the range rates. This new map is essentially identical to the first map and we have produced a composite map by combining all 90 passes of SST data. The resolution of the map is at worst about 5 deg and much better in most places. A comparison of this map with conventional GEM models shows very good agreement. A reduction of the SEASAT altimeter data has also been carried out for an additional comparison. Although the SEASAT geoid contains much more high frequency information, it agrees very well with both the SST and GEM fields. The maps are dominated (especially in the east) by a pattern of roughly east-west anomalies with a transverse wavelength of about 2000 km. A further comparison with regional bathymetric data shows a remarkably close correlation with plate age.

Marsh, B. D.↗

Controller reduction by preserving impulse response energy

A model order reduction algorithm based on a Krylov recurrence formulation is developed to reduce order of controllers. The reduced-order controller is obtained by projecting the full-order LQG controller onto a Krylov subspace in which either the controllability or the observability grammian is equal to the identity matrix. The reduced-order controller preserves the impulse response energy of the full-order controller and has a parameter-matching property. Two numerical examples drawn from other controller reduction literature are used to illustrate the efficacy of the proposed reduction algorithm.

Craig, Roy R., Jr.↗

Research in robust control for hypersonic aircraft

The research during the third reporting period focused on fixed order robust control design for hypersonic vehicles. A new technique was developed to synthesize fixed order H(sub infinity) controllers. A controller canonical form is imposed on the compensator structure and a homotopy algorithm is employed to perform the controller design. Various reduced order controllers are designed for a simplified version of the hypersonic vehicle model used in our previous studies to demonstrate the capabilities of the code. However, further work is needed to investigate the issue of numerical ill-conditioning for large order systems and to make the numerical approach more reliable.

Calise, A. J.↗

Drive Rig Mufflers for Model Scale Engine Acoustic Testing

Testing of air breathing propulsion systems in the 9x15 foot wind tunnel at NASA Glenn Research Center depends on compressed air turbines for power. The drive rig turbines exhaust directly to the wind tunnel test section, and have been found to produce significant unwanted noise that reduces the quality of the acoustic measurements of the model being tested. In order to mitigate this acoustic contamination, a muffler can be attached downstream of the drive rig turbine. The modern engine designs currently being tested produce much less noise than traditional engines, and consequently a lower noise floor is required of the facility. An acoustic test of a muffler designed to mitigate this extraneous noise is presented, and a noise reduction of 8 dB between 700 Hz and 20 kHz was documented, significantly improving the quality of acoustic measurements in the facility.

Stephens, David↗