Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic 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.

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

At least 901 records · Page 50

Self-Aligned Focusing Schlieren at the 0.3-M Transonic Cryogenic Tunnel and the National Transonic Facility

The implementation of a self-aligned focusing schlieren (SAFS) system at two cryogenic wind tunnels at NASA Langley Research Center is discussed. Risk-reduction testing of the SAFS system was first performed at the 0.3-M Transonic Cryogenic Tunnel to evaluate the system’s operation in a small-scale characteristic cryogenic facility. Testing was conducted with three models: a three-quarter span 25.4-mm-diameter cylinder, a semi-span 65A006 tapered unswept airfoil, and a full-span SC(3)-0712 airfoil. Testing with the cylinder and semi-span airfoil revealed a highly dynamic shock environment, whereas the shock on the full-span airfoil was stationary, solidifying the usage of this model for a pre-/post-shock particle tracking velocimetry measurement. Temperature-induced polarization-altering window stresses were encountered during low-temperature testing and were mitigated using a “non-ideal” quartz/quartz Rochon prism that had largely been neglected since the SAFS system’s first introduction in favor of the more favorable “ideal” glass/quartz Rochon prism. The size of the SAFS system was then decreased in order to fit inside an environmentally controlled camera can enclosure at the National Transonic Facility (NTF) for testing of a sting-mounted aircraft model. The SAFS system was demonstrated to be effective at filtering out the large density gradient flow in the0.3-M plenum, and the thick, high density turbulent boundary layers on the wind tunnel walls at the NTF. Results of the testing campaigns and improvements to future systems are discussed.

Joshua M Weisberger↗

Model reduction by weighted Component Cost Analysis

Component Cost Analysis considers any given system driven by a white noise process as an interconnection of different components, and assigns a metric called 'component cost' to each component. These component costs measure the contribution of each component to a predefined quadratic cost function. A reduced-order model of the given system may be obtained by deleting those components that have the smallest component costs. The theory of Component Cost Analysis is extended to include finite-bandwidth colored noises. The results also apply when actuators have dynamics of their own. Closed-form analytical expressions of component costs are also derived for a mechanical system described by its modal data. This is very useful to compute the modal costs of very high order systems. A numerical example for MINIMAST system is presented.

Kim, Jae H.↗

Validation of the ERO2.0 code using W7-X and JET experiments and predictions for ITER operation

Abstract The paper provides an overview of recent modelling of global material erosion and deposition in the fusion devices Wendelstein 7-X (W7-X), JET and ITER using the Monte-Carlo code ERO2.0. For validating the modelling tool in a three-dimensional environment, W7-X simulations are performed to describe carbon erosion from the graphite test divertor units, which were equipped in operational phase OP 1.2 and analysed post-mortem. Synthetic spectroscopy of carbon line emission is compared with experimental results from the divertor spectrometer measurement system, showing a good agreement in the e-folding lengths in the radial intensity profiles of carbon. In the case of metallic wall materials, earlier modelling of the Be/W environment in JET and ITER is revisited and extended with an updated set of sputtering and reflection data, as well as including the mixing model for describing the Be/W dynamics in the divertor. Motivated by recent H/D/T isotope experiments in JET, limited and diverted configuration pulses are modelled, showing the expected trend of both Be and W erosion increasing with isotope mass. For the JET diverted configuration pulses, it is shown that Be migrates predominantly to the upper part of the inner divertor where it initially leads to strong W erosion. With longer exposure time, the growth of a Be deposited layer leads to a reduction of W erosion in that region. A similar trend is observed in simulations of the ITER baseline Q = 10 scenario, however with a more symmetric Be migration pattern leading to deposition also on the outer divertor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Breakups are complicated: an efficient representation of collisional breakup in the superdroplet method

Abstract. A key constraint of particle-based methods for modeling cloud microphysics is the conservation of total particle number, which is required for computational tractability. The process of collisional breakup poses a particular challenge to this framework, as breakup events often produce many droplet fragments of varying sizes, which would require creating new particles in the system. This work introduces a representation of collisional breakup in the so-called “superdroplet” method which conserves the total number of superdroplets in the system. This representation extends an existing stochastic collisional-coalescence scheme and samples from a fragment size distribution in an additional Monte Carlo step. This method is demonstrated in a set of idealized box model and single-column warm-rain simulations. We further discuss the effects of the breakup dynamic and fragment size distribution on the particle size distribution, hydrometeor population, and microphysical process rates. Box model experiments serve to characterize the impacts of properties such as coalescence efficiency and fragmentation function on the relative roles of collisional breakup and coalescence. The results demonstrate that this representation of collisional breakup can produce a stationary particle size distribution, in which breakup and coalescence rates are approximately equal, and that it recovers expected behavior such as a reduction in precipitate-sized particles in the column model. The breakup algorithm presented here contributes to an open-source pythonic implementation of the superdroplet method, PySDM, which will facilitate future research using particle-based microphysics.

de Jong, Emily (ORCID:0000000253104554)↗

Rotordynamic Characteristics of the HPOTP (High Pressure Oxygen Turbopump) of the SSME (Space Shuttle Main Engine)

Rotational stability of turbopump components in the space shuttle main engine was studied via analysis of component and structural dynamic models. Subsynchronous vibration caused unacceptable migration of the rotor/housing unit with unequal load sharing of the synchronous bearings that resulted in the failure of the High Pressure Oxygen Turbopump. Linear analysis shows that a shrouded inducer eliminates the second critical speed and the stability problem, a stiffened rotor improves the rotordynamic characteristics of the turbopump, and installing damper boost/impeller seals reduces bearing loads. Nonlinear analysis shows that by increasing the "dead band' clearances, a marked reduction in peak bearing loads occurs.

Childs, D. W.↗

Damage-Mitigating Control of Space Propulsion Systems for High Performance and Extended Life

A major goal in the control of complex mechanical system such as spacecraft rocket engine's advanced aircraft, and power plants is to achieve high performance with increased reliability, component durability, and maintainability. The current practice of decision and control systems synthesis focuses on improving performance and diagnostic capabilities under constraints that often do not adequately represent the materials degradation. In view of the high performance requirements of the system and availability of improved materials, the lack of appropriate knowledge about the properties of these materials will lead to either less than achievable performance due to overly conservative design, or over-straining of the structure leading to unexpected failures and drastic reduction of the service life. The key idea in this report is that a significant improvement in service life could be achieved by a small reduction in the system dynamic performance. The major task is to characterize the damage generation process, and then utilize this information in a mathematical form to synthesize a control law that would meet the system requirements and simultaneously satisfy the constraints that are imposed by the material and structural properties of the critical components. The concept of damage mitigation is introduced for control of mechanical systems to achieve high performance with a prolonged life span. A model of fatigue damage dynamics is formulated in the continuous-time setting, instead of a cycle-based representation, for direct application to control systems synthesis. An optimal control policy is then formulated via nonlinear programming under specified constraints of the damage rate and accumulated damage. The results of simulation experiments for the transient upthrust of a bipropellant rocket engine are presented to demonstrate efficacy of the damage-mitigating control concept.

Ray, Asok↗

Insights on the Effects of Magnetic Forces on the Efficiency of Vibration Energy Harvesting Absorbers in Controlling Dynamical Systems

This study investigates the effects of magnetic constraints on a piezoelectric energy harvesting absorber while simultaneously controlling a primary structure and harnessing energy. An accurate forcing representation of the magnetic force is investigated and developed. A reduced-order model is derived using the Euler–Lagrange principle, and the impact of the magnetic force is evaluated on the absorber’s static position and coupled natural frequency of the energy harvesting absorber and the coupled primary absorber system. The results show that attractive magnet configurations cannot improve the system substantially before pull-in occurs. A rigorous eigenvalue problem analysis is performed on the absorber’s substrate thickness and tip mass to effectively design an energy harvesting absorber for multiple initial gap sizes for the repulsive configurations. Then, the effects of the forcing amplitude on the primary structure absorber are studied and characterized by determining an effective design of the system for a simultaneous reduction in the primary structure’s motion and improvement in the harvester’s efficiency.

42 ENGINEERING↗

Large global variations in the carbon dioxide removal potential of seaweed farming due to biophysical constraints

Abstract Estimates suggest that over 4 gigatons per year of carbon dioxide (Gt-CO 2 year −1 ) be removed from the atmosphere by 2050 to meet international climate goals. One strategy for carbon dioxide removal is seaweed farming; however its global potential remains highly uncertain. Here, we apply a dynamic seaweed growth model that includes growth-limiting mechanisms, such as nitrate supply, to estimate the global potential yield of four types of seaweed. We estimate that harvesting 1 Gt year −1 of seaweed carbon would require farming over 1 million km 2 of the most productive exclusive economic zones, located in the equatorial Pacific; the cultivation area would need to be tripled to attain an additional 1 Gt year −1 of harvested carbon, indicating dramatic reductions in carbon harvest efficiency beyond the most productive waters. Improving the accuracy of annual harvest yield estimates requires better understanding of biophysical constraints such as seaweed loss rates (e.g., infestation, disease, grazing, wave erosion).

54 ENVIRONMENTAL SCIENCES↗

Tenth NASTRAN User's Colloquium

The development of the NASTRAN computer program, a general purpose finite element computer code for structural analysis, was discussed. The application and development of NASTRAN is presented in the following topics: improvements and enhancements; developments of pre and postprocessors; interactive review system; the use of harmonic expansions in magnetic field problems; improving a dynamic model with test data using Linwood; solution of axisymmetric fluid structure interaction problems; large displacements and stability analysis of nonlinear propeller structures; prediction of bead area contact load at the tire wheel interface; elastic plastic analysis of an overloaded breech ring; finite element solution of torsion and other 2-D Poisson equations; new capability for elastic aircraft airloads; usage of substructuring analysis in the get away special program; solving symmetric structures with nonsymmetric loads; evaluation and reduction of errors induced by Guyan transformation.

Source record↗

Optical imaging spectroscopy

During the recent solar maximum the combination of imaging and spectroscopy in the visible part of the spectrum became a powerful tool for observational study of flares primarily because of the development of two-dimensional charge-coupled-device (CCD) arrays. In combination with appropriate new operational methods, this has led to the ability to observe, for the first time, the preflare and impulsive-phase physical processes associated with spatially resolved features of flare loops. As a result of concurrent theoretical developments, modeling progressed from an empirical to a physical level. This made it possible to interpret imaging spectra in terms of coronal pressure and heat flux, particle beam heating, chromospheric evaporation, and explosive chromospheric dynamics at the footpoints of flare loops. There is clear potential for further advances in the near future, taking advantage of improvements in digital recording speed (approx. 10-fold), number of photosensitive elements per array (approx. 10-fold), real-time data pre-reduction (potentially 10- to 100-fold), and using multiple CCD arrays. By the time of the next solar maximum imaging spectroscopy is expected to achieve spatial resolution or approx. arc 1 arc s, temporal resolution or approx. 5 s, and simultaneous critically-sampled spectroscopy of several lines and continua. As a result, continued increase in our understanding of the physical processes and configurations of solar flares in the chromosphere, temperature minimum region, and photosphere can be anticipated. Even greater progress toward a more global understanding of flares will obviously come about when simultaneous optical, X-ray, and gamma-ray imaging spectroscopy are possible.

Canfield, R. C.↗

Active control of structurally-coupled sound fields in elastic cylinders by vibrational force inputs

Active control of structurally-coupled sound fields in elastic cylinders is analytically and experimentally studied. The primary (noise) field in the cylinder model is generated by the coupled dynamic response of the shell under loading by a single exterior acoustic source. Control of the interior sound field is achieved by applying vibrational force inputs directly to the shell wall. Action of the point controllers serve to increase the input impedance of select structural modes of the shell which are well-coupled to the interior acoustic cavity, thus substantially reducing sound transmission into the cavity. Spatially-averaged noise reductions in excess of 30 dB are demonstrated for acoustic resonant conditions within the cavity. Twin controller configurations are presented which demonstrate the ability to independently control orthogonal modes of the interior acoustic space. Benefits and drawbacks of this new methodology for noise control are discussed and clearly demonstrated.

Jones, J. D.↗

Ionospheric response to the sustained high geomagnetic activity during the March '89 great storm

A simulation was conducted to model the high-latitude ionospheric to the sustaied level of high geomagnetic activity for the great magnetic storm period of March 13-14, 1989. The geomagnetic and solar activity indices and the Defense Meterological Satellite Program (DMSP) F8 and F9 satellite data for particle precipitation and high-latitude convection were used as inputs to a time-dependent ionospheric model (TDIM). The results of the TDIM were compared to both DMSP plasma density data and ground-based total electron content (TEC) measurements for the great storm period as well as with earlier storm observations. The comparisons show that the overall structure of the high-latitude ionosphere was dominated by an increased convection speed within the polar cap that led to increased ion temperatures. In turn, this enhanced the NO(+) density, raised the atomic-to-molecular ion transition height to over 300 km, decreased N(sub m)F(sub 2), increased h(sub m)F(sub 2), and in places either increased n(sub e) at 800 km or slightly decreased it. The morphology of the ionosphere under these extreme conditions was considerably different than that modeled for less distributed intervals. These differences included the character of the dayside tongue of ionization that no longer extended deep into the polar cap. Instead, as a result of the ion heating and consequent reduction in N(sub m)F(sub 2), a large polar hole occupied much of the polar region. This polar hole extended beyond the auroral oval and merged with the night sector midatitude trough. The limitaions associated with the applicability of the TDIM to the geomagnetic conditions present on March 13 and 14 are discussed. The primary limitations of the TDIM derive from the limited temporal resolution of the model input parameters and the lack of suitably dynamic thermospheric specification for the great storm conditions. These limitations leads to midlatitude ionospheric storm phases that do no follow those observed.

Sojka, J. J.↗

Neptune Aerocapture Systems Analysis

A Neptune Aerocapture Systems Analysis is completed to determine the feasibility, benefit and risk of an aeroshell aerocapture system for Neptune and to identify technology gaps and technology performance goals. The high fidelity systems analysis is completed by a five center NASA team and includes the following disciplines and analyses: science; mission design; aeroshell configuration screening and definition; interplanetary navigation analyses; atmosphere modeling; computational fluid dynamics for aerodynamic performance and database definition; initial stability analyses; guidance development; atmospheric flight simulation; computational fluid dynamics and radiation analyses for aeroheating environment definition; thermal protection system design, concepts and sizing; mass properties; structures; spacecraft design and packaging; and mass sensitivities. Results show that aerocapture can deliver 1.4 times more mass to Neptune orbit than an all-propulsive system for the same launch vehicle. In addition aerocapture results in a 3-4 year reduction in trip time compared to all-propulsive systems. Aerocapture is feasible and performance is adequate for the Neptune aerocapture mission. Monte Carlo simulation results show 100% successful capture for all cases including conservative assumptions on atmosphere and navigation. Enabling technologies for this mission include TPS manufacturing; and aerothermodynamic methods and validation for determining coupled 3-D convection, radiation and ablation aeroheating rates and loads, and the effects on surface recession.

Lockwood, Mary Kae↗

Radiative Heating of Large Meteoroids During Atmospheric Entry

A high-fidelity approach for simulating the aerothermodynamic environments of meteor entries was developed, which allows the commonly assumed heat transfer coefficient of 0.1 to be assessed. This model uses chemically reacting computational fluid dynamics (CFD), coupled with radiation transport and surface ablation. Coupled radiation accounts for the impact of radiation on the flowfield energy equations, while coupled ablation explicitly models the injection of ablation products within the flowfield and radiation simulations. For a meteoroid with a velocity of 20 km/s, coupled radiation is shown to reduce the stagnation point radiative heating by over 60%. The impact of coupled ablation (with coupled radiation) is shown to provide at least a 70% reduction in the radiative heating relative to cases with only coupled radiation. This large reduction is partially the result of the low ionization energies of meteoric ablation products relative to air species. The low ionization energies of ablation products, such as Mg and Ca, provide strong photoionization and atomic line absorption in regions of the spectrum that air species do not. MgO and CaO are also shown to provide significant absorption. Turbulence is shown to impact the distribution of ablation products through the shock-layer, which results in up to a 100% increase in the radiative heating downstream of the stagnation point. To create a database of heat transfer coefficients, the developed model was applied to a range of cases. This database considered velocities ranging from 14 to 20 km/s, altitudes ranging from 20 to 50 km, and nose radii ranging from 1 to 100 m. The heat transfer coefficients from these simulations are below 0.045 for the range of cases, for both laminar and turbulent, which is significantly lower than the canonical value of 0:1. When the new heat transfer model is applied to a Tunguska-like 15 Mt entry, the effect of the new model is to lower the height of burst by up to 2 km, depending on assumed entry angle. This, in turn, results in a significantly larger ground damage footprint than when the canonical heating assumption is used.

Christopher O Johnston↗

Anisotropic Mach Cone Aligned Mesh Adaptation for Low Boom Simulations

An automated off-body Mach cone aligned structured curvilinear grid generation procedure is presented for near-field computational fluid dynamics simulations. This procedure combines output-based indicators and mesh redistribution to perform anisotropic mesh adaptation while maintaining Mach cone alignment. Automation is achieved through a novel direction-based adaptation indicator formulation. The adaptation procedure is demonstrated on the JAXA Wing Body geometry and X-59 C608 demonstrator model from the Second and Third AIAA Sonic Boom Prediction Workshops, respectively. It is demonstrated that anisotropic mesh adaptation may result in a greater than fifty percent reduction in resource usage required to achieve the same level of accuracy as uniform and user constructed Mach cone aligned grids for near-field pressure signatures, ground-level overpressure signatures, and loudness metrics.

CST↗

Anisotropic Mach Cone Aligned Mesh Adaptation for Low Boom Simulations

An automated off-body Mach cone aligned structured curvilinear grid generation procedure is presented for near-field computational fluid dynamics simulations. This procedure combines output-based indicators and mesh redistribution to perform anisotropic mesh adaptation while maintaining Mach cone alignment. Automation is achieved through a novel direction-based adaptation indicator formulation. The adaptation procedure is demonstrated on the JAXA Wing Body geometry and X-59 C608 demonstrator model from the Second and Third AIAA Sonic Boom Prediction Workshops, respectively. It is demonstrated that anisotropic mesh adaptation may result in a greater than fifty percent reduction in resource usage required to achieve the same level of accuracy as uniform and user constructed Mach cone aligned grids for near-field pressure signatures, ground-level overpressure signatures, and loudness metrics.

CST↗

Assessing the Impact of a Novel TBC Material on Heat Transfer in a Spark Ignition Engine through 3D CFD-FEA Co-Simulation Routine

Thermal barrier coatings (TBCs) have been of interest since the 1970s for application in internal combustion (IC) engines. Thin TBCs exhibit a temperature swing phenomenon wherein wall temperatures dynamically respond to the transient working-gas temperature throughout the engine cycle, thus reducing the temperature difference driving the heat transfer. Determining these varying wall temperatures is necessary to evaluate and study the effect of coatings on wall heat transfer. This study focuses on developing a 3D computational fluid dynamics (CFD)-finite element analysis (FEA) coupled simulation, or co-simulation, routine to determine the wall temperatures of a piston coated with a thin TBC layer subject to spark ignition combustion heat flux. A CONVERGE 3D-CFD model was used to simulate the combustion process in a single-cylinder, light-duty experimental spark ignition (SI) engine. Transient piston heat transfer analysis was conducted using ABAQUS, a FEA package, under the simulated combustion heat flux load. The effect of the temperature swing phenomenon due to this TBC layer was observed in a CFD simulation by implementing the FEA results as the piston thermal boundary conditions. The boundary conditions were passed between the CFD and FEA tools until a quasi-steady state solution was achieved. Furthermore, a reduction in wall heat transfer was observed due to a reduced temperature difference between the wall and the working gas.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

An Empirical Cumulus Parameterization Scheme for a Global Spectral Model

Realistic vertical heating and drying profiles in a cumulus scheme is important for obtaining accurate weather forecasts. A new empirical cumulus parameterization scheme based on a procedure to improve the vertical distribution of heating and moistening over the tropics is developed. The empirical cumulus parameterization scheme (ECPS) utilizes profiles of Tropical Rainfall Measuring Mission (TRMM) based heating and moistening derived from the European Centre for Medium- Range Weather Forecasts (ECMWF) analysis. A dimension reduction technique through rotated principal component analysis (RPCA) is performed on the vertical profiles of heating (Q1) and drying (Q2) over the convective regions of the tropics, to obtain the dominant modes of variability. Analysis suggests that most of the variance associated with the observed profiles can be explained by retaining the first three modes. The ECPS then applies a statistical approach in which Q1 and Q2 are expressed as a linear combination of the first three dominant principal components which distinctly explain variance in the troposphere as a function of the prevalent large-scale dynamics. The principal component (PC) score which quantifies the contribution of each PC to the corresponding loading profile is estimated through a multiple screening regression method which yields the PC score as a function of the large-scale variables. The profiles of Q1 and Q2 thus obtained are found to match well with the observed profiles. The impact of the ECPS is investigated in a series of short range (1-3 day) prediction experiments using the Florida State University global spectral model (FSUGSM, T126L14). Comparisons between short range ECPS forecasts and those with the modified Kuo scheme show a very marked improvement in the skill in ECPS forecasts. This improvement in the forecast skill with ECPS emphasizes the importance of incorporating realistic vertical distributions of heating and drying in the model cumulus scheme. This also suggests that in the absence of explicit models for convection, the proposed statistical scheme improves the modeling of the vertical distribution of heating and moistening in areas of deep convection.

Rajendran, K.↗