Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model initialization”

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 1,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
NASA NTRS2015

AST Launch Vehicle Acoustics

The liftoff phase induces acoustic loading over a broad frequency range for a launch vehicle. These external acoustic environments are then used in the prediction of internal vibration responses of the vehicle and components which result in the qualification levels. Thus, predicting these liftoff acoustic (LOA) environments is critical to the design requirements of any launch vehicle. If there is a significant amount of uncertainty in the predictions or if acoustic mitigation options must be implemented, a subscale acoustic test is a feasible pre-launch test option to verify the LOA environments. The NASA Space Launch System (SLS) program initiated the Scale Model Acoustic Test (SMAT) to verify the predicted SLS LOA environments and to determine the acoustic reduction with an above deck water sound suppression system. The SMAT was conducted at Marshall Space Flight Center and the test article included a 5% scale SLS vehicle model, tower and Mobile Launcher. Acoustic and pressure data were measured by approximately 250 instruments. The SMAT liftoff acoustic results are presented, findings are discussed and a comparison is shown to the Ares I Scale Model Acoustic Test (ASMAT) results.

Houston, Janice↗
NASA NTRS2015

Application of the Refined Zigzag Theory to the Modeling of Delaminations in Laminated Composites

The Refined Zigzag Theory is applied to the modeling of delaminations in laminated composites. The commonly used cohesive zone approach is adapted for use within a continuum mechanics model, and then used to predict the onset and propagation of delamination in five cross-ply composite beams. The resin-rich area between individual composite plies is modeled explicitly using thin, discrete layers with isotropic material properties. A damage model is applied to these resin-rich layers to enable tracking of delamination propagation. The displacement jump across the damaged interfacial resin layer is captured using the zigzag function of the Refined Zigzag Theory. The overall model predicts the initiation of delamination to within 8% compared to experimental results and the load drop after propagation is represented accurately.

Groh, Rainer M. J.↗
NASA NTRS2016

A New Tropical Cyclone Dynamic Initialization Technique Using High Temporal and Spatial Density Atmospheric Motion Vectors and Airborne Field Campaign Data

Background: Initialization of tropical cyclones in numerical weather prediction (NWP) systems is a great challenge: Mass-wind field balance; Secondary circulation and heating; Asymmetries. There can be large adjustments in structure and intensity in the first 24 hours if the initial vortex is not in balance: Spurious gravity waves; Spin-up (model and physics). Existing mesoscale NWP model TC (Tropical Cyclone) initialization strategies: Bogus vortex, cold start from global analyses; 3DVAR or 4DVAR, possibly with synthetic observations; EnKF (Ensemble Kalman Filter); Dynamic initialization. Dynamic initialization allows vortex to have improved balance and physics spin-up at the initial time (e.g., Hendricks et al. 2013, 2011; Nguyen and Chen 2011; Fiorino and Warner 1981; Hoke and Anthes 1976). Himawari-8 geostationary satellite has capability of continuous imagery (10-minutes) over the full disk: New GOES-R satellites will have same capability. This will allow for unprecedented observations of tropical cyclones. However, current data assimila1on systems are not capable of ingesting such high temporal observations (Atmospheric Mo1on Vectors - AMVs). Hourly AMVs are produced, and thinned to 100-kilometer spacing in the horizontal. An entirely new data assimilation concept is required to utilize these observations.

remote sensing↗
NASA NTRS2017

Environmental Barrier Coating Fracture, Fatigue and High-Heat-Flux Durability Modeling and Stochastic Progressive Damage Simulation

Advanced environmental barrier coatings will play an increasingly important role in future gas turbine engines because of their ability to protect emerging light-weight SiC/SiC ceramic matrix composite (CMC) engine components, further raising engine operating temperatures and performance. Because the environmental barrier coating systems are critical to the performance, reliability and durability of these hot-section ceramic engine components, a prime-reliant coating system along with established life design methodology are required for the hot-section ceramic component insertion into engine service. In this paper, we have first summarized some observations of high temperature, high-heat-flux environmental degradation and failure mechanisms of environmental barrier coating systems in laboratory simulated engine environment tests. In particular, the coating surface cracking morphologies and associated subsequent delamination mechanisms under the engine level high-heat-flux, combustion steam, and mechanical creep and fatigue loading conditions will be discussed. The EBC compostion and archtechture improvements based on advanced high heat flux environmental testing, and the modeling advances based on the integrated Finite Element Analysis Micromechanics Analysis Code/Ceramics Analysis and Reliability Evaluation of Structures (FEAMAC/CARES) program will also be highlighted. The stochastic progressive damage simulation successfully predicts mud flat damage pattern in EBCs on coated 3-D specimens, and a 2-D model of through-the-thickness cross-section. A 2-parameter Weibull distribution was assumed in characterizing the coating layer stochastic strength response and the formation of damage was therefore modeled. The damage initiation and coalescence into progressively smaller mudflat crack cells was demonstrated. A coating life prediction framework may be realized by examining the surface crack initiation and delamination propagation in conjunction with environmental degradation under high-heat-flux and environment load test conditions.

Mechanical creep and fatigue loading condition↗
NASA NTRS2018

Earth-Based Analogs & Modeling for Exercise Biomechanics in Space

The biomechanics of exercise in space is difficult to study and there are unknowns surrounding exercise performance on future space exploration countermeasures systems. These issues are beginning to be addressed through enhanced modeling techniques fueled initially by human-in-the-loop data collections in ground-based environments. The presentation will focus on an effort completed at the University of South Florida to apply the Computer Assisted Rehabilitation Environment (CAREN) system to address a human spaceflight need. The research explored the interaction between a human and a moving platform while exercise was completed.

Kaitlin Lostroscio↗
NASA NTRSDate not supplied

Aerodynamic and Acoustic Interactions Associated with Inboard Propeller-Wing Configurations

A series of aerodynamic performance and acoustic measurements have been made on a range of inboard propeller-wing interaction configurations in the NASA Langley Low Speed Aeroacoustic Wind Tunnel (LSAWT). The results presented in this paper are part of a more expansive testing campaign encompassing both single propeller-wing and multipropeller-wing interactions, the former of which is discussed in the present work. The primary testing parameters of interest to this study are the axial and vertical positioning of the wing relative to the propeller slipstream under a constant propeller advance ratio. A multi-faceted computational effort was also employed in an effort to identify reflection and scattering effects imposed by both the wing geometry as well as the primary components of the facility test setup. This effort consisted of aerodynamic predictions using high-fidelity computational fluid dynamics (CFD), acoustic predictions using an impermeable Ffowcs Williams and Hawkings (FW-H) solver, and acoustic scattering predictions. Acoustic measurements reveal variations in the acoustic directivity behavior of the propeller blade passage frequency for even modest variations in wing position. CFD-based acoustic predictions reveal discrepancies relative to the experimental data, which is believed to be due to complex acoustic scattering behavior within the test section. Initial attempts at modeling the scattered acoustic field showed functional dependency of the acoustic amplitude variations on the wing position relative to the propeller disk, however discrepancies with experimental data remain.

Nikolas S. Zawodny↗
NASA NTRS2020

Optimized Profile Retrievals of Aerosol Microphysical Properties from Simulated Spaceborne Multiwavelength Lidar

This work is an expanded study of one previously published onretrievals of aerosol microphysical properties from space-borne multiwavelengthlidarmeasurements. The earlier studiesand this one weredone in the framework of the NASA Aerosol-Clouds-Ecosystems (now the Aerosol Clouds Convection and Precipitation) NASA mission. The focus here is on the capabilities of a simulated spacebornemultiwavelengthlidar system for retrieving aerosol complex refractive index (m = mr+ imi) and spectral single scattering albedo (SSA(λ)), although other bulk parameters such as effective (reff) radius and particle volume (V) and surface (S) concentrations are also studied. The novelty presented here is the use of recently published, case dependent optimized-constraints on the microphysical retrievals using three backscattering coefficients (β) at 355, 532 and 1064 nm and two extinction coefficients (α) at 355 and 532 nm, typically known as the stand-alone 3β+2α lidar inversion. Case-dependent optimized-constraints (CDOC) limit the ranges of refractive index, both real (mr) and imaginary (mi) parts, and of radii that are permitted in the retrievals. Such constraints are selected directly from the 3β+2α41measurements through an analysis of the relationship between spectral dependence of aerosol extinction-to-backscatter ratios (LR) and the Ångström exponent of extinction. The analyses presented here for different sets ofsize distributions and refractive indices reveal that the direct determination of CDOCareonly feasible for cases where the uncertaintiesin the input optical data areless than 15 %.Forthe same simulated spacebornesystem and yield than in Whiteman et al., (2018), we demonstrated that the use of CDOC as essential for the retrievals of refractive index and also largely improved retrieval of bulk parameters. A discussion of the global representativeness of CDOC is presented using simulated lidar data from a 24-hour satellite track using GEOS model output to initialize the lidar simulator.We found that CDOCare representative of many aerosol mixtures in spite of some outliers (e.g. highly hydrated particles) associatedwith the assumptions of bimodal size distributions and of the same refractive index for fine and coarse modes. Moreover, sensitivity tests performed using synthetic data reveal that retrievals of imaginary refractive index (mi) and SSA are extremely sensitive to β(355).

NASA Aerosol-Clouds-Ecosystems↗
NASA NTRS2018

Is MBSE Helping? Measuring Value on Europa Clipper

JPL’s Integrated Model Centric Engineering Initiative (IMCE) has led the infusion of MBSE at JPL. In 2009 I authored for IMCE a set of challenges that confronted systems engineering at JPL, based on examination of issues and lessons from several recent projects (the “Five System Engineering Challenges”). In 2012 I augmented this work to describe specific areas where it appeared MBSE could help address these challenges. Around the same time, the Europa Clipper Mission became the first project at JPL to attempt, in close collaboration with IMCE, a widespread adoption of MBSE. This adoption has been highly successful so far. Europa Clipper completed Phase A (Formulation) with a successful System Requirements Review / Mission Definition Review, and is nearing the end of Phase B, with a successful Flight System Preliminary Design Review in October 2017. This is an opportune time to take stock of the progress and the work yet to go. This paper proposes a practical method for measuring value added through use of MBSE, applies this method to the Europa Clipper Project, and suggests how it might be applied to other projects and organizations.

Bayer, Todd↗
NASA NTRSDate not supplied

Pragmatic Stress Prediction on Additively Manufactured Coupons

Prediction of residual stresses from process parameters for additively manufactured large metal parts is computationally expensive. NASA is currently developing meter-scale parts with direct energy deposition. Practically, the predictive computational methods need to efficiently scale-up to meter-scale parts. Coupled thermal-mechanical multi-physics simulations have been developed with the pragmatic method using ABAQUS, COMSOL Multiphysics, ALE3D software. The residual stresses are a result of the manufacturing process which creates thermal cycling of the build layers. The pragmatic method uses lumped thermal layers for stress predictions to reduce computational costs. The stress predictions as well as deformations of the different codes are compared with each other and with ANSYS Additive using identical material models, boundary and initial conditions. The codes were used to simulate three different geometries: a thin wall, hollow cylinder and twin-cantilever part. The coupon parts were then manufactured with Inconel-625. The residual stresses in these parts were measured using X-ray diffraction as well as neutron beam diffraction at NIST. The stress measurements for the two technologies are compared. The pragmatic stress prediction method enabled predictions of the multi-centimeter scale parts using desktop computer workstations in only a few hours for each coupon. The results of the simulated stress predictions compared favorably with the measured stresses even though thermally lumped layers were employed. Finally, a two-meter scale nozzle was simulated using ANSYS Additive. The simulations were used to examine the build orientation trade-space with respect to resulting geometric deformation. The predicted deformations were compared to measurements of an actual subscale part manufactured with direct energy deposition.

Additive Manufacturing↗
NASA NTRSDate not supplied

Pragmatic Stress Prediction on Additively Manufactured Coupons

Prediction of residual stresses from process parameters for additively manufactured large metal parts is computationally expensive. NASA is currently developing meter-scale parts with direct energy deposition. Practically, the predictive computational methods need to efficiently scale-up to meter-scale parts. Coupled thermal-mechanical multi-physics simulations have been developed with the pragmatic method using ABAQUS, COMSOL Multiphysics, ALE3D software. The residual stresses are a result of the manufacturing process which creates thermal cycling of the build layers. The pragmatic method uses lumped thermal layers for stress predictions to reduce computational costs. The stress predictions as well as deformations of the different codes are compared with each other and with ANSYS Additive using identical material models, boundary and initial conditions. The codes were used to simulate three different geometries: a thin wall, hollow cylinder and twin-cantilever part. The coupon parts were then manufactured with Inconel-625. The residual stresses in these parts were measured using X-ray diffraction as well as neutron beam diffraction at NIST. The stress measurements for the two technologies are compared. The pragmatic stress prediction method enabled predictions of the multi-centimeter scale parts using desktop computer workstations in only a few hours for each coupon. The results of the simulated stress predictions compared favorably with the measured stresses even though thermally lumped layers were employed. Finally, a two-meter scale nozzle was simulated using ANSYS Additive. The simulations were used to examine the build orientation trade-space with respect to resulting geometric deformation. The predicted deformations were compared to measurements of an actual subscale part manufactured with direct energy deposition.

pragmatic method↗
NASA NTRSDate not supplied

Conceptual Design of Propulsors for the SUSAN Electrofan Transport Aircraft

Conceptual designs of the propulsor modules for the SUSAN electro-fan aircraft are sought after. Shaft power requirement is reduced by using boundary layer ingesting propulsion technology. There are several obstacles to designing feasible BLI propulsion systems such as inlet distortion, flow blockage from ingesting low-momentum flow, strong coupling between airframe and propulsion systems. Thus, the high fidelity CFD tool is indispensable to assess the performance of the propulsion systems and evaluate the inlet profiles during the design process. Consequently, the CFD data is used to update the inputs for the NPSS model from the initial stage of the system design. Hence, the inlet sizing, fan design, and estimation of the power saving are carried out. Various installation concepts of the mail-slot nacelle, such as under-/over-wing and trailing edge configurations, are investigated for the wing-mounted turbo-electric distributed propulsor module, and their power-saving is evaluated. As for the tail-mounted turbofan engine, the fan diameter and inlet captured area is determined based on the CFD profiles. An appropriate area ratio of the bypass and core ducts is derived from meeting the target bypass ratio from the system design. The baseline fan stage is analyzed by full annulus URANS CFD to assess the efficiency penalty due to the ingested pressure and swirl distortions.

SUSAN Electro-Fan Aircraft↗
NASA NTRSDate not supplied

Conceptual Design of Propulsors for the SUSAN Electrofan Transport Aircrafts

Conceptual designs of the propulsor modules for the SUSAN electro-fan aircraft are sought after. Shaft power requirement is reduced by using boundary layer ingesting propulsion technology. There are several obstacles to designing feasible BLI propulsion systems such as inlet distortion, flow blockage from ingesting low-momentum flow, strong coupling between airframe and propulsion systems. Thus, the high fidelity CFD tool is indispensable to assess the performance of the propulsion systems and evaluate the inlet profiles during the design process. Consequently, the CFD data is used to update the inputs for the NPSS model from the initial stage of the system design. Hence, the inlet sizing, fan design, and estimation of the power saving are carried out. Various installation concepts of the mail-slot nacelle, such as under-/over-wing and trailing edge configurations, are investigated for the wing-mounted turbo-electric distributed propulsor module, and their power-saving is evaluated. As for the tail-mounted turbofan engine, the fan diameter and inlet captured area is determined based on the CFD profiles. An appropriate area ratio of the bypass and core ducts is derived from meeting the target bypass ratio from the system design. The baseline fan stage is analyzed by full annulus URANS CFD to assess the efficiency penalty due to the ingested pressure and swirl distortions.

Electric Aircraft↗
NASA NTRS2023

Seasonal Forecasting Skill for the High Mountain Asia Region in the Goddard Earth Observing System

Seasonal variability of the global hydrologic cycle directly impacts human activities, including hazard assessment and mitigation, agricultural decisions, and water resources management. This is particularly true across the High Mountain Asia (HMA) region, where availability of water resources can change depending on local seasonality of the hydrologic cycle. Forecasting the atmospheric states and surface conditions, including hydrometeorological relevant variables, at subseasonal-to-seasonal (S2S) lead times of weeks-to-months is an area of active research and development. NASA’s 15 Goddard Earth Observing System (GEOS) S2S prediction system has been developed with this research goal in mind. Here, we benchmark the forecast skill of GEOS-S2S (version 2) hydrometeorological forecasts at 1-3 month lead times in the HMA region, including a portion of the Indian Subcontinent, during the retrospective forecast period, 1981-2016. To assess forecast skill, we evaluate 2-m air temperature, total precipitation, fractional snow cover, snow water equivalent, surface soil moisture, and terrestrial water storage forecasts against the Modern-Era Retrospective analysis for Research and 20 Applications, Version 2 (MERRA-2) and independent reanalysis data, satellite observations, and data fusion products. Anomaly correlation is highest when the forecasts are evaluated against MERRA-2 and particularly in variables with long memory in the climate system, likely due to similar initial conditions and model architecture used in GEOS-S2S and MERRA-2. When compared to MERRA-2, results for the 1-month forecast skill range from anomaly correlation of R anom =0.18 for precipitation to R anom =0.62 for soil moisture. Anomaly correlations are consistently lower when forecasts are 25 evaluated against independent observations; results for the 1-month forecast skill range from R anom =0.13 for snow water equivalent to R anom =0.24 for fractional snow cover. We find that, generally, hydrometeorological forecast skill is dependent on the forecast lead time, the memory of the variable within the physical system, and the validation dataset used. Overall, these results benchmark the GEOS-S2S system’s ability to forecast HMA hydrometeorology.

GEOS↗
NASA NTRSDate not supplied

Expanding the Time Dimension of Hyperspectral Infrared Sounding Observations: Designing the NUCAPS-Forecast System

Low Earth orbit (LEO) hyper-spectral infrared (IR) sounders have significant yet untapped potential for characterizing thermodynamic environments of convective initiation and ongoing convection. While LEO soundings are of enormous value to weather forecasters, they cannot provide the temporal resolution needed to resolve the rapidly evolving thermodynamics of the convective environment. A novel methodology was developed to extend snapshots of LEO soundings forward in time up to six hours into the future and has been made available within the NOAA National Weather Service Advanced Weather Interactive Processing System (AWIPS) for testing in a quasi-operational setting. This methodology is based on parcel forward-trajectory calculations from the satellite observing time to recreate future soundings of temperature and moisture at regularly gridded intervals in space and time. The soundings are based on NOAA-Unique Combined Atmospheric Processing System (NUCAPS) retrievals from the Suomi NPP and NOAA-20 satellite platforms to allow for rapid near-real time product processing. Science code was transitioned to operational production in 2019 and the product, NUCAPS-Forecast, was tested at the NOAA Hazardous Weather Testbed Spring Experiment in 2019 and 2021. Forecaster feedback has led to improving the product and processing system to address product limitations and improve its applicability for pre-convective forecasting. This presentation highlights the design of the NUCAPS-Forecast system. Solutions were tested and implemented to remove excess moisture and conjoin forward-trajected parcels into vertical profiles for calculating grids of stability indices. Additionally, use cases that explore the tendencies of derived convective available potential energy (CAPE) and convective inhibition (CIN) were investigated. When evaluated against gridded, hourly accumulated rainfall obtained from the Multi-Radar Multi-Sensor (MRMS) observational data set areas with elevated CAPE (reduced CIN) are shown to be associated with areas of precipitation and increase confidence the LEO soundings advected forward in time physically represent atmospheric processes. These results suggest that adiabatic parcel advection of LEO satellite sounding snapshots forward in time exhibits skill in discriminating locations and times that are more likely to lead to convective initiation. Leveraging trajectory modeling to increase the spatial and temporal resolution of LEO IR sounding observations today demonstrates the value of these observations from the future constellation of JPSS satellites and GeoXO mission.

weather↗
NASA NTRSDate not supplied

Atmospheric Formaldehyde Trend and Its Source Attributions in the Recent Decades

Formaldehyde (CH 2 O) is one of the most important reactive trace gases in the atmosphere with an important role in tropospheric chemistry and in the control of surface air quality. Sources of CH 2 O includes direct emissions and in situ production from the oxidation of volatile organic compounds (VOCs). In this study, we investigate spatial-temporal variations in atmospheric CH 2 O and how changes in CH 2 O couple with tropospheric ozone variability during the past two decades using recently available retrievals from the Ozone Monitoring Instrument (OMI) aboard NASA’s Aura satellite and the RefD1 simulation from the NASA Goddard Earth Observing System Chemistry Climate Model (GEOSCCM). Our initial analysis provides evidence of significant increases in atmospheric CH 2 O concentrations over several of the world’s major anthropogenic regions, including China, India, and the Middle East, and several biomass burning regions, including South America, southern Africa, boreal regions and Indonesia. We furthermore will identify regional source contributions to CH 2 O trends, including anthropogenic VOCs emissions, biogenic isoprene emissions and open fires. In particular, we will focus on anthropogenic emission regions to examine the role of changing anthropogenic VOC emissions on CH 2 O abundance.

Atmospheric Formaldehyde↗
NASA NTRSDate not supplied

Summary of the 6th AIAA Propulsion Aerodynamics Workshop Nozzle Test Case: Heated Supersonic Axisymmetric Jets

This paper summarizes findings from the sixth AIAA Propulsion Aerodynamics Workshop (PAW) nozzle test case. The focus of the workshop participants was turbulent supersonic round jets with variation in jet Mach number and jet heating. Three nozzles were utilized with on-design jet Mach numbers of 1.36, 1.63, and 2.00. A heated jet for each of these three nozzles at on-design conditions, a temperature matched on-design condition at Mach 1.63, and an off-design heated condition for the Mach 1.63 nozzle were the five cases examined. Workshop participants utilized Reynolds-averaged Navier-Stokes (RANS) and Large Eddy Simulation (LES) based techniques. Computational fluid dynamics (CFD) solutions were compared with non-intrusive measurements of velocities and temperatures, including both mean values and turbulence statistics. In general, the RANS approaches were unable to capture the trend of reduced mixing as the jet Mach number increased. Incorporation of compressibility corrections improved the agreement with experimental measurements in terms of trends with varying jet Mach number, but also increased the jet potential core lengths to be in worse agreement with data. RANS methods were able to reasonably capture trends of jet heating. The LES-based approaches as a whole were able to capture the effects of compressibility more accurately. Effects of jet heating were captured well by LES. None of the LES approaches introduced disturbances from the wall boundary layers of the nozzle into the jet mixing layer, which would affect the jet flow physics modeling in the initial part of the jet shear layer, but would be very computationally demanding.

supersonic jet↗
NASA NTRSDate not supplied

On the Prediction of Surface Melt Over the Greenland Ice Sheet at NWP and S2S Timescales

Greenland surface melt “events” – encompassing large swaths of the ice sheet over approximately 4-10 day periods – have occurred with increasing frequency in the 21st Century. Surface processes are significant to the ice sheet’s overall contribution to sea level. In addition, the sudden onset of widespread surface melt produces hazardous conditions for researchers and for communities adjacent to the ice sheet. The mechanisms producing melt events, their associated large-scale forcing, and timescales of interest are highly variable and poorly understood. Here we examine the predictability of melt events in a numerical weather prediction system (NWP) and a subseasonal-to-seasonal system (S2S). The NASA Global Modeling and Assimilation Office forward processing (GMAO-FP) is an NWP system that incorporates an evolving, global hybrid 4D ensemble–variational (4DEnVar) data assimilation system and a state-of-the-art model that is run twice-daily with forecasts out to five and ten days. The GMAO S2S system is a hindcast-anomaly, coupled-model prediction system initialized with an atmospheric analysis, and in situ and satellite ocean data assimilation. Both systems represent ice sheet surface hydrology and snowpack processes. We examine the performance of the systems over the period 2018-2023, a span that includes approximately nine significant events, against analyses and satellite-derived melt extent. In addition to extent, we examine related conditions and their role in the melt forecast, including the onset of both omega and Rex atmospheric blocking patterns, sea surface temperatures, sea ice extent, water vapor transport, precipitation, and surface preconditioning. We further examine the tendency for false-alarm events in each system and assess reasons for the forecast performance. North Atlantic atmospheric blocking patterns leading to surface melt are well forecast by the NWP system out to around eight days, while the S2S system generally underestimates the magnitude of such events at monthly lead times. These systems suggest that while short term forecasts of ice sheet melt and melt events are relatively robust, S2S predictions of surface conditions are highly dependent on the ability of a coupled system to generate the short-term atmospheric conditions associated with surface melt events.

Richard Cullather↗
NASA NTRSDate not supplied

Construction of an Exact Pressure-Equilibrium Scheme for the Five-Equation Two-Phase Flow Model With Thermal Relaxation

Numerical simulation of compressible multiphase flows based on the four-equation (homogeneous relaxation) model is known to suffer from two fundamental difficulties with (a) wave propagation and (b) pressure equilibrium preservation. First, the mixture sound speed exhibits non-monotonic dependency with respect to the volume fraction, which leads to robustness issues in the resolution of shocks and acoustic wave propagation across two-phase regions. This difficulty can be mitigated by solving Allaire’s five-equation model augmented with infinitely fast phasic temperature equilibrium, from which solutions of the four-equation model can be recovered. However, when temperature is non-uniform, this augmented five-equation formulation still fails to preserve pressure equilibrium across material interfaces. In this work, we propose a fully conservative numerical scheme that exactly preserves pressure equilibrium at the discrete level for the augmented five-equation model, for arbitrary initial distributions of temperature and volume fraction. Combined with the monotonic sound speed property of the five-equation formulation, the proposed pressure-equilibrium preserving scheme significantly improves robustness in the presence of strong multiphase interactions, including shock–interface interactions and advection of material interfaces.

ESG↗