Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Coupled Simulation”

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 271 records · Page 15

Vibrations of three-dimensional pipe systems with acoustic coupling

A general algorithm is developed to calculate the beam-type dynamic response of three dimensional multiplane finite length pipe systems, consisting of elbow and straight ducts with continuous interfaces. Emphasis is on secondary acoustic wave effects giving rise to coupling mechanisms; and the simulation accounts for one-dimensional elastoacoustic coupling from a plane acoustic wave and secondary loads resulting from wave asymmetries. The transfer matrix approach is adopted in modeling the elastodynamics of each duct, with allowance for distribution loads. Secondary loads from plane wave distortion are considered with a solution of the Helmholtz equation in an equivalent rigid waveguide, and effects of path imperfection are introduced as a perturbation from the hypothetical perfectly straight pipe. Computations indicate that the one-dimensional acoustic assumption is valid for frequencies below one-half the first cut-off frequency, and the three-dimensional acoustic effects produce an increase in response levels near and above cut-off.

El-Raheb, M.↗

Numerical Study of Combustor-Turbine Interactions Considering a Two-Stage High-Pressure Turbine

In this study, preliminary results are presented from the Open National Combustion Code (OpenNCC) applied to the two-stage high-pressure turbine (HTP) from the energy efficient engine (E3) program. Although the combustor is designed for efficient fuel–air mixing and a relatively uniform exit temperature profile, temporal and spatial variation in combustor exit conditions inevitably lead to hot-streaks in the HPT. These hot streaks can cause local hot spots on the HPT blade surfaces that reduce blade life. To improve understanding of these effects and the overall design of combustor and HPT, we are developing a fully-coupled combustor-turbine computational fluid dynamics (CFD) capability. Previous work has simulated the E^3 combustor at the simulated Sea-Level Take-off (SLTO) condition (P(41) = 2.4 [atm]) with and without the first-stage HPT stator included. This work addresses HPT simulations using inflow conditions based on previous simulation results of the E^3 combustor, as well as investigating the effects of spatial and temporal temperature non-uniformity of the combustor-exit flow on the HPT performance for a more realistic SLTO (P(41) = 27.4 [atm])). First, we analyze existing results from the combustor and the HPT 1st-stage stator simulation, where the unsteady hot gas coming out of the combustor exit travels between the stationary stators and then directly exits the computational domain (i.e., no effect from the rotors). Second, we simulate the HPT including the 1st and 2nd stage stators and rotors with the time-averaged, spatially-nonuniform inflow to the HPT taken from the combustor (with HPT 1st-stage stators) solution. Differences in simulations results for the flow field in the 1st-stage stator region are presented. The results of these two simulations (combustor and HPT simulated in a sequential manner) will be compared to a fully-coupled EEE combustor/HPT simulation to understand the impacts of fully-coupled simulations on combustor and HPT performance for two different operating conditions.

LES↗

Numerical Study of Combustor-Turbine Interactions Considering a Two-Stage High-Pressure Turbine

In this study, preliminary results are presented from the Open National Combustion Code (OpenNCC) applied to the two-stage high-pressure turbine (HTP) from the energy efficient engine (E3) program. Although the combustor is designed for efficient fuel–air mixing and a relatively uniform exit temperature profile, temporal and spatial variation in combustor exit conditions inevitably lead to hot-streaks in the HPT. These hot streaks can cause local hot spots on the HPT blade surfaces that reduce blade life. To improve understanding of these effects and the overall design of combustor and HPT, we are developing a fully-coupled combustor-turbine computational fluid dynamics (CFD) capability. Previous work has simulated the E3 combustor at the simulated Sea-Level Take-off (SLTO) condition (P41 = 2.4 [atm]) with and without the first-stage HPT stator included. This work addresses HPT simulations using inflow conditions based on previous simulation results of the E3 combustor, as well as investigating the effects of spatial and temporal temperature non-uniformity of the combustor-exit flow on the HPT performance for a more realistic SLTO (P41 = 27.4 [atm])). First, we analyze existing results from the combustor and the HPT 1st-stage stator simulation, where the unsteady hot gas coming out of the combustor exit travels between the stationary stators and then directly exits the computational domain (i.e., no effect from the rotors). Second, we simulate the HPT including the 1st and 2nd stage stators and rotors with the time-averaged, spatially-nonuniform inflow to the HPT taken from the combustor (with HPT 1st-stage stators) solution. Differences in simulations results for the flow field in the 1st-stage stator region are presented. The results of these two simulations (combustor and HPT simulated in a sequential manner) will be compared to a fully-coupled EEE combustor/HPT simulation to understand the impacts of fully-coupled simulations on combustor and HPT performance for two different operating conditions.

LES↗

Coupling of Charged Particles Via Coulombic Interactions: Numerical Simulations and Resultant Kappa-Like Velocity Space Distribution Functions

A parametric study is performed using the electrostatic simulations of Randol and Christian (2014) in which the number density, n, and initial thermal speed, theta, are varied. The range of parameters covers an extremely broad plasma regime, all the way from the very weak coupling of space plasmas to the very strong coupling of solid plasmas. The first result is that simulations at the same Lambda(sub D), where Lambda(sub D) is the plasma coupling parameter, but at different combinations of n and theta, behave exactly the same. As a function of Lambda(sub D), the form of p(v), the velocity distribution function of v, the magnitude of v, the velocity vector, is studied. For intermediate to high D, heating is observed in p(v) that obeys conservation of energy, and a suprathermal tail is formed, with a spectral index that depends on Lambda(sub D). For strong coupling (Lambda(sub D) much > 1), the form of the tail is v5, consistent with the findings of Randol and Christian (2014). For weak coupling (Lambda(sub D much <1), no acceleration or heating occurs, as there is no free energy. The dependence on N, the number of particles in the simulation, is also explored. There is a subtle dependence in the index of the tail, such that v5 appears to be the N approaches infinity limit.

Randol, Brent M.↗

Precision Landing Performance of a Human-Scale Lunar Lander Using a Generalized Simulation Framework

NASA has established goals of returning humans to Moon with an initial landing by 2024 and a subsequent sustained presence by 2028, which will require technological advances in spacecraft navigation to enable precision landing. The ability to assess the navigation performance of these new and existing technologies is critical to identifying areas of risk reduction and investment. To that end, the Safe and Precise Landing Integrated Capabilities Evolution (SPLICE) project has demonstrated that a detailed six degree-of-freedom integrated performance simulation framework can provide information on and assessment of expected navigation performance. This framework incorporates engineering models of the on-board spacecraft guidance, navigation, and control systems at varying levels of fidelity. Recent advances in the development of this integrated performance simulation permit running these systems “in-the-loop,” rather than assuming perfect knowledge of the spacecraft states. This development, coupled with fast simulation time and modularization of the various system models, enables a wide variety of system trades to be assessed at once. This paper presents a summary of the advances in the SPLICE simulation framework, updates to the spacecraft navigation models, and an application of the framework to characterize the precision landing performance of a human-scale lunar lander. A series of trade studies examining effects of ground state update qualities shows that given all other assumptions, sufficiently accurate Deep Space Network (DSN) measurements can enable safe and precise human-scale Lunar landings.

Spacecraft navigation↗

Precision Landing Performance and Technology Assessments of a Human-Scale Lunar Lander Using a Generalized Simulation Framework

NASA has established goals of returning humans to Moon with an initial landing by 2024 and a subsequent sustained presence by 2028, which will require technological advances in spacecraft navigation to enable precision landing. The ability to assess the navigation performance of these new and existing technologies is critical to identifying areas of risk reduction and investment. To that end, the Safe and Precise Landing Integrated Capabilities Evolution (SPLICE) project has demonstrated that a detailed six degree-of-freedom integrated performance simulation framework can provide information on and assessment of expected navigation performance. This framework incorporates engineering models of the on-board spacecraft guidance, navigation, and control systems at varying levels of fidelity. Recent advances in the development of this integrated performance simulation permit running these systems “in-the-loop,” rather than assuming perfect knowledge of the spacecraft states. This development, coupled with fast simulation time and modularization of the various system models, enables a wide variety of system trades to be assessed at once. This paper presents a summary of the advances in the SPLICE simulation framework, updates to the spacecraft navigation models, and an application of the framework to characterize the precision landing performance of a human-scale lunar lander. A series of trade studies examining effects of ground state update qualities shows that given all other assumptions, sufficiently accurate Deep Space Network (DSN) measurements can enable safe and precise human-scale Lunar landings.

"D'Souza, Sarah", 'Pensado, Alegandro R.↗

Coupled Fluid-Structure Interaction Analysis of Solid Rocket Motor with Flexible Inhibitors

Flexible inhibitors are generally used in solid rocket motors (SRMs) as a means to control the burning of propellant. Vortices generated by the flow of propellant around the flexible inhibitors have been identified as a driving source of instabilities that can lead to thrust oscillations in launch vehicles. Potential coupling between the SRM thrust oscillations and structural vibration modes is an important risk factor in launch vehicle design. As a means to predict and better understand these phenomena, a multidisciplinary simulation capability that couples the NASA production CFD code, Loci/CHEM, with CFDRC's structural finite element code, CoBi, has been developed. This capability is crucial to the development of NASA's new space launch system (SLS). This paper summarizes the efforts in applying the coupled software to demonstrate and investigate fluid-structure interaction (FSI) phenomena between pressure waves and flexible inhibitors inside reusable solid rocket motors (RSRMs). The features of the fluid and structural solvers are described in detail, and the coupling methodology and interfacial continuity requirements are then presented in a general Eulerian-Lagrangian framework. The simulations presented herein utilize production level CFD with hybrid RANS/LES turbulence modeling and grid resolution in excess of 80 million cells. The fluid domain in the SRM is discretized using a general mixed polyhedral unstructured mesh, while full 3D shell elements are utilized in the structural domain for the flexible inhibitors. Verifications against analytical solutions for a structural model under a steady uniform pressure condition and under dynamic modal analysis show excellent agreement in terms of displacement distribution and eigenmode frequencies. The preliminary coupled results indicate that due to acoustic coupling, the dynamics of one of the more flexible inhibitors shift from its first modal frequency to the first acoustic frequency of the solid rocket motor. This insight could have profound implications for SRM and flexible inhibitor designs for current and future launch vehicles including SLS.

Yang, H. Q.↗

Radiation-Spray Coupling for Realistic Flow Configurations

Three Large Eddy Simulations (LES) for a lean-direct injection (LDI) combustor are performed and compared. In addition to the cold flow simulation, the effect of radiation coupling with the multi-physics reactive flow is analyzed. The flame let progress variable approach is used as a subgrid combustion model combined with a stochastic subgrid model for spray atomization and an optically thin radiation model. For accurate chemistry modeling, a detailed Jet-A surrogate mechanism is utilized. To achieve realistic inflow, a simple recycling technique is performed at the inflow section upstream of the swirler. Good comparison is shown with the experimental data mean and root mean square profiles. The effect of combustion is found to change the shape and size of the central recirculation zone. Radiation is found to change the spray dynamics and atomization by changing the heat release distribution and the local temperature values impacting the evaporation process. The simulation with radiation modeling shows wider range of droplet size distribution by altering the evaporation rate. The current study proves the importance of radiation modeling for accurate prediction in realistic spray combustion configurations, even for low pressure systems.

El-Asrag, Hossam↗

Use of Data to Improve El Nino Southern Oscillation Simulations and Forecasts

Various data over the period 1980-1998 have been used to improve El Nino Southern Oscillation (ENSO) simulations and forecasts performed by coupled ocean/atmosphere models. For such a study, oceanic and atmospheric satellite data are necessary. Models assume a baroclinic ocean with a mixed layer coupled to a baroclinic atmosphere to simulate the anomalous fields. The various data have been used to revisit the parameterizations of the subsurface temperature in the mixed layer, atmospheric convection, and friction in the baroclinic ocean. This new model, named Tsub.Conv, simulates 4-year oscillations like the model with the original parameterization, but with anomalous thermocline and wind anomalies that agree a lot better with reality. The atmospheric model still presents some deficiencies in the eastern Pacific and when it is replaced by a statistical atmosphere, the coupled model (Tsub.Astat) performs even better. It is found that the off-equatorial ocean and wind curl anomalies play a crucial role for all these models and that Tsub.Astat simulates an ocean heat content recharge coming from the north prior to a warm event like the ones observed during the 1985-1992 period. The ENSO forecasts over 1980-1998 are considerably improved with Tsub.Astat. One needs to add the simulation of westerly wind bursts to be able to forecast the 1997 event one year in advance. Additional information is contained in the original.

Perigaud, Claire↗

One-Dimensional Coupled Ecosystem-Carbon Flux Model for the Simulation of Biogeochemical Parameters at Ocean Weather Station P

In this Technical Publication, we describe the model functionality and analyze its application to the seasonal and interannual variations of phytoplankton, nutrients, pCO2 and CO2 concentrations in the eastern subarctic Pacific at Ocean Weather Station P (OWSP, 50 deg. N 145 deg. W). We use a verified one-dimensional ecosystem model, coupled with newly incorporated carbon flux and carbon chemistry components, to simulate 22 years (1958-1980) of pCO2 and CO2 variability at Ocean Weather Station P (OWS P). This relatively long period of simulation verifies and extends the findings of previous studies using an explicit approach for the biological component and realistic coupling with the carbon flux dynamics. The slow currents and the horizontally homogeneous ocean in the subarctic Pacific make OWS P one of the best available candidates for modeling the chemistry of the upper ocean in one dimension. The chlorophyll and ocean currents composite for 1998 illustrates this premise. The chlorophyll concentration map was derived from SeaWiFS data and the currents are from an OGCM simulation (from R. Murtugudde).

Signorini, S.↗

Evaluation of Cloud Microphysics Simulated using a Meso-Scale Model Coupled with a Spectral Bin Microphysical Scheme through Comparison with Observation Data by Ship-Borne Doppler and Space-Borne W-Band Radars

Equivalent radar reflectivity factors (Ze) measured by W-band radars are directly compared with the corresponding values calculated from a three-dimensional non-hydrostatic meso-scale model coupled with a spectral-bin-microphysical (SBM) scheme for cloud. Three case studies are the objects of this research: one targets a part of ship-borne observation using 95 GHz Doppler radar over the Pacific Ocean near Japan in May 2001; other two are aimed at two short segments of space-borne observation by the cloud profiling radar on CloudSat in November 2006. The numerical weather prediction (NWP) simulations reproduce general features of vertical structures of Ze and Doppler velocity. A main problem in the reproducibility is an overestimation of Ze in ice cloud layers. A frequency analysis shows a strong correlation between ice water contents (IWC) and Ze in the simulation; this characteristic is similar to those shown in prior on-site studies. From comparing with the empirical correlations by the prior studies, the simulated Ze is overestimated than the corresponding values in the studies at the same IWC. Whereas the comparison of Doppler velocities suggests that large-size snowflakes are necessary for producing large velocities under the freezing level and hence rules out the possibility that an overestimation of snow size causes the overestimation of Ze. Based on the results of several sensitivity tests, we conclude that the source of the overestimation is a bias in the microphysical calculation of Ze or an overestimation of IWC. To identify the source of the problems needs further validation research with other follow-up observations.

Iguchi, T.↗

Goal-Oriented Discretization Error Control in Coupled Nearfield-Farfield Low-Boom Simulations

The method of adjoint weighted residuals is used to determine the level of discretization error in loudness predictions of sonic booms on the ground. We analyze the standard nearfied-farfield domain decomposition approach. In the nearfield domain, the three-dimensional Euler equations are solved to obtain a pressure signature generated by the aircraft. In the farfield, this waveform is propagated through the atmosphere to the ground by solving the augmented Burgers’ equation. Loudness is characterized using weighted sound-exposure-level metrics. We formulate discretization error estimates for the ground signature and loudness metrics for this one-way coupled system. Although the nearfield solution is independent of the farfield, the adjoint formulation for the coupled system provides feedback from the farfield to identify high-error regions in the nearfield. The results demonstrate that the discrete adjoint implementation is asymptotically consistent and provides reliable error estimates. Furthermore, we show how the error can be controlled through adaptive refinement of the nearfield mesh. The approach is evaluated on two- and three-dimensional problems, including the X-59 flight demonstration aircraft.

CST↗

Rheoceptive mediators of graviperception in a water flea: Morphological implications of antennal-socket setae in daphnia magna

Aquatic microcrustaceans of the genus Daphnia are known to orient to light during the day. At night, in the absence of visual cues, daphnids were suspected of maintaining equilibrium by monitoring the direction of gravity through their swimming antennae. Recent investigations using simulated, weightlessness conditions coupled with absence of illumination revealed hair like structures or setae on the basal, articulating socket of the antennae that, when surgically removed, resulted in disorientation. Given the simulated weightlessness or neutrally buoyant condition that eliminated sinking of the normally negatively buoyant Daphnia, it was proposed that the antennal socket setae function as rheoceptors stimulated by the upward rush of water currents during gravity induced, sinking phase of daphnid swimming movements. This rheoceptively mediated, gravity perception hypothesis is further supported by morphological investigations. Scanning electron micrographs indicate that antennal socket setae are anatomically similar to proprioceptors used by higher crustaceans to monitor gravitational direction.

Meyers, D. G.↗

Cable coupling lightning transient qualification

Simulated lightning strike testing of instrumentation cabling on the redesigned solid rocket motor was performed. Testing consisted of subjecting the lightning evaluation test article to simulated lightning strikes and evaluating the effects of instrumentation cable transients on cables within the system tunnel. The maximum short-circuit current induced onto a United Space Boosters, Inc., operational flight cable within the systems tunnel was 92 A, and the maximum induced open-circuit voltage was 316 V. These levels were extrapolated to the worst-case (200 kA) condition of NASA specification NSTS 07636 and were also scaled to full-scale redesigned solid rocket motor dimensions. Testing showed that voltage coupling to cables within the systems tunnel can be reduced 40 to 90 dB and that current coupling to cables within the systems tunnel can be reduced 30 to 70 dB with the use of braided metallic sock shields around cables that are external to the systems tunnel. Testing also showed that current and voltage levels induced onto cables within the systems tunnel are partially dependant on the cables' relative locations within the systems tunnel. Results of current injections to the systems tunnel indicate that the dominant coupling mode on cables within the systems tunnel is not from instrumentation cables but from coupling through the systems tunnel cover seam apertures. It is recommended that methods of improving the electrical bonding between individual sections of the systems tunnel covers be evaluated. Further testing to better characterize redesigned solid rocket motor cable coupling effects as an aid in developing methods to reduce coupling levels, particularly with respect to cable placement within the systems tunnel, is also recommended.

Cook, M.↗

LDEXPT, an intelligent database system for the Composite Load Spectra project

The Composite Load Spectra project develops probabilistic models to simulate the probabilistic loads for selected components of a generic space propulsion system. Tremendous information such as engine load variables and their distributions is needed by the simulation program. An intelligent data base system was constructed and integrated with the probabilistic load simulation program to manage and maintain the knowledge base of the Composite Load Spectra project. The intelligent data base system takes care of the data retrieval and storage functions and has expert knowledge on engine load models and associated engine variables. The integration of the intelligent data base into the load simulation program achieves a smooth coupling between the numeric processing (load simulation calculation) and the symbolic processing (intelligent load information management).

Ho, H.↗

Internal Variability and Disequilibrium Confound Estimates of Climate Sensitivity from Observations

An emerging literature suggests that estimates of equilibrium climate sensitivity (ECS) derived from recent observations and energy balance models are biased low because models project more positive climate feedback in the far future. Here we use simulations from the Coupled Model Intercomparison Project Phase 5 (CMIP5) to show that across models, ECS inferred from the recent historical period (1979-2005) is indeed almost uniformly lower than that inferred from simulations subject to abrupt increases in CO2-radiative forcing. However, ECS inferred from simulations in which sea surface temperatures are prescribed according to observations is lower still. ECS inferred from simulations with prescribed sea surface temperatures is strongly linked to changes to tropical marine low clouds. However, feedbacks from these clouds are a weak constraint on long-term model ECS. One interpretation is that observations of recent climate changes constitute a poor direct proxy for long-term sensitivity.

Marvel, Kate↗

A Flexible Multibody Framework for Space Launch System Liftoff with Umbilical Separations

ICPSU clearances to the lifting SLS vehicle were computed utilizing a coupled system nonlinear simulation approach in order to accurately capture the impact of system flexibilities, SLS forcing functions, pad separation (twang), ML dynamics and timings • i.e., integrate ICPSU model into a SLS nonlinear Coupled Loads Analysis (CLA) This required utilizing a flexible multibody framework for the CLA and the integration of large displacements/rotations capabilities into that CLA • Gantry rotation • Disconnect and retraction of flexible hoses

Arya Majed↗

SST Variation Due to Interactive Convective-Radiative Processes

The recent linking of Cloud-Resolving Models (CRMs) to Ocean-Mixed Layer (OML) models has provided a powerful new means of quantifying the role of cloud systems in ocean-atmosphere coupling. This is due to the fact that the CRM can better resolve clouds and cloud systems and allow for explicit cloud-radiation interaction. For example, Anderson (1997) applied an atmospheric forcing associated with a CRM simulated squall line to a 3-D OML model (one way or passive interaction). His results suggested that the spatial variability resulting from the squall forcing can last at least 24 hours when forced with otherwise spatially uniform fluxes. In addition, the sea surface salinity (SSS) variability continuously decreased following the forcing, while some of the SST variability remained when a diurnal mixed layer capped off the surface structure. The forcing used in the OML model, however, focused on shorter time (8 h) and smaller spatial scales (100-120 km). In this study, the 3-D Goddard Cumulus Ensemble Model (GCE; 512 x 512 x 23 cu km, 2-km horizontal resolution) is used to simulate convective active episodes occurring in the Western Pacific warm pool and Eastern Atlantic regions. The model is integrated for seven days, and the simulated results are coupled to an OML model to better understand the impact of precipitation and changes in the planetary boundary layer upon SST variation. We will specifically examine and compare the results of linking the OML model with various spatially-averaged outputs from GCE simulations (i.e., 2 km vs. 10-50 km horizontal resolutions), in order to help understand the SST sensitivity to multi-scale influences. This will allow us to assess the importance of explicitly simulated deep and shallow clouds, as well as the subgrid-scale effects (in coarse-model runs) upon SST variation. Results using both 1-D and 2-D OML models will be evaluated to assess the effects of horizontal advection.

Tao, W.-K.↗