SEARCH · Engineering Papers

Results for “thermodynamic coupling”

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.

SourceResearch connectionFrom yearThrough year

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

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

NASA NTRS1988

A note on chimney formation in ice edge regions

This paper investigates a coupled mesoscale ice-ocean system forced by winds and cooling, leading to the phenomenon of 'chimney' formation (i.e., the formation of a rapid deepening of the mixed layer due to a strong surface heat loss) in ice-edge regions. A numerical model, based on Hakkinen's (1987) coupled dynamic-thermodynamic model of an ice-ocean system, was used that included, in addition to wind stirring, the entrainment parameterization. The results are for the most part similar to the results of Hakkinen, indicating the overwhelming relevance of wind mixing, with penetrative convection having additional effects. The model calculations show that the area affected by ice edge upwelling is most prone to entrainment; the convectively induced entrainment is necessary for the production of a chimney.

Hakkinen, Sirpa↗

NASA NTRS1989

Interactions among radiation, convection, and large-scale dynamics in a general circulation model

The effects of radiatively active clouds on the climate simulated by a GCM are studied, focusing on the effects of upper tropospheric stratiform clouds associated with deep cumulus convection and the interactions of these clouds with convection and the large-scale circulation. The isolated effects of atmospheric cloud radiative forcing are analyzed using a model climate of an ocean-covered earth. It is found that the upper tropospheric cloudiness associated with cumulus convection has a powerful influence on both the convection and the large-scale circulation. Consideration is given to the importance of developing a GCM which includes a direct representation of the radiative, thermodynamical, and dynamical couplings between the convection parameterization and the large-scale saturation parameterization.

Randall, David A.↗

NASA NTRS1996

Surface Turbulent Fluxes Over Pack Ice Inferred from TOVS Observations

A one-dimensional, atmospheric boundary layer model is coupled to a thermodynamic ice model to estimate the surface turbulent fluxes over thick sea ice. The principal forcing parameters in this time-dependent model are the air temperature, humidity, and wind speed at a specified level (either at 2 m or at 850 mb) and the downwelling surface radiative fluxes. The free parameters. are the air temperature, humidity, and wind speed profiles below the specified level, the surface skin temperature, the ice temperature profile, and the surface turbulent fluxes. The goal is to determine how well we can estimate the turbulent surface heat and momentum fluxes using forcing parameters from atmospheric temperatures and radiative fluxes retrieved from the TIROS-N Operational Vertical Sounder (TOVS) data. Meteorological observations from the Lead Experiment (LeadEx, April 1992) ice camp are used to validate turbulent fluxes computed with the surface observations and the results are used to compare with estimates based on radio-sonde observations or with estimates based on TOVS data. We find that the TOVS-based estimates of the stress are significantly more accurate than those found with a constant geostrophic drag coefficient, with a root-mean-square error about half as large. This improvement is due to stratification effects included in the boundary layer model. The errors in the sensible heat flux estimates, however, are large compared to the small mean values observed during the field experiment.

Lindsay, R. W.↗

NASA NTRS2004

Ventless pressure control of two-phase propellant tanks in microgravity

This work studies pressurization and pressure control of a large liquid hydrogen storage tank. A finite element model is developed that couples a lumped thermodynamic formulation for the vapor region with a complete solution of the Navier-Stokes and energy equations for the flow and temperature fields in the liquid. Numerical results show that buoyancy effects are strong, even in microgravity, and can reposition a vapor bubble that is initially at the center of the tank to a region near the tank wall in a relatively short time. Long-term tank pressurization with the vapor bubble at the tank wall shows that after an initial transient lasting about a week, the final rate of pressure increase agrees with a purely thermodynamic analysis of the entire tank. However, the final pressure levels are quite different from thermodynamic predictions. Numerical results also show that there is significant thermal stratification in the liquid due to the effects of natural convection. A subcooled jet is used to provide simultaneous cooling and mixing in order to bring the tank pressure back down to its initial value. Three different jet speeds are examined. Although the lowest jet speed is ineffective at controlling the pressure because of insufficient penetration into the liquid region, the highest jet speed is shown to be quite effective at disrupting thermal stratification and reducing the tank pressure in reasonable time.

Weightlessness↗

NASA NTRS2006

Simulations of Sea-Ice Dynamics Using the Material-Point Method

In recent years, the availability of large volumes of recorded ice motion derived from high-resolution SAR data has provided an amazingly detailed look at the deformation of the ice cover. The deformation is dominated by the appearance of linear kinematic features that have been associated with the presence of leads. These remarkable data put us in a position to begin detailed evaluation of current coupled mechanical and thermodynamic models of sea ice. This presentation will describe the material point method (MPM) for solving these model equations. MPM is a numerical method for continuum mechanics that combines the best aspects of Lagrangian and Eulerian discretizations. The material points provide a Lagrangian description of the ice that models convection naturally. Thus, properties such as ice thickness and compactness are computed in a Lagrangian frame and do not suffer from errors associated with Eulerian advection schemes, such as artificial diffusion, dispersion, or oscillations near discontinuities. This desirable property is illustrated by solving transport of ice in uniform, rotational and convergent velocity fields. Moreover, the ice geometry is represented by unconnected material points rather than a grid. This representation facilitates modeling the large deformations observed in the Arctic, as well as localized deformation along leads, and admits a sharp representation of the ice edge. MPM also easily allows the use of any ice constitutive model. The versatility of MPM is demonstrated by using two constitutive models for simulations of wind-driven ice. The first model is a standard viscous-plastic model with two thickness categories. The MPM solution to the viscous-plastic model agrees with previously published results using finite elements. The second model is a new elastic-decohesive model that explicitly represents leads. The model includes a mechanism to initiate leads, and to predict their orientation and width. The elastic-decohesion model can provide similar overall deformation as the viscous-plastic model; however, explicit regions of opening and shear are predicted. Furthermore, the efficiency of MPM with the elastic-decohesive model is competitive with the current best methods for sea ice dynamics. Simulations will also be presented for an area of the Beaufort Sea, where predictions can be validated against satellite observations of the Arctic.

Sulsky, D.↗

NASA NTRS2013

A Model to Assess the Risk of Ice Accretion Due to Ice Crystal Ingestion in a Turbofan Engine and its Effects on Performance

The occurrence of ice accretion within commercial high bypass aircraft turbine engines has been reported under certain atmospheric conditions. Engine anomalies have taken place at high altitudes that were attributed to ice crystal ingestion, partially melting, and ice accretion on the compression system components. The result was one or more of the following anomalies: degraded engine performance, engine roll back, compressor surge and stall, and flameout of the combustor. The main focus of this research is the development of a computational tool that can estimate whether there is a risk of ice accretion by tracking key parameters through the compression system blade rows at all engine operating points within the flight trajectory. The tool has an engine system thermodynamic cycle code, coupled with a compressor flow analysis code, and an ice particle melt code that has the capability of determining the rate of sublimation, melting, and evaporation through the compressor blade rows. Assumptions are made to predict the complex physics involved in engine icing. Specifically, the code does not directly estimate ice accretion and does not have models for particle breakup or erosion. Two key parameters have been suggested as conditions that must be met at the same location for ice accretion to occur: the local wet-bulb temperature to be near freezing or below and the local melt ratio must be above 10%. These parameters were deduced from analyzing laboratory icing test data and are the criteria used to predict the possibility of ice accretion within an engine including the specific blade row where it could occur. Once the possibility of accretion is determined from these parameters, the degree of blockage due to ice accretion on the local stator vane can be estimated from an empirical model of ice growth rate and time spent at that operating point in the flight trajectory. The computational tool can be used to assess specific turbine engines to their susceptibility to ice accretion in an ice crystal environment.

Jorgenson, Philip C. E.↗

NASA NTRSDate not supplied

Development of an Additive Manufactured Cryogenic TVS Augmented Injector

Cryogenic fluid transfer is critical in ground and space systems. Transfer typically requires venting a receiver vessel during the chill and fill process to maintain a pressure favorable to fluid flow. Venting in space is problematic due to low-gravity where fluid position is indeterminate; possibly allowing fluid to flow out the vent. To control receiver tank venting the vented chill/no-vent fill methodology is attractive. An additively manufactured Injector coupled with a Thermodynamic Venting System (TVS) augmentation can enable the vented chill/no-vent fill approach with simplified operations. The TVS augmented injector includes an internal heat exchanger, fluid injector spray head, and external surface condensation heat exchange; combined with integrated flow paths in one part. Operations using the injector allow the tank vent to be closed very early in the transfer and not reopened. The TVS circuit eliminates the potential for stalled fill by maintaining cooling flow on the injector and condenser elements, reducing tank pressure. The injector can also be employed for long-term cryo fluid storage as a means of controlling tank pressure if integrated with a recirculation pump. A cryo-cooler can be used in place of/or augmenting the TVS flow circuit for a zero loss pressure control system. The Injector was developed using the principles of design for additive manufacture and printed as one part, resulting in part count reduction, improved reproducibility, shorter lead times, and reduced cost. The TVS augmented injector is compact and highly adaptable simplifying its integration. A key benefit is that this methodology enables simplified operations in reaching 100% fill conditions.

Omar Mireles↗

NASA NTRS2017

Approach to Modeling Boundary Layer Ingestion Using a Fully Coupled Propulsion-RANS Model

Airframe-propulsion integration concepts that use boundary layer ingestion have the potential to reduce aircraft fuel burn. One concept that has been recently explored is NASA's Starc-ABL aircraft configuration, which offers the potential for 12% mission fuel burn reduction by using a turbo-electric propulsion system with an aft-mounted electrically driven boundary layer ingestion propulsor. This large potential for improved performance motivates a more detailed study of the boundary layer ingestion propulsor design, but to date, analyses of boundary layer ingestion have used uncoupled methods. These methods account for only aerodynamic effects on the propulsion system or propulsion system effects on the aerodynamics, but not both simultaneously. This work presents a new approach for building fully coupled propulsive-aerodynamic models of boundary layer ingestion propulsion systems. A 1D thermodynamic cycle analysis is coupled to a RANS simulation to model the Starc-ABL aft propulsor at a cruise condition and the effects variation in propulsor design on performance are examined. The results indicates that both propulsion and aerodynamic effects contribute equally toward the overall performance and that the fully coupled model yields substantially different results compared to uncoupled. The most significant finding is that boundary layer ingestion, while offering substantial fuel burn savings, introduces throttle dependent aerodynamics effects that need to be accounted for. This work represents a first step toward the multidisciplinary design optimization of boundary layer ingestion propulsion systems.

Boundary Layer Ingestion↗

DOE OSTI2025

Anticorrelation between electron-phonon coupling strength and stability of ternary metal diborides

This paper endeavors to expand the exploration of prospective superconductivity within the 𝑃⁢6/𝑚⁢𝑚⁢𝑚−MgB 2 structural framework by substituting Mg atoms with a random selection of two distinct cations from a pool of 24 candidates. Employing a high-throughput screening methodology to scrutinize zone-center electron-phonon couplings, we found an inverse relationship between thermodynamic stability and electron-phonon coupling strength in the resulting ternary metal borides. Notably, within the MgNiB 4 system, we demonstrated that electron-phonon superconductivity is stronger in less dynamically stable configurations. Further, the qualitative reasons for the emerging instabilities when electron-phonon coupling becomes strong are discussed. This insight indicates that the quest for superior superconducting materials should concentrate on the intricate balance between electronic and phononic characteristics in structures that are on the brink of instability.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS2004

Thermodynamic Modeling of a Solid Oxide Fuel Cell to Couple with an Existing Gas Turbine Engine Model

The Controls and Dynamics Technology Branch at NASA Glenn Research Center are interested in combining a solid oxide fuel cell (SOFC) to operate in conjunction with a gas turbine engine. A detailed engine model currently exists in the Matlab/Simulink environment. The idea is to incorporate a SOFC model within the turbine engine simulation and observe the hybrid system's performance. The fuel cell will be heated to its appropriate operating condition by the engine s combustor. Once the fuel cell is operating at its steady-state temperature, the gas burner will back down slowly until the engine is fully operating on the hot gases exhausted from the SOFC. The SOFC code is based on a steady-state model developed by The U.S. Department of Energy (DOE). In its current form, the DOE SOFC model exists in Microsoft Excel and uses Visual Basics to create an I-V (current-voltage) profile. For the project's application, the main issue with this model is that the gas path flow and fuel flow temperatures are used as input parameters instead of outputs. The objective is to create a SOFC model based on the DOE model that inputs the fuel cells flow rates and outputs temperature of the flow streams; therefore, creating a temperature profile as a function of fuel flow rate. This will be done by applying the First Law of Thermodynamics for a flow system to the fuel cell. Validation of this model will be done in two procedures. First, for a given flow rate the exit stream temperature will be calculated and compared to DOE SOFC temperature as a point comparison. Next, an I-V curve and temperature curve will be generated where the I-V curve will be compared with the DOE SOFC I-V curve. Matching I-V curves will suggest validation of the temperature curve because voltage is a function of temperature. Once the temperature profile is created and validated, the model will then be placed into the turbine engine simulation for system analysis.

Brinson, Thomas E.↗

NASA NTRS2012

Simulation of Oxygen Disintegration and Mixing With Hydrogen or Helium at Supercritical Pressure

The simulation of high-pressure turbulent flows, where the pressure, p, is larger than the critical value, p(sub c), for the species under consideration, is relevant to a wide array of propulsion systems, e.g. gas turbine, diesel, and liquid rocket engines. Most turbulence models, however, have been developed for atmospheric-p turbulent flows. The difference between atmospheric-p and supercritical-p turbulence is that, in the former situation, the coupling between dynamics and thermodynamics is moderate to negligible, but for the latter it is very significant, and can dominate the flow characteristics. The reason for this stems from the mathematical form of the equation of state (EOS), which is the perfect-gas EOS in the former case, and the real-gas EOS in the latter case. For flows at supercritical pressure, p, the large eddy simulation (LES) equations consist of the differential conservation equations coupled with a real-gas EOS. The equations use transport properties that depend on the thermodynamic variables. Compared to previous LES models, the differential equations contain not only the subgrid scale (SGS) fluxes, but also new SGS terms, each denoted as a correction. These additional terms, typically assumed null for atmospheric pressure flows, stem from filtering the differential governing equations, and represent differences between a filtered term and the same term computed as a function of the filtered flow field. In particular, the energy equation contains a heat-flux correction (q-correction) that is the difference between the filtered divergence of the heat flux and the divergence of the heat flux computed as a function of the filtered flow field. In a previous study, there was only partial success in modeling the q-correction term, but in this innovation, success has been achieved by using a different modeling approach. This analysis, based on a temporal mixing layer Direct Numerical Simulation database, shows that the focus in modeling the q-correction should be on reconstructing the primitive variable gradients rather than their coefficients, and proposes the approximate deconvolution model (ADM) as an effective means of flow field reconstruction for LES heat flux calculation. Further, results for a study conducted for temporal mixing layers initially containing oxygen in the lower stream, and hydrogen or helium in the upper stream, show that, for any LES, including SGS-flux models (constant-coefficient Gradient or Scale-Similarity models, dynamic-coefficient Smagorinsky/Yoshizawa or mixed Smagorinsky/Yoshizawa/Gradient models), the inclusion of the q-correction in the LES leads to the theoretical maximum reduction of the SGS heat-flux difference. The remaining error in modeling this new subgrid term is thus irreducible.

Bellan, Josette↗

DOE OSTI2025

Single-Electron Transfer Stabilizes Metastable Alane in a Bipyridine-Functionalized MOF Nanopore

Nanoconfinement of metastable hydrides such as alane (AlH 3 ) in functionalized porous frameworks represents a promising yet largely untapped strategy for high-capacity energy storage. Despite its potential, the underlying mechanisms responsible for the thermodynamic stabilization of metastable hydrides are poorly understood. Here, in this study, concepts from solution Lewis acid–base chemistry were applied to a crystalline metal–organic framework (MOF) to stabilize AlH 3 . The long-range order and synthetically versatile pore chemistry enabled us to reveal the intimate details of the hydride-host interactions. Powder X-ray diffraction (PXRD) of AlH 3 -infiltrated UiO-67bpy (Zr 6 O 4 (OH) 4 (bpydc) 6 ; bpydc 2– = 2,2′-bipyridine-5,5′-dicarboxylate) confirms that the framework maintains its crystallinity up to 250 °C and is stable under 70 MPa H 2 pressure. We find that thermodynamic stabilization arises from coupled single-electron and hydrogen-atom transfer from AlH 3 to bipyridine-containing linkers. Electron paramagnetic resonance detects a bipyridyl radical with an anisotropic g-tensor (g values of 2.00132, 2.00215, and 2.00275), consistent with Bader charge analysis predicting 0.82 e – transferred from the hydride to the MOF. These findings establish critical structure–property relationships underpinning AlH 3 stabilization and suggest general strategies for tuning the thermodynamics and kinetics of metastable materials.

Shivanna, Mohana [Sandia National Laboratories (SN↗

DOE OSTI2025

Inverse Thermodynamics: Designing Interactions for Targeted Phase Behavior

The traditional goal of inverse self-assembly is to design interactions that drive particles toward a desired target structure. However, achieving successful self-assembly also requires tuning the thermodynamic conditions under which the structure is stable. In this work, we extend the inverse design paradigm to explicitly address this challenge by developing a framework for inverse thermodynamics, i.e., the design of interaction potentials that realize specific thermodynamic behavior. As a step in this direction, using patchy particle mixtures as a model system, we demonstrate how precise control over both bonding topology and bond energetics enables the programming of targeted phase behavior. In particular, we establish design principles for azeotropic demixing and show how to create mixtures that exhibit azeotropy at any prescribed composition. Our predictions are validated through Gibbs-ensemble simulations [Panagiotopoulos, Mol. Phys. 1987, 61, 813−826]. These results highlight the necessity of coupling structural design with thermodynamic engineering, and provide a blueprint for controlling complex phase behavior in multicomponent systems.

Azeotropes↗

DOE OSTI2024

Thermodynamic stability of a spin microemulsion in Rashba spin-orbit-coupled bosons

Recent finite-temperature numerical simulations have unveiled a quantum “spin” microemulsion analog, found by raising the temperature of a stripe supersolid phase in a Rashba spin-orbit coupled Bose gas. This microemulsion state is a highly correlated, isotropic normal fluid where atoms self-arrange based on their internal pseudospin into patterns that resemble bicontinuous microemulsions. This finding leaves several open questions regarding the broader accessibility of this phase in experiments. Here, we use equilibrium finite-temperature numerical simulations based on a coherent-state path integral representation to perform a computational investigation into the thermodynamic stability of the spin microemulsion state. Numerical simulations emphasize the requirement of a nearly, but not perfectly, isotropic spin-orbit coupling in order to achieve the microemulsion phase in cold-atom experiments. Moreover, the microemulsion state exists independent of miscibility of the pseudospin components and for a wide range of pseudospin population imbalance, suggesting a high degree of flexibility in choosing the atom and hyperfine states in an experimental realization. Lastly, we demonstrate this feasibility by mimicking a Rashba spin-orbit-coupled 87 Rb experiment in an isotropic harmonic trap, where we confirm the microemulsion's existence via its density profile and equilibrium quasimomentum distribution.

Complex Langevin dynamics↗

DOE OSTI2026

Wind and Temperature Consensus at Horn Point, HU-Beltsville, Piney Run (Maryland) in support of CoURAGE

The Maryland Department of the Environment (MDE) operates a ground-based atmospheric profiling network consisting of collocated radar wind profilers (RWP) and radio acoustic sounding systems (RASS) as part of its Ambient Air Monitoring Program. This network provides continuous observations of wind and temperature structure in the lower troposphere to support air quality forecasting, regulatory analysis, and atmospheric research. The network currently includes three fixed sites across Maryland: Horn Point (HP, lower eastern shore) [38.587525°,-76.141006°], Howard University-Beltsville (HUB, central Maryland) [39.055277°, -76.878632°], and Piney Run (PR, western Maryland) [39.705950°, -79.012000°] The network is designed to capture regional variability in atmospheric transport and boundary-layer processes. These systems measure vertical profiles of horizontal wind speed and direction using Doppler radar techniques, with observations typically spanning from ~100 m above ground level up to approximately 2.5–4 km. Measurements are derived from the Doppler shift of backscattered electromagnetic signals, enabling retrieval of wind vectors at multiple altitudes with high temporal resolution (e.g., 30-minute averages reported every 6 minutes). Each radar wind profiler is paired with a Radio Acoustic Sounding System (RASS) to provide profiles of virtual temperature in the lower atmosphere (~100–200 m AGL) by measuring the propagation speed of acoustic waves. Together, the RWP/RASS system yields a coupled data set of thermodynamic and kinematic atmospheric structure, including additional parameters such as vertical velocity, radial velocity, signal-to-noise ratio, and spectral width for advanced analysis. There are two types of files for each station: wind data (files with a "w" prefix) and virtual temperature RASS data (files with a "t" prefix). The wind data files are in the format wYYDDD.cns, where YY is the 2-digit year and DDD is the day of the year. The RASS virtual temperature data files are in the format tYYDDD.cns. Each record has the following header structure: Line 1 : Station Name RASS files Line 2 : RASS rev DeTect_2.0, WINDS files Line 2 : WINDS rev ATI 5.1 Line 3 : N latitude, W longitude, and site elevation (m) Line 4 : Date and begin time of consensus: yy mm dd hh mn ss plus # minutes to add to get UTC Line 5 : Consensus averaging time (minutes); number of beams; number of range gates Line 6 : Number of records required to make consensus (num) total number of records (tot) and the consensus window size (m/s) in the format: num:tot (window) RASS files Line 7 : no. of coded cells, no. of spec, pulse width (ns), and inter-pulse period (µs), WINDS files Line 7 : No. of coded cells, no. of spectra, pulse width (ns), and inter-pulse period (µs), each with a pair of values: first value is for oblique beams, second for vertical RASS files Line 8 : Full scale Doppler value (m/s) Delay to first gate (ns) Number of gates Spacing of gates (ns), WINDS files Line 8 : Full scale Doppler velocity (m/s), oblique and vertical Vertical correction applied to oblique beams? (0 = no, 1 = yes) Delay to first gate (ns), oblique and vertical Number of gates, oblique and vertical Spacing of gates (ns), oblique and vertical Line 9 : Azimuth and elevation (9s indicate vertical beam not used) RASS files Line 10, values : HT = Height above ground (km), T = Uncorrected virtual temperature consensus (deg C), Tc = Corrected virtual temperature consensus (deg C), W = Vertical wind consensus (9s indicate vertical beam not used, w-component, positive upward, m/s), CNT = Number of records that made consensus (for the 3 values in same order), SNR = Average signal to noise ratio (dB) of records in consensus (same order) WINDS files Line 10, values : HT = Height above ground (km), SPD = Wind speed (m/s), DIR = Wind direction (deg E of N from N), RAD = Radial velocities for each beam (m/s) in order given in azimuth and elevation line (positive toward radar; 9s indicate vertical beam not used, CNT = Number of records that made consensus, SNR = Average signal to noise ratio (dB) of records in consensus

{"wind speed and direction",temperature}↗

NASA NTRS1987

The 1982-1983 El Nino Atlas: Nimbus-7 microwave radiometer data

Monthly maps of sea surface temperature, atmospheric water vapor, and surface level wind speed as measured by the Scanning Multichannel Microwave Radiometer (SMMR) on the Nimbus-7 satellite for the tropical Pacific from June 1982 to October 1983, during one of the most intense El Nino Southern Oscillations (ENSO) episodes, are presented. The non-ENSO annual cycle was compiled by averaging the 1980 and 1981 data for each calendar month and was removed from monthly fields of 1982 and 1983 to reveal the anomalous distributions. The anomaly fields and part of the non-ENSO annual cycle are also presented. This study and earlier evaluations demonstrate that the Nimbus/SMMR can be used to monitor large scale and low frequency variabilities in the tropical ocean. The SMMR data support and extend conventional measurements. The variabilities of the three parameters are found to represent various aspects of ENSO related through ocean atmosphere interaction. Their simultaneous and quantitative descriptions pave the way for the derivation of ocean atmosphere latent heat exchange and further the understanding of the coupled atmospheric and oceanic thermodynamics.

Liu, W. Timothy↗

NASA NTRS1987

Directional solidification in an imperfect furnace

The upward directional solidification of a binary alloy in which the rejected solute is more dense than the solvent is considered. The system is stably stratified both thermally and solutally. If heat is lost through the sidewalls, horizontal thermal gradients drive a steady double-diffusive convection. Solidification cells with small aspect ratios (width/length) A are considered, and asymptotic solutions to the steady double-diffusive system are obtained as A tends to 0 under the assumption heat losses are also small. When the thermal Grashof number is large, buoyancy boundary layers are at the sidewalls. The shape of the crystal-melt interface is determined by the coupled effects of local thermodynamic equilibrium, diffusion and convection.

Brattkus, K.↗

NASA NTRS2003

Towards an Automated Full-Turbofan Engine Numerical Simulation

The objective of this study was to demonstrate the high-fidelity numerical simulation of a modern high-bypass turbofan engine. The simulation utilizes the Numerical Propulsion System Simulation (NPSS) thermodynamic cycle modeling system coupled to a high-fidelity full-engine model represented by a set of coupled three-dimensional computational fluid dynamic (CFD) component models. Boundary conditions from the balanced, steady-state cycle model are used to define component boundary conditions in the full-engine model. Operating characteristics of the three-dimensional component models are integrated into the cycle model via partial performance maps generated automatically from the CFD flow solutions using one-dimensional meanline turbomachinery programs. This paper reports on the progress made towards the full-engine simulation of the GE90-94B engine, highlighting the generation of the high-pressure compressor partial performance map. The ongoing work will provide a system to evaluate the steady and unsteady aerodynamic and mechanical interactions between engine components at design and off-design operating conditions.

Reed, John A.↗