Search NASA⌕ Search

SEARCH · Search NASA

Results for “High dimensional 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 379 records · Page 21

Instabilities in the Wake of a Pseudorandom Roughness on a Supersonic Flat Plate

A roughness patch defined by a two-dimensional Fourier series was installed in a highly-polished flat plate model immersed in a supersonic quiet flow. The NASA Langley Research Center Supersonic Low Disturbance Tunnel is capable of providing a Mach 3.5 flow with low freestream acoustic noise and low turbulence levels. Measurements made with a hot-wire probe downstream of the roughness element showed the presence of instability modes that did not cause transition in the measurement region at a freestream unit Reynolds number of Re =12.6×10^6/m. The largest amplitude mode appeared to be similar to a symmetric mode. The same roughness patch immersed in a higher freestream acoustic noise environment showed little to no evidence of the instability modes. Instead, disturbances found in the boundary layer matched closely to the freestream spectra measured in the empty tunnel under the same conditions, indicating some evidence of freestream forcing in the boundary layer.

Roughness induced transition↗

Instabilities in the Wake of a Pseudorandom Roughness on a Supersonic Flat Plate

A roughness patch defined by a two-dimensional Fourier series was installed in a highly-polished flat plate model immersed in a supersonic quiet flow. The NASA Langley Research Center Supersonic Low Disturbance Tunnel is capable of providing a Mach 3.5 flow with low freestream acoustic noise and low turbulence levels. Measurements made with a hot-wire probe downstream of the roughness element showed the presence of instability modes that did not cause transition in the measurement region at a freestream unit Reynolds number of Re =12.6×10^6/m. The largest amplitude mode appeared to be similar to a symmetric mode. The same roughness patch immersed in a higher freestream acoustic noise environment showed little to no evidence of the instability modes. Instead, disturbances found in the boundary layer matched closely to the freestream spectra measured in the empty tunnel under the same conditions, indicating some evidence of freestream forcing in the boundary layer.

Roughness induced transition↗

Convective dynamics - Panel report

Aspects of highly organized forms of deep convection at midlatitudes are reviewed. Past emphasis in field work and cloud modeling has been directed toward severe weather as evidenced by research on tornadoes, hail, and strong surface winds. A number of specific issues concerning future thrusts, tactics, and techniques in convective dynamics are presented. These subjects include; convective modes and parameterization, global structure and scale interaction, convective energetics, transport studies, anvils and scale interaction, and scale selection. Also discussed are analysis workshops, four-dimensional data assimilation, matching models with observations, network Doppler analyses, mesoscale variability, and high-resolution/high-performance Doppler. It is also noted, that, classical surface measurements and soundings, flight-level research aircraft data, passive satellite data, and traditional photogrammetric studies are examples of datasets that require assimilation and integration.

Carbone, Richard↗

Attitude determination of a high altitude balloon system. Part 1: Development of the mathematical model

A mathematical model for predicting the three dimensional motion of the balloon system is developed, which includes the effects of bounce, pendulation and spin of each subsystem. Boundary layer effects are also examined, along with the aerodynamic forces acting on the balloon. Various simplified forms of the system mathematical model were developed, based on an order of magnitude analysis.

Nigro, N. J.↗

Separable physics-informed DeepONet: Breaking the curse of dimensionality in physics-informed machine learning

The deep operator network (DeepONet) has shown remarkable potential in solving partial differential equations (PDEs) by mapping between infinite-dimensional function spaces using labeled datasets. However, in scenarios lacking labeled data, the physics-informed DeepONet (PI-DeepONet) approach, which utilizes the residual loss of the governing PDE to optimize the network parameters, faces significant computational challenges, particularly due to the curse of dimensionality. This limitation has hindered its application to high-dimensional problems, making even standard 3D spatial with 1D temporal problems computationally prohibitive. Additionally, the computational requirement increases exponentially with the discretization density of the domain. Here, to address these challenges and enhance scalability for high-dimensional PDEs, we introduce the Separable physics-informed DeepONet (Sep-PI-DeepONet). This framework employs a factorization technique, utilizing sub-networks for individual one-dimensional coordinates, thereby reducing the number of forward passes and the size of the Jacobian matrix required for gradient computations. By incorporating forward-mode automatic differentiation (AD), we further optimize computational efficiency, achieving linear scaling of computational cost with discretization density and dimensionality, making our approach highly suitable for high-dimensional PDEs. We demonstrate the effectiveness of Sep-PI-DeepONet through three benchmark PDE models: the viscous Burgers’ equation, Biot’s consolidation theory, and a parameterized heat equation. Our framework maintains accuracy comparable to the conventional PI-DeepONet while reducing training time by two orders of magnitude. Notably, for the heat equation solved as a 4D problem, the conventional PI-DeepONet was computationally infeasible (estimated 289.35 h), while the Sep-PI-DeepONet completed training in just 2.5 h. These results underscore the potential of Sep-PI-DeepONet in efficiently solving complex, high-dimensional PDEs, marking a significant advancement in physics-informed machine learning.

Neural operator↗

Flood Hazard Assessment from Storm Tides, Rain and Sea Level Rise for a Tidal River Estuary

Cities and towns along the tidal Hudson River are highly vulnerable to flooding through the combination of storm tides and high streamflows, compounded by sea level rise. Here a three-dimensional hydrodynamic model, validated by comparing peak water levels for 76 historical storms, is applied in a probabilistic flood hazard assessment. In simulations, the model merges streamflows and storm tides from tropical cyclones (TCs), offshore extratropical cyclones (ETCs) and inland "wet extratropical" cyclones (WETCs). The climatology of possible ETC and WETC storm events is represented by historical events (1931-2013), and simulations include gauged streamflows and inferred ungauged streamflows (based on watershed area) for the Hudson River and its tributaries. The TC climatology is created using a stochastic statistical model to represent a wider range of storms than is contained in the historical record. TC streamflow hydrographs are simulated for tributaries spaced along the Hudson, modeled as a function of TC attributes (storm track, sea surface temperature, maximum wind speed) using a statistical Bayesian approach. Results show WETCs are important to flood risk in the upper tidal river (e.g., Albany, New York), ETCs are important in the estuary (e.g., New York City) and lower tidal river, and TCs are important at all locations due to their potential for both high surge and extreme rainfall. The raising of floods by sea level rise is shown to be reduced by approximately 30-60 percent at Albany due to the dominance of streamflow for flood risk. This can be explained with simple channel flow dynamics, in which increased depth throughout the river reduces frictional resistance, thereby reducing the water level slope and the upriver water level.

Tidal river↗

CELSS system control: issues, methods, and directions

In the general control perspective, the CELSS concept implies a very complex system and presents challenges at every level. These challenges are generated by: (1) the prospect that the system will be inherently unstable, (2) the prospective difficulty of establishing an adequate mathematical model of the system for the purpose of control law synthesis (dimensionality is high, and the dynamics and interactive processes of some of the subsystems are not understood well), (3) assuring control law robustness (assuring that the resulting control law(s) will be effective over the domain of the specified uncertainties), (4) hardware realization of the control law, (5) hardware system robustness ("fault tolerance") and (6) achieving the logistics of the automation (or "management") aspects of the problem. A suggested organization of the problem, a sketch of the issues related to perceived difficulties, a commentary/evaluation of the issues, a review of methods available to address the issues, and a suggested strategy to address the broad CELSS systems control problem are presented.

NASA Center ARC↗

Design of a High-Energy, Two-Stage Pulsed Plasma Thruster

Design details of a proposed high-energy (approx. 50 kJ/pulse), two-stage pulsed plasma thruster are presented. The long-term goal of this project is to develop a high-power (approx. 500 kW), high specific impulse (approx. 7500 s), highly efficient (approx. 50%),and mechanically simple thruster for use as primary propulsion in a high-power nuclear electric propulsion system. The proposed thruster (PRC-PPT1) utilizes a valveless, liquid lithium-fed thermal plasma injector (first stage) followed by a high-energy pulsed electromagnetic accelerator (second stage). A numerical circuit model coupled with one-dimensional current sheet dynamics, as well as a numerical MHD simulation, are used to qualitatively predict the thermal plasma injection and current sheet dynamics, as well as to estimate the projected performance of the thruster. A set of further modelling efforts, and the experimental testing of a prototype thruster, is suggested to determine the feasibility of demonstrating a full scale high-power thruster.

Markusic, T. E.↗

Computation of Space Shuttle high-pressure cryogenic turbopump ball bearing two-phase coolant flow

A homogeneous two-phase fluid flow model, implemented in a three-dimensional Navier-Stokes solver using computational fluid dynamics methodology is described. The application of the model to the analysis of the pump-end bearing coolant flow of the high-pressure oxygen turbopump of the Space Shuttle main engine is studied. Results indicate large boiling zones and hot spots near the ball/race contact points. The extent of the phase change of the liquid oxygen coolant flow due to the frictional and viscous heat fluxes near the contact areas has been investigated for the given inlet conditions of the coolant.

Chen, Yen-Sen↗

Multidisciplinary Design Optimization and Analysis of an Open Rotor Stage: Part 1

Successful design of open rotor propulsors requires effective analysis across multiple disciplines, including aerodynamics, acoustics, and structures. A viable design must not only be efficient but must also produce an acceptable level of noise and meet all static and dynamic structural requirements. For design and optimization, this is especially challenging because running high fidelity analyses is resource-intensive, and optimizing a design may require many thousands of cases to be analyzed. For this reason, the NASA team has applied design methodology that utilizes low-cost aerodynamic methods, machine learning models, and high-fidelity analyses when necessary. This includes standard two-dimensional methods such as throughflow analysis and quasi-3D blade-to-blade CFD, as well as some newly developed methods. Optimization using 3D CFD is necessary to maximize performance, and this is considered as well. All optimizations are carried out subject to structural constraints evaluated using finite element analysis. Doing this accurately requires a robust trunnion design, capable of pitching the blade between cruise and takeoff conditions while maintaining acceptable factor of safety. Hot to cold analysis must also be applied in order to correctly determine the as-manufactured shape. For acoustics, the unsteady pressures on the blade surfaces must be predicted, and this can be done either through full-annulus unsteady CFD or through a nonlinear harmonic method (NLH). NLH can run much faster, allowing some acoustic considerations to be made earlier in the design process. The design process is ongoing, and this presentation will review the current status and planned next steps. This part of the talk will focus on aerodynamic performance and be followed by a talk on structures and acoustics.

Design↗

High Reynolds number analysis of flat plate and separated afterbody flow using non-linear turbulence models

The ability of the three-dimensional Navier-Stokes method, PAB3D, to simulate the effect of Reynolds number variation using non-linear explicit algebraic Reynolds stress turbulence modeling was assessed. Subsonic flat plate boundary-layer flow parameters such as normalized velocity distributions, local and average skin friction, and shape factor were compared with DNS calculations and classical theory at various local Reynolds numbers up to 180 million. Additionally, surface pressure coefficient distributions and integrated drag predictions on an axisymmetric nozzle afterbody were compared with experimental data from 10 to 130 million Reynolds number. The high Reynolds data was obtained from the NASA Langley 0.3m Transonic Cryogenic Tunnel. There was generally good agreement of surface static pressure coefficients between the CFD and measurement. The change in pressure coefficient distributions with varying Reynolds number was similar to the experimental data trends, though slightly over-predicting the effect. The computational sensitivity of viscous modeling and turbulence modeling are shown. Integrated afterbody pressure drag was typically slightly lower than the experimental data. The change in afterbody pressure drag with Reynolds number was small both experimentally and computationally, even though the shape of the distribution was somewhat modified with Reynolds number.

Carlson, John R.↗

Towards dynamical low-rank approximation for neutrino kinetic equations. Part I: Analysis of an idealized relaxation model

Dynamical low-rank approximation (DLRA) is an emerging tool for reducing computational costs and provides memory savings when solving high-dimensional problems. Here, in this work, we propose and analyze a semi-implicit dynamical low-rank discontinuous Galerkin (DLR-DG) method for the space homogeneous kinetic equation with a relaxation operator, modeling the emission and absorption of particles by a background medium. Both DLRA and the discontinuous Galerkin (DG) scheme can be formulated as Galerkin equations. To ensure their consistency, a weighted DLRA is introduced so that the resulting DLR-DG solution is a solution to the fully discrete DG scheme in a subspace of the standard DG solution space. Similar to the standard DG method, we show that the proposed DLR-DG method is well-posed. We also identify conditions such that the DLR-DG solution converges to the equilibrium. Numerical results are presented to demonstrate the theoretical findings.

97 MATHEMATICS AND COMPUTING↗

High-frequency ECG

The standard ECG is by convention limited to 0.05-150 Hz, but higher frequencies are also present in the ECG signal. With high-resolution technology, it is possible to record and analyze these higher frequencies. The highest amplitudes of the high-frequency components are found within the QRS complex. In past years, the term "high frequency", "high fidelity", and "wideband electrocardiography" have been used by several investigators to refer to the process of recording ECGs with an extended bandwidth of up to 1000 Hz. Several investigators have tried to analyze HF-QRS with the hope that additional features seen in the QRS complex would provide information enhancing the diagnostic value of the ECG. The development of computerized ECG-recording devices that made it possible to record ECG signals with high resolution in both time and amplitude, as well as better possibilities to store and process the signals digitally, offered new methods for analysis. Different techniques to extract the HF-QRS have been described. Several bandwidths and filter types have been applied for the extraction as well as different signal-averaging techniques for noise reduction. There is no standard method for acquiring and quantifying HF-QRS. The physiological mechanisms underlying HF-QRS are still not fully understood. One theory is that HF-QRS are related to the conduction velocity and the fragmentation of the depolarization wave in the myocardium. In a three-dimensional model of the ventricles with a fractal conduction system it was shown that high numbers of splitting branches are associated with HF-QRS. In this experiment, it was also shown that the changes seen in HF-QRS in patients with myocardial ischemia might be due to the slowing of the conduction velocity in the region of ischemia. This mechanism has been tested by Watanabe et al by infusing sodium channel blockers into the left anterior descending artery in dogs. In their study, 60 unipolar ECGs were recorded from the entire ventricular surface and were signal-averaged and filtered in the 30-250 Hz frequency range. The results showed that the decrease noted in the HF-QRS correlated linearly with the local conduction delay. The results suggest that HF-QRS is a potent indicator of disturbed local conduction. An alternative theory is that HF-QRS reflect the shape of the original electrocardiographic signal. Bennhagen et al showed that root mean square (RMS) voltage values of the depolarization signal correlate poorly with the signal amplitude but highly with the first and second derivatives, i.e. the velocity and the acceleration of the signal. It has also been suggested that the autonomic nervous system affects HF-QRS. For example, sitting up causes significant changes in HF-QRS in some leads compared to the supine position [Douglas et al., 2006]. Unpublished results indicate that familial dysautonomic patients (both vagal and sympathetic degeneration) have very little Reduced Amplitude Zones (RAZ) formation . Athletic individuals, especially elite athletes, who have vagally-mediated changes on the conventional ECG (i.e. early repolarization, bradycardia) have increased RAZ formation. Further electrophysiological studies are needed, however, to better understand the underlying mechanisms of HF-QRS. Several investigators have studied HF-QRS in different cardiac conditions, including acute myocardial ischemia and myocardial infarction (MI). However, in order for clinicians to confidently use HF-QRS as an adjunct to standard ECG, more knowledge about the characteristics of HF-QRS is needed.

Tragardh, Elin↗

Improved Cell Culture Method for Growing Contracting Skeletal Muscle Models

An improved method for culturing immature muscle cells (myoblasts) into a mature skeletal muscle overcomes some of the notable limitations of prior culture methods. The development of the method is a major advance in tissue engineering in that, for the first time, a cell-based model spontaneously fuses and differentiates into masses of highly aligned, contracting myotubes. This method enables (1) the construction of improved two-dimensional (monolayer) skeletal muscle test beds; (2) development of contracting three-dimensional tissue models; and (3) improved transplantable tissues for biomedical and regenerative medicine applications. With adaptation, this method also offers potential application for production of other tissue types (i.e., bone and cardiac) from corresponding precursor cells.

Marquette, Michele L.↗

Refractory-based thermal energy storage for industrial process heat: one-dimensional modeling, control, and optimization

The variable and weather-dependent output of wind and solar power plants present a substantial challenge for planning and operating electricity-systems, particularly in the absence of cost-effective and dispatchable energy storage technologies. This study investigates a high-temperature, electrically heated, refractory-based thermal energy storage (RTES) system that stores electrical energy as sensible heat in dense ceramic bricks over the 950–1800 °C range. The stored heat can be discharged as a controlled hot-gas stream for industrial heating, fuel substitution in high-temperature processes, or electricity generation. The main novelty is a comprehensive modelling, control, mapping, and optimization framework that integrates one-dimensional transient gas–solid heat transfer, fan-assisted discharge, bypass-flow regulation, reheating logic, fan-power evaluation, insulation-loss assessment, and genetic-algorithm-based design optimization. The model uses feedback from outlet temperature and delivered power to regulate discharge, while a two-stage genetic algorithm optimizes brick-channel geometry, gas-flow operation, and multilayer insulation thicknesses. Storage capacities below 50 MWh and discharge powers of 5–30 MW are analyzed to evaluate hold time, thermal delivery, fan-power penalty, heat loss, state-of-charge evolution, and indicative capital cost. Results demonstrate that optimized and well-insulated refractory-based thermal energy storage units can provide stable, efficient, and repeatable heat delivery over multiple discharge cycles. The generated performance and cost maps support modular refractory thermal energy storage as a practical option for large-scale integration of wind and solar generation and for high-temperature industrial process heat.

25 ENERGY STORAGE↗

Aeroelastic Calculations of Quiet High- Speed Fan Performed

An advanced high-speed fan was recently designed under a cooperative effort between the NASA Glenn Research Center and Honeywell Engines & Systems. The principal design goals were to improve performance and to reduce fan noise at takeoff. Scale models of the Quiet High-Speed Fan were tested for operability, performance, and acoustics. During testing, the fan showed significantly improved noise characteristics, but a self-excited aeroelastic vibration known as flutter was encountered in the operating range. Flutter calculations were carried out for the Quiet High-Speed Fan using a three-dimensional, unsteady aerodynamic, Reynolds-averaged Navier-Stokes turbomachinery code named "TURBO." The TURBO code can accurately model the viscous flow effects that can play an important role in various aeroelastic problems such as flutter with flow separation, flutter at high loading conditions near the stall line (stall flutter), and flutter in the presence of shock and boundary-layer interaction. Initially, calculations were performed with no blade vibrations. These calculations were at a constant rotational speed and a varying mass flow rate. The mass flow rate was varied by changing the backpressure at the exit boundary of the computational domain. These initial steady calculations were followed by aeroelastic calculations in which the blades were prescribed to vibrate harmonically in a natural mode, at a natural frequency, and with a fixed interblade phase angle between adjacent blades. The AE-prep preprocessor was used to interpolate the in-vacuum mode shapes from the structural dynamics mesh onto the computational fluid dynamics mesh and to smoothly propagate the grid deformations from the blade surface to the interior points of the grid. The aeroelastic calculations provided the unsteady aerodynamic forces on the blade surface due to blade vibrations. These forces were vector multiplied with the structural dynamic mode shape to calculate the work done on the blade during one vibration period, then this result was converted to an aerodynamic damping. Flutter occurs when the aerodynamic damping becomes negative, if structural damping is ignored. The results of these aeroelastic calculations are summarized in this plot of aerodynamic damping versus mass flow rate at a constant rotational speed. As the backpressure is increased, the mass flow rate through the fan decreases and the fan operating point moves towards the stall line. The aeroelastic calculations showed that the aerodynamic damping decreases as the stall line is approached, as observed during testing. In addition, the aeroelastic calculations with the TURBO code correctly predicted the aeroelastic parameters: the most unstable vibration mode and interblade phase angle, as observed during testing. The Quiet High-Speed Fan demonstrated significant noise reductions during testing, but flutter imposed limits on its operating range. The accurate calculation of the aeroelastic characteristics using the TURBO code is a significant step toward eliminating flutter from the operating range and toward realizing the benefits of reduced fan noise. The aeroelastic calculations described here were performed under a grant by University of Toledo researchers in collaboration with Glenn's researchers.

Bakhle, Milind A.↗

Scale-Similar Models for Large-Eddy Simulations

Scale-similar models employ multiple filtering operations to identify the smallest resolved scales, which have been shown to be the most active in the interaction with the unresolved subgrid scales. They do not assume that the principal axes of the strain-rate tensor are aligned with those of the subgrid-scale stress (SGS) tensor, and allow the explicit calculation of the SGS energy. They can provide backscatter in a numerically stable and physically realistic manner, and predict SGS stresses in regions that are well correlated with the locations where large Reynolds stress occurs. In this paper, eddy viscosity and mixed models, which include an eddy-viscosity part as well as a scale-similar contribution, are applied to the simulation of two flows, a high Reynolds number plane channel flow, and a three-dimensional, nonequilibrium flow. The results show that simulations without models or with the Smagorinsky model are unable to predict nonequilibrium effects. Dynamic models provide an improvement of the results: the adjustment of the coefficient results in more accurate prediction of the perturbation from equilibrium. The Lagrangian-ensemble approach [Meneveau et al., J. Fluid Mech. 319, 353 (1996)] is found to be very beneficial. Models that included a scale-similar term and a dissipative one, as well as the Lagrangian ensemble averaging, gave results in the best agreement with the direct simulation and experimental data.

Sarghini, F.↗

Asymptotic high frequency analysis of the electromagnetic backscattering from an inlet model consisting of piecewise linearly tapered sections

Electromagnetic backscattering from an open ended three dimensional inlet model is analyzed and computed patterns are compared with results of experimental measurements. The model is comprised of two sections. The first section consists of a linearly tapered waveguide with a rectangular opening at one end and the other end is connected to the second section which is a uniform rectangular waveguide with a planar perfectly conducting termination. The model is electrically large so that many propagating modes are excited. The method of analysis contains conventional aperture integration and modal techniques combined with high frequency techniques, which employ concepts such as modal rays, geometrical theory of diffraction and equivalent currents. For the cases considered, it is shown that only a few of the many propagating modes contribute appreciably to the backscattered field. These modes are selected according to their modal ray angle directions.

Altintas, A.↗