Search NASA⌕ Search

SEARCH · Search NASA

Results for “COUPLE MULTIPLE METHOD”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

A fourth order sharp immersed method for the incompressible Navier-Stokes equations with stationary and moving boundaries and interfaces

We propose a fourth order Navier-Stokes solver based on the immersed interface method (IIM), for flow problems with stationary and one-way coupled moving boundaries and interfaces. Our algorithm employs a Runge-Kutta-based projection method that maintains high-order temporal accuracy in both velocity and pressure for steady and unsteady velocity boundary conditions. Fourth order spatial accuracy is achieved through a novel fifth order IIM discretization scheme for the advection term, as well as existing high-order interface-corrected finite difference schemes for the other differential operators. Using a set of manufactured flow problems with stationary and moving boundaries, we demonstrate fourth order convergence of velocity and pressure in the infinity norm, both inside the domain and on the immersed boundaries. The solver’s performance is further validated through a range of practical flow simulations, highlighting its efficiency over a second order scheme. Finally, we showcase the ability of our immersed discretization scheme to handle interface-coupled multiphysics problems by solving a conjugate heat transfer problem with multiple immersed solids. Overall, the proposed approach robustly combines the efficiency of high order discretization schemes with the flexibility of immersed discretizations for flow problems with complex, moving boundaries and interfaces.

42 ENGINEERING↗
NASA NTRS1987

Robust eigenstructure assignment by a projection method: Application using multiple optimization criteria

New ideas which lead to feedback control laws for large flexible structures which are insensitive to model uncertainty are presented. A pole placement method is presented which leads to near-unitary closed loop eigenvectors, and a new method is introduced to design the control while simultaneously considering three competing measures of optimality. Robustness versus Integral algorithms are applicable to at least moderately high-dimensioned systems. In the present discussion, controls for two coupled flexible bodies are considered. A 6x24 gain matrix is designed to control a 12 modes system using 6 actuators. Researchers also developed control laws for the R2P2 simulator at Martin Marietta; in this case 3 actuators are used to control a 12th order system. Simulation studies indicate that researchers indeed achieved robust designs without significant difficulties associated with spillover into the uncontrolled modes. Here, several key ideas and numerical results are given. In the references, details of the formulation, discussions of salient features, and connection to the availiable literature are given.

Junkins, J. L.↗
NASA NTRS1986

The Helmholtz resonance behavior of single and multiple rooms

The results of exploratory measurements of the noise fields inside rooms which are excited to resonance either acoustically or mechanically are presented. The data illustrate the nature and extent of the sound pressure level enhancements in single rooms and also how multiple rooms may resonate by means of either acoustic or mechanical coupling. Sound pressure level enhancements of about 5 dB were measured during resonance of rooms having flexible walls. For such conditions the sound pressure levels in the room were essentially uniform and in phase. Variability of up to 20 dB was measured in a room hallway complex having significant acoustic interactions. Resonant frequency prediction methods which work well at model scale, give only fair results for rooms.

Hubbard, H. H.↗
NASA NTRS1996

Estimation of the longitudinal and lateral-directional aerodynamic parameters from flight data for the NASA F/A-18 HARV

This progress report presents the results of an investigation focused on parameter identification for the NASA F/A-18 HARV. This aircraft was used in the high alpha research program at the NASA Dryden Flight Research Center. In this study the longitudinal and lateral-directional stability derivatives are estimated from flight data using the Maximum Likelihood method coupled with a Newton-Raphson minimization technique. The objective is to estimate an aerodynamic model describing the aircraft dynamics over a range of angle of attack from 5 deg to 60 deg. The mathematical model is built using the traditional static and dynamic derivative buildup. Flight data used in this analysis were from a variety of maneuvers. The longitudinal maneuvers included large amplitude multiple doublets, optimal inputs, frequency sweeps, and pilot pitch stick inputs. The lateral-directional maneuvers consisted of large amplitude multiple doublets, optimal inputs and pilot stick and rudder inputs. The parameter estimation code pEst, developed at NASA Dryden, was used in this investigation. Results of the estimation process from alpha = 5 deg to alpha = 60 deg are presented and discussed.

Napolitano, Marcello R.↗
NASA NTRS2002

System and Method of Locating Lightning Strikes

A system and method of determining locations of lightning strikes has been described. The system includes multiple receivers located around an area of interest, such as a space center or airport. Each receiver monitors both sound and electric fields. The detection of an electric field pulse and a sound wave are used to calculate an area around each receiver in which the lighting is detected. A processor is coupled to the receivers to accurately determine the location of the lighting strike. The processor can manipulate the receiver data to compensate for environmental variables such as wind, temperature, and humidity. Further, each receiver processor can discriminate between distant and local lightning strikes.

Medelius, Pedro J.↗
NASA NTRS2004

Investigation of Navier-Stokes Code Verification and Design Optimization

With rapid progress made in employing computational techniques for various complex Navier-Stokes fluid flow problems, design optimization problems traditionally based on empirical formulations and experiments are now being addressed with the aid of computational fluid dynamics (CFD). To be able to carry out an effective CFD-based optimization study, it is essential that the uncertainty and appropriate confidence limits of the CFD solutions be quantified over the chosen design space. The present dissertation investigates the issues related to code verification, surrogate model-based optimization and sensitivity evaluation. For Navier-Stokes (NS) CFD code verification a least square extrapolation (LSE) method is assessed. This method projects numerically computed NS solutions from multiple, coarser base grids onto a freer grid and improves solution accuracy by minimizing the residual of the discretized NS equations over the projected grid. In this dissertation, the finite volume (FV) formulation is focused on. The interplay between the xi concepts and the outcome of LSE, and the effects of solution gradients and singularities, nonlinear physics, and coupling of flow variables on the effectiveness of LSE are investigated. A CFD-based design optimization of a single element liquid rocket injector is conducted with surrogate models developed using response surface methodology (RSM) based on CFD solutions. The computational model consists of the NS equations, finite rate chemistry, and the k-6 turbulence closure. With the aid of these surrogate models, sensitivity and trade-off analyses are carried out for the injector design whose geometry (hydrogen flow angle, hydrogen and oxygen flow areas and oxygen post tip thickness) is optimized to attain desirable goals in performance (combustion length) and life/survivability (the maximum temperatures on the oxidizer post tip and injector face and a combustion chamber wall temperature). A preliminary multi-objective optimization study is carried out using a geometric mean approach. Following this, sensitivity analyses with the aid of variance-based non-parametric approach and partial correlation coefficients are conducted using data available from surrogate models of the objectives and the multi-objective optima to identify the contribution of the design variables to the objective variability and to analyze the variability of the design variables and the objectives. In summary the present dissertation offers insight into an improved coarse to fine grid extrapolation technique for Navier-Stokes computations and also suggests tools for a designer to conduct design optimization study and related sensitivity analyses for a given design problem.

Vaidyanathan, Rajkumar↗
NASA NTRS2017

Remote Sensing of Multiple Cloud Layer Heights Using Multi-Angular Measurements

Cloud top height (CTH) affects the radiative properties of clouds. Improved CTH observations will allow for improved parameterizations in large-scale models and accurate information on CTH is also important when studying variations in freezing point and cloud microphysics. NASAs airborne Research Scanning Polarimeter (RSP) is able to measure cloud top height using a novel multi-angular contrast approach. For the determination of CTH, a set of consecutive nadir reflectances is selected and the cross-correlations between this set and co-located sets at other viewing angles are calculated for a range of assumed cloud top heights, yielding a correlation profile. Under the assumption that cloud reflectances are isotropic, local peaks in the correlation profile indicate cloud layers. This technique can be applied to every RSP footprint and we demonstrate that detection of multiple peaks in the correlation profile allow retrieval of heights of multiple cloud layers within single RSP footprints. This paper provides an in-depth description of the architecture and performance of the RSPs CTH retrieval technique using data obtained during the Studies of Emissions and Atmospheric Composition, Clouds and Climate Coupling by Regional Surveys (SEAC(exp. 4)RS) campaign. RSP retrieved cloud heights are evaluated using collocated data from the Cloud Physics Lidar (CPL). The method's accuracy associated with the magnitude of correlation, optical thickness, cloud thickness and cloud height are explored. The technique is applied to measurements at a wavelength of 670 nm and 1880 nm and their combination. The 1880-nm band is virtually insensitive to the lower troposphere due to strong water vapor absorption.

Sinclair, Kenneth↗
NASA NTRS2018

Simulation and Analysis of Opportunistic MSPA for Multiple Cubesat Deployments

We describe a software approach for simultaneous demodulation and decoding of multiple frequency- multiplexed spacecraft across multiple ground stations. The approach involves the use of a single ground antenna with wide enough aperture to receive multiple angularly adjacent spacecraft simultaneously. This technique uses a wide-band RF signal digitizer coupled with easily and cheaply duplicated software-receiver modules to independently process the frequency-channelized downlink signals from the spacecraft. This approach relaxes the need for realizable modems at each antenna, for example in the Deep-SpaceNetwork (DSN), and can also function as a delayed data retrieval method for cubesat missions with routine science data return. Thus, we hope that this solution can enable more efficient utilization of DSN assets. We concentrate on a simulation similar to the proposed Exploration Mission 1 (EM-1) mission from a geometric perspective. To that end, the simulated scenario involves ten secondary cubesats deployed and tracked for a span of four days (the modeled motion over the first four days of the mission does not include any Trajectory Correction Maneuvers (TCMs) that may be necessary as the cubesats approach the moon). Transmissions from the cubesats may be received through a combination of DSN ground sites based on visibility. The cubesats are assumed to utilize typical cubesat transmit powers and a waveform similar to that of the Iris radio. One cubesat is chosen as the “target” and is considered tracked by all ground-stations throughout the simulation (i.e., it is consistently at the center of the main lobe of each ground station antenna when the ground station is in view of the cubesat in terms of elevation angle). The simulation effort involves synthesizing a wideband signal that includes the ten cubesats across a large bandwidth due to each cubesat having its own center-frequency in X-band, near 8.4 GHz. Each signal is characterized by its own Doppler-frequency shift and free-space path-loss computed through the underlying geometry of the simulation as well as the cubesat transmit power. Finally, each cubesat signal experiences a unique antenna gain at each ground station due to the underlying antenna pattern and spacecraft-ground station geometry. We establish two interesting findings: First, the link-budget for the EM-1-like scenario is almost completely limited by the angle between the center of the antenna main beam and the cubesat, which means that the signal-to-noise ratio is wholly adequate for demodulation and decoding at the simulated bitrates unless the cubesat exits the main beam. Secondly, and fortunately due to the use of a software radio architecture, we show that it is possible to successfully receive signals from most cubesats for the entire 4-day simulation. Due to the signal-to-noise ratio being sufficient for demodulation even after the side-lobe’s 17dB reduction in SNR, the software radio can still demodulate such cubesats as long as the radio is capable of re-establishing carrier lock as the cubesats leave the main beam and enter the side-lobe. Considering that the target application utilizes offline processing, a loss of lock can be detected and lock can be re-established by iterating over the data. Extensive simulations demonstrate these results.

Abraham, Douglas S.↗
NASA NTRS1993

Time-accurate simulation of a self-excited oscillatory supersonic external flow with a multi-block solution-adaptive mesh algorithm

Results are presented of an investigation of the time-accurate simulation of supersonic unsteady flow oscillations over spike-tipped bodies using the multistage Runge-Kutta scheme coupled with a dynamic solution-adaptive grid algorithm modified for multiblock capabilities. The inviscid fluxes are described by a modified advective upwind split method to obviate the need for artificial dissipation. If a time-varying, solution-adaptive mesh algorithm is incorporated, resolution of the details of the unsteady spike-tipped body flow is improved. The adaptive algorithm is also shown to resolve multiple and diverse features of the flow simultaneously, with the adapted regions in the mesh convecting with these features as they translate.

Ingram, Clint L.↗
NASA NTRS2023

Microstructural Changes and Pb Mobility During the Zircon to Reidite Transformation: Implications for Planetary Impact Chronology

Impact events modify and leave behind a complex history of rock metamorphism on terrestrial planets. Evidence for an impact event may be recorded in physical changes to minerals, such as mineral deformation and formation of high P-T polymorphs, but also in the form of chemical fingerprints, such as enhanced elemental diffusion and isotopic mixing. Here we explore laboratory shock-induced physical and chemical changes to zircon and feldspar, the former of which is of interest because its trace elements abundances and isotope ratios are used extensively in geochemistry and geochronology. To this end, a granular mixture of Bishop Tuff sanidine and Kuehl Lake zircon, both with well characterized Pb isotope compositions, was prepared and then shocked via a flat plate accelerator. The peak pressure of the experiment, as calculated by the impedance matching method, was ~24 GPa although a broader range of P-T conditions is anticipated due to starting sample porosity. Unshocked and shocked materials were characterized via scanning electron microscopy (SEM), electron backscatter diffraction (EBSD), and Raman spectroscopy. These methods show that the starting zircon material had abundant metamict regions, and the conversion of the feldspar to glass in the post-shock material. Analyses of the shocked product also yielded multiple occurrences of the high-pressure ZrSiO 4 polymorph reidite, with some domains up to 300 μm across. The possibility of U-Pb system disturbance was evaluated via laser ablation-inductively coupled plasma-mass spectrometry (LA-ICP-MS) and secondary ion mass spectrometry (SIMS). The isotopic data reveal that disturbance of the U-Pb geochronometer in the reidite was minimal (<2% for the main U-Pb geochronometers). To better constrain the P-T conditions during the shock experiment, we complement impedance matching pressure calculations with iSALE2D impact simulations. The simulated results yield a range of P-T conditions experienced during the experiment and show that much of the sample may have reached >30 GPa, which is consistent with formation of reidite. In the recovered shocked material, we identified lamellae of reidite, some of which interlock with zircon lamellae. Reidite {112} twins were identified, which we interpret to have formed to reduce stress between the crystal structure of the host zircon and reidite. These two findings support the interpretation that shear transformation enabled the transition of zircon to reidite. The size and presence of reidite found here indicate that this phase is probably common in impact-shocked crustal rocks that experienced ~25 to ~35 GPa, especially when the target material has porosity. Additionally, shock loading of the zircon and transformation to reidite at these pressures in porous materials is unlikely to significantly disturb the U-Pb system in zircon and that the reidite inherits the primary U and Pb elemental and isotopic ratios from the zircon.

Bishop Tuff↗
NASA NTRSDate not supplied

On the Implementation and Further Validation of a Time Domain Boundary Element Method Broadband Impedance Boundary Condition

A time domain boundary integral equation with Burton-Miller reformulation is presented for acoustic scattering by surfaces with liners in a uniform mean flow. The Ingard-Myers impedance boundary condition is implemented using a broadband multipole impedance model and converted into time domain differential equations to augment the boundary integral equation. The coupled integral-differential equations are solved numerically by a March-On-in-Time (MOT) scheme. While the Ingard-Myers condition is known to support Kelvin-Helmholtz instability due to its use of a vortex sheet interface between the flow and the liner surface, it is found that by neglecting a second derivative term in the current time domain impedance boundary condition formulation, the instability an be effectively suppressed in computation. The proposed formulation and implementation are validated using a NASA Langley Research Center Grazing Flow Impedance Tube (GFIT) experimental dataset with satisfactory results. Moreover, a minimization procedure for finding the poles and coefficients of the broadband multiple impedance model is formulated in this paper by which, unlike the commonly used vector-fitting method, passivity of the model is ensured. Numerical tests show the proposed minimization approach is effective for modeling liners that are commonly used in aeroacoustic applications.

Time Domain Impedance Boundary Condition↗
DOE OSTI2025

Collective Modes in Multilayer Graphene/𝛼−RuCl 3 Heterostructures

Collective modes in multilayer graphene, such as plasmons and phonons, exhibit sensitivity to displacement fields and interlayer coupling, distinguishing them from their counterparts in single-layer graphene. Here, we engineer collective modes in charge-transfer heterostructures composed of multilayer graphene and 𝛼−RuCl 3 . In heterostructures with a single 𝛼−RuCl 3 interface, the charge transfer generates displacement fields up to 7 V/nm at the interface between 𝛼−RuCl 3 and the adjacent graphene layer—the highest value achieved through charge-transfer methods. As a result of the broken inversion symmetry, we discover enhanced nonlinear optical response and modified phonon selection rules. Conversely, we find that multilayer graphene sandwiched between two 𝛼−RuCl 3 flakes causes displacement fields to cancel. There, we achieve carrier densities as high as 8 ×10 13 cm −2 in multilayer graphene and restore the phonon selection rules to their unperturbed state. Meanwhile, we demonstrate that plasmonic properties derive from the depletion of multiple valence bands. As a result of the quasilinear band dispersion, these “Dirac multiband plasmons” are relatively unaffected by displacement fields. On the other hand, the inverted heterostructure sequence—two multilayer graphene sheets encapsulating 𝛼−RuCl 3 —activates significant alteration of the plasmons via interlayer plasmon-plasmon coupling. Hence, multilayer graphene and 𝛼−RuCl 3 heterostructures offer a gate-free platform for engineering collective modes derived from inversion symmetry and interlayer coupling.

Moore, Samuel L. [Columbia University, New York, N↗
NASA NTRS2025

PSE-Based Aerodynamic Flow Transition Prediction Using Automated Unstructured CFD Integration

Accurate, robust, and efficient prediction of transition in viscous flows is a significant challenge in computational fluid dynamics. We present a coupled, high-fidelity, iterative approach that leverages the FUN3D flow solver and the LASTRAC stability code to predict transition in low-disturbance environments, initiated by the linear growth of boundary-layer instability modes. Our method integrates the ability of FUN3D to compute mixed laminar-transitional-turbulent mean flows via transition-sensitized Reynolds-Averaged Navier-Stokes equations with the ability of LASTRAC to perform linear stability analysis, all within an automated framework that requires no intermediate user involvement. Unlike conventional frameworks that rely on classical stability theory or reduced-order metamodels, our approach employs the parabolized stability equations to provide more accurate and reliable estimates of disturbance growth for multiple instability mechanisms, including Tollmien-Schlichting, Kelvin-Helmholtz, and crossflow modes. By accounting for the effects of mean-flow nonparallelism as well as the surface curvature, this approach lays the foundation for improved N-factor correlations for transition onset prediction in a broad class of flows. We apply this method to three distinct flow configurations: 1) flow over a zero-pressure-gradient flat plate, 2) the NLF-0416 airfoil with both natural and separation-induced transitions, and 3) a 6:1 prolate spheroid, where transition is primarily driven by crossflow instability. For the two-dimensional cases, a formulated intermittency distribution is used to model the transition zone in between the laminar and fully turbulent flows. The results include comparisons with experimental measurements, similar numerical approaches, and transport-equations-based models, demonstrating good agreement in surface pressure coefficients, transition onset locations, and skin-friction coefficients for all three configurations. Besides contributing a couple of new insights into boundary-layer transition in these canonical cases, this study provides a powerful tool for transition modeling in both research and design applications in aerodynamics.

Transition Prediction↗
NASA NTRS1972

Liquid-propellant dynamics and suppression

Technology related to liquid-propellant interactions with space shuttle vehicles is reviewed. Potential problems unique to the shuttle include liquid-structure interactions resulting from coupled lateral and longitudinal deformations, traveling wave phenomena at shallow propellant levels, and liquid impact during abort, staging, or docking. Technology efforts to define the slosh dynamics under shuttle operating conditions are described with emphasis on analytical representations of the liquid by finite-element and marker-and-cell methods. In addition, slosh suppression is discussed and includes baffle damping and pressure loads for tanks fitted with multiple baffles both above and below the undisturbed liquid surface.

Stephens, D. G.↗
NASA NTRS1975

Calibration of large area scintillation counters on cosmic rays

Construction and calibration are reported for a 6-sq m scintillation counter consisting of two pieces of clear Lucite optically joined along a common edge, sixteen squares of plastic scintillator optically coupled on top of the Lucite, and eight photomultiplier tubes for light collection. Calibration of this counter using an anticoincidence counter to exclude multiple particle events is described. Theoretical calculations for purely electromagnetic extensive air showers indicate that this counter should detect an average of 2.5 particles per event, while the calibration data yield about 1.8 particles per event. It is shown that the anticoincidence-counter calibration method must be used for counters with areas of more than a few square meters if calibration on the response from a single particle of an extensive air shower is desired.

Saterlie, S. F.↗
NASA NTRS1984

Two-dimensional computer simulation of EMVJ and grating solar cells under AMO illumination

A computer program, SCAP2D (Solar Cell Analysis Program in 2-Dimensions), is used to evaluate the Etched Multiple Vertical Junction (EMVJ) and grating solar cells. The aim is to demonstrate how SCAP2D can be used to evaluate cell designs. The cell designs studied are by no means optimal designs. The SCAP2D program solves the three coupled, nonlinear partial differential equations, Poisson's Equation and the hole and electron continuity equations, simultaneously in two-dimensions using finite differences to discretize the equations and Newton's Method to linearize them. The variables solved for are the electrostatic potential and the hole and electron concentrations. Each linear system of equations is solved directly by Gaussian Elimination. Convergence of the Newton Iteration is assumed when the largest correction to the electrostatic potential or hole or electron quasi-potential is less than some predetermined error. A typical problem involves 2000 nodes with a Jacobi matrix of order 6000 and a bandwidth of 243.

Gray, J. L.↗
NASA NTRS2006

Electromagnetic Considerations for Planar Bolometer Arrays in the Single Mode Limit

Filled arrays of planar bolometers are finding astronomical applications at wavelengths as long as several millimeters. In an effort to keep focal planes to a reasonable size while maintaining large numbers of detectors, a common strategy is to push these arrays to operate close to or at the single mode limit. Doing so introduces several new challenges that are not experienced in the multi-mode case of far-infrared detectors having similar pixel sizes. First, diffractive effects of the pixels themselves are no longer insignificant and will ultimately contribute to the resolution limit of the optical system in which they reside. We use the method of Withlngton et al. (2003) to model the polarized diffraction in this limit. Second, it is necessary to re-examine the coupling between the radiation and the absorbing element that is thermally connected to the bolometers. The small f-numbers that are often employed to make use of large focal planes makes backshort construction problematic. We introduce a new strategy to increase detector efficiency that uses an antireflective layer on the front side of the detector array. In addition, typical methods for stray light control that rely on multiple reflections in a lossy medium fail due to physical size constraints. For this application, we find that resonant absorbers are a more effective strategy that can be implemented in the space available.

Wollack, Edward J.↗
NASA NTRSDate not supplied

Advanced Modeling of Control-Structure Interaction in Thrust Vector Control Systems

The Space Launch System (SLS) Core Stage (CS) Thrust Vector Control (TVC) system is comprised of 8 mechanical feedback Shuttle heritage Type III TVC actuators and four RS-25 engines, each attached to a Shuttle heritage gimbal block/bearing. Two actuators are used to move each engine in two planes perpendicular to one another (i.e., pitch and yaw). The TVC system design leverages hardware from the Space Shuttle program as well as new hardware designed specifically for the Core Stage. During the development of the SLS TVC system, a family of advanced dynamics models were developed to extend and compliment the simplified quasi-linear “simplex” model historically used for flight control design and stability analysis. The importance of these advanced models became increasingly evident after ambient and hot fire testing of the Core Stage, which revealed a number of findings associated with the dynamic response of the TVC integrated system. Test responses suggested that the TVC did not meet its performance specifications and its step and frequency responses exhibited unexpected departures from prior lab tests and modeled behavior. One driving factor for these results was a higher-than-expected degree of coupling between the TVC system, the engine dynamics, and the Core Stage structure. This paper is the third installment in a seven-paper series surveying the design, engineering, test validation, and flight performance of the Core Stage Thrust Vector Control system. In this paper, a new method of modeling rocket vehicle thrust vectoring servoelastic dynamics is presented. In this approach, the load dynamics are replaced by a detailed finite element model containing both the rigid body and elastic modes. A partitioning technique is used to compute the effective compliance from the modal data and obtain accurate simulation results using a reduced number of generalized coordinates. Coupled backup structure and nozzle attach compliance effects on multiple engines are captured in higher fidelity than with a spring approximation, eliciting novel effects due to the complex load paths involved in the Core Stage structure. Validation of the model is demonstrated using a variety of structural/modal, laboratory, and full-scale hot fire test data.

Launch Vehicles↗