Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “coupling constant”

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 487 records · Page 27

Room Temperature Thin Film Ba(x)Sr(1-x)TiO3 Ku-Band Coupled MicrostripPhase Shifters: Effects of Film Thickness, Doping, Annealing and Substrate Choice

We report on measurements taken on over twenty Ku-band coupled microstrip phase shifters (CMPS) using thin ferroelectric films of Ba(x)Sr(1-x)TiO3. This CMPS design is a recent innovation designed to take advantage of the high tunability and tolerate the high dielectric constant of ferroelectric films at Ku- and K-band frequencies. These devices are envisioned as a component in low-cost steerable beam phased area antennas, Comparisons are made between devices with differing film thickness, annealed vs unannealed, Mn-doped vs. undoped, and also substrates of LaAlO3 and MgO. A comparison between the CMPS structure and a CPW phase shifter was also made oil the same ferroelectric film.

VanKeuls, F. W.↗

Architecture for Mitigating Short-Term Warning Cosmic Threats: READI Project

Earth is being constantly bombarded by a large variety of celestial bodies and has been since its formation 4.5 billion years ago. Among those bodies, mainly asteroids and comets, there are those that have the potential to create large scale destruction upon impact. The only extinction-level impact recorded to date was 65 million years ago, during the era of dinosaurs. The probability of another extinction-level, or even city-killer, impact may be negligible, but the consequences can be severe for the biosphere and for our species. Therefore it is highly imperative for us to be prepared for such a devastating impact in the near future, especially since humanity is at the threshold of wielding technologies that allow us to do so. Majority of scientists, engineers, and policymakers have focused on long-term strategies and warning periods for Earth orbit crossing Near-Earth Objects (NEOs), and have suggested methods and policies to tackle such problems. However, less attention has been paid to short warning period NEO threats. Such NEOs test current technological and international cooperation capabilities in protecting ourselves, and can create unpredictable devastation ranging from local to global scale. The most recent example is the Chelyabinsk incident in Russia. This event has provided a wakeup call for space agencies and governments around the world towards establishing a Planetary Defense Program. The Roadmap for EArth Defense Initiative (READI) is a project by a team of international, intercultural, and interdisciplinary participants of the International Space University's Space Studies Program 2015 hosted by Ohio University, Athens, OH proposing a roadmap for space agencies, governments, and the general public to tackle NEOs with a short warning before impact. Taking READI as a baseline, this paper presents a technical description of methodologies proposed for detection and impact mitigation of a medium-sized comet (up to 800m across) with a short-warning period of two years on a collision course with Earth. The hypothetical comet is on a highly-inclined orbit having a high probability for Earth impact after its perihelion. For detection, we propose a space-based infrared detection system consisting of two satellites located at the Earth-Moon Lagrange points L1 and L2 coupled with space observatories, like the James Webb telescope and the Centennial telescope. These telescopes are supported by ground-based telescopes, like the Arecibo and Green Bank telescope, in the search for NEOs. Upon detection, the comet is tracked constantly using space- and ground-based telescopes. The deflection system is two-pronged, firstly involving the use of a high energy Directed Energy Laser Terminals (DELT) placed at Sun-Earth Lagrange points L4 and L5 so as to initiate and increase the ablation rate of the comet and deviate it from its collision trajectory, and secondly by the Hypervelocity Comet Intercept Vehicle (HCIV), a space-borne system combining a kinetic impactor with a thermonuclear device. The policy and international collaboration aspects to implement these methods are also outlined in the paper. The techniques mentioned could also be applied to mitigate medium-to-large sized asteroids (up to 2km across).

Laser Applications↗

Planetary Defense Architecture for Mitigating Short-Term Warning Cosmic Threats: READI Project

Earth is being constantly bombarded by a large variety of celestial bodies and has been since its formation 4.5 billion years ago. Among those bodies, mainly asteroids and comets, there are those that have the potential to create large scale destruction upon impact. The only extinction-level impact recorded to date was 65 million years ago, during the era of dinosaurs. The probability of another extinction-level, or even city-killer, impact may be negligible, but the consequences can be severe for the biosphere and for our species. Therefore it is highly imperative for us to be prepared for such a devastating impact in the near future, especially since humanity is at the threshold of wielding technologies that allow us to do so. Majority of scientists, engineers, and policymakers have focused on long-term strategies and warning periods for Earth orbit crossing Near-Earth Objects (NEOs), and have suggested methods and policies to tackle such problems. However, less attention has been paid to short warning period NEO threats. Such NEOs test current technological and international cooperation capabilities in protecting ourselves, and can create unpredictable devastation ranging from local to global scale. The most recent example is the Chelyabinsk incident in Russia. This event has provided a wakeup call for space agencies and governments around the world towards establishing a Planetary Defense Program. The Roadmap for EArth Defense Initiative (READI) is a project by a team of international, intercultural, and interdisciplinary participants of the International Space University's Space Studies Program 2015 hosted by Ohio University, Athens, OH proposing a roadmap for space agencies, governments, and the general public to tackle NEOs with a short warning before impact. Taking READI as a baseline, this paper presents a technical description of methodologies proposed for detection and impact mitigation of a medium-sized comet (up to 800m across) with a short-warning period of two years on a collision course with Earth. The hypothetical comet is on a highly-inclined orbit having a high probability for Earth impact after its perihelion. For detection, we propose a space-based infrared detection system consisting of two satellites located at the Earth-Moon Lagrange points L1 and L2 coupled with space observatories, like the James Webb telescope and the Centennial telescope. These telescopes are supported by ground-based telescopes, like the Arecibo and Green Bank telescope, in the search for NEOs. Upon detection, the comet is tracked constantly using space- and ground-based telescopes. The deflection system is two-pronged, firstly involving the use of a high energy Directed Energy Laser Terminals (DELT) placed at Sun-Earth Lagrange points L4 and L5 so as to initiate and increase the ablation rate of the comet and deviate it from its collision trajectory, and secondly by the Hypervelocity Comet Intercept Vehicle (HCIV), a space-borne system combining a kinetic impactor with a thermonuclear device. The policy and international collaboration aspects to implement these methods are also outlined in the paper. The techniques mentioned could also be applied to mitigate medium-to-large sized asteroids (up to 2km across).

Laser Applications↗

Multiwavelength study of the magnetically active T Tauri star HD 283447

We observed the luminous T Tauri star HD 283447 = V773 Tauri simultaneously at X-ray, ultraviolet, optical photometric and spectroscopic, and radio wavelengths for several hours on UT 1992 September 11. ROSAT, IUE, Very Large Array (VLA) and an intercontinental Very Long Baseline Interferometry (VLBI) network, and three optical observatories participated in the campaign. The star is known for its unusually high and variable nonthermal radio continuum emission. High levels of soft X-ray and Mg II line emission are discovered, with luminosity L(sub x) = 5.5 x 10(exp 30) ergs/s (0.2 - 2 keV) and L(sub Mg II) = 1 x 10(exp 29) ergs/s, respectively. Optically, the spectrum exhibits rather weak characteristics of `classical' T Tauri stars. A faint, broad emission line component, probably due to a collimated wind or infall, is present. During the campaign, the radio luminosity decreased by a factor of 4, while optical/UV lines and X-ray emission remained strong but constant. The large gyrosynchrotron-emitting regions are therefore decoupled from the chromospheric and coronal emission. Five models for the magnetic geometry around the star are discussed; solar-type activity, dipole magnetosphere, star-disk magnetic coupling, disk magnetic fields, and close binary interaction. The data suggest that two magnetic geometries are simultaneously present: complex multipolar fields like those on the Sun, and a large-scale field possibly associated with the circumstellar disk.

Feigelson, Eric D.↗

Combustion Of Moving Droplets And Of Droplets Suspended Within A Convective Environment: Transient Numerical Results

The problem considered is that of a single-component liquid fuel (n-heptane) droplet undergoing evaporation and combustion in a hot, convective, low pressure, zero-gravity environment of infinite expanse. For a moving droplet, the relative velocity (U(sub infinity)) between the droplet and freestream is subject to change due to the influence of the drag force on the droplet. For a suspended droplet, the relative velocity is kept constant. The governing equations for the gas-phase and the liquid-phase consist of the unsteady, axisymmetric equations of mass, momentum, species (gas-phase only) and energy conservation. Interfacial conservation equations are employed to couple the two phases. Variable properties are used in the gas- and liquid-phase. Multicomponent diffusion in the gas-phase is accounted for by solving the Stefan-Maxwell equations for the species diffusion velocities. A one-step overall reaction is used to model the combustion. The governing equations are discretized using the finite volume and SIMPLEC methods. A colocated grid is adopted. Hyperbolic tangent stretching functions are used to concentrate grid points near the fore and aft lines of symmetry and at the droplet surface in both the gas- and liquid-phase. The discretization equations are solved using the ADI method with the TDMA used on each line of the two alternating directions. Iterations are performed within each time-step until convergence is achieved. The grid spacing, size of the computational domain and time-step were tested to ensure that all solutions are independent of these parameters. A detailed discussion of the numerical model is given.

Gogos, George↗

Changes in the Temperature of Saturn's Stratosphere from 2002 to 2004 and Direct Evidence of a Mesopause

The state of a planetary stratosphere is intricately linked to the radiation it receives from its host star. The stellar flux impinging on the the planetary stratosphere affects the temperature and chemistry occurring there. The revolution of the planet about its orbit coupled to the planets axial tilt can cause variations in the amount of flux captured by the planet at a given latitude over the period of a year even if the flux output of the star is constant. This variation of stellar flux at a given latitude during the orbit of a planet is what causes seasons. Saturn, with its axial tilt of 27 degrees, experiences seasonal forcing much like the Earth, albeit with far less solar flux and a year equal to 30 Earth years. In an attempt to measure the effect seasonal forcing has on the temperature structure of Saturn's stratosphere, we have started a long term observing program. By observing Saturn with a high-resolution mid-infrared spectrograph in the v(sub 4) band of methane, we can infer the stratospheric temperature between the pressure levels of 10 and 0.01 mbar. We present latitudinally resolved observations of Saturn's southern hemisphere from 2002 and 2004. We have analyzed the data using a line-by-line radiative transfer model and inferred stratospheric temperatures and winds. We compare the temperatures from the two observing runs and the predictions of a stratospheric seasonal climate model. The zonal winds in the southern hemisphere, between the 10 and 0.01 mbar pressure levels, will be shown and discussed.

Greathouse, Thomas K.↗

Correlating CFD Simulation with Wind Tunnel Test for the Full-Scale UH-60A Airloads Rotor

Data from the recent UH-60A Airloads Test in the National Full-Scale Aerodynamics Complex 40- by 80- Foot Wind Tunnel at NASA Ames Research Center are presented and compared to predictions computed by a loosely coupled Computational Fluid Dynamics (CFD)/Comprehensive analysis. Primary calculations model the rotor in free-air, but initial calculations are presented including a model of the tunnel test section. The conditions studied include a speed sweep at constant lift up to an advance ratio of 0.4 and a thrust sweep at constant speed into deep stall. Predictions show reasonable agreement with measurement for integrated performance indicators such as power and propulsive but occasionally deviate significantly. Detailed analysis of sectional airloads reveals good correlation in overall trends for normal force and pitching moment but pitching moment mean often differs. Chord force is frequently plagued by mean shifts and an overprediction of drag on the advancing side. Locations of significant aerodynamic phenomena are predicted accurately although the magnitude of individual events is often missed.

Romandr, Ethan↗

Highlights from 4STAR Sky-Scanning Retrievals of Aerosol Intensive Optical Properties from Multiple Field Campaigns with Detailed Comparisons of SSA Reported During SEAC4RS

The 4STAR (Spectrometer for Sky-Scanning, Sun-Tracking Atmospheric Research) instrument combines airborne sun tracking capabilities of the Ames Airborne Tracking Sun Photometer (AATS-14) with AERONET (Aerosol Robotic Network)-like sky-scanning capability and adds state-of-the-art fiber-coupled grating spectrometry to yield hyperspectral measurements of direct solar irradiance and angularly resolved sky radiance. The combination of sun-tracking and sky-scanning capability enables retrievals of wavelength-dependent aerosol optical depth (AOD), mode-resolved aerosol size distribution (SD), asphericity, and complex refractive index, and thus also the scattering phase function, asymmetry parameter, single-scattering albedo (SSA), and absorption aerosol optical thickness (AAOT). From 2012 to 2014 4STAR participated in four major field campaigns: the U.S. Dept. of Energy's TCAP (Two-Column Aerosol Project) I & II campaigns, and NASA's SEAC4RS (Studies of Emissions, Atmospheric Composition, Clouds and Climate Coupling by Regional Surveys) and ARISE (Arctic Radiation - IceBridge Sea & Ice Experiment) campaigns. Establishing a strong performance record, 4STAR operated successfully on all flights conducted during each of these campaigns. Sky radiance spectra from scans in either constant azimuth (principal plane) or constant zenith angle (almucantar) were interspersed with direct beam measurements during level legs. During SEAC4RS and ARISE, 4STAR airborne measurements were augmented with flight-level albedo from the collocated Shortwave Spectral Flux Radiometer (SSFR) providing improved specification of below-aircraft radiative conditions for the retrieval. Calibrated radiances and retrieved products will be presented with particular emphasis on detailed comparisons of ambient SSA retrievals and measurements during SEAC4RS from 4STAR, AERONET, HSRL2 (High Spectral Resolution Lidar), and from in situ measurements.

Radiance↗

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO). PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Material Response↗

On the question of adequate hingeless rotor modeling in flight dynamics.

The somewhat controversial question of which elastic blade modes are essential in the flight mechanics of hingeless rotorcraft is studied on the basis of quasi-steady linear aerodynamics including reversed flow effects and uniform inflow. The modes are for the rotating blade, and intermode aerodynamic coupling terms are retained. The criteria for judging elastic mode effects include 19 hub moment and force derivatives, rotor trim data, rotor stability charts for lagged hub moment feedback, step gust and random gust responses. Fixed hub and constant chord blades with widely differing elasticity and inertia and with moderate twist are assumed.

Hohenemser, K. H.↗

Intermode energy transfer in vibrationally excited O3

The laser-excited fluorescence method has been employed to determine the rate constants for vibrational relaxation of the O3 (010), O3 (100) and O3 (001) levels at 298 K. The fluorescence observations from the O3 (010) level provide direct measurements of the rate for intermode vibrational energy transfer from the coupled nu sub 1 and nu sub 3 modes to the nu sub 2 mode. The slowness of this process indicates the likelihood that the nu sub 1 and/or nu sub 3 modes (rather than the nu sub 2 mode) play a predominant role in the laser-enhanced reaction between vibrationally excited O3 and NO at 298 K.

Hui, K.-K.↗

Estimation of the viscosities of silicate liquids at high pressure from measurements of oxygen diffusivity

The dynamics and evolution of a magma ocean depend significantly on the structure and physical properties of silicate liquids at high pressure. Viscosity is particularly important because of its effect on the dynamics of convection and crystal setting. Traditionally, the viscosities of silicate liquids have been measured at high pressure by falling sphere viscometry but, due to some limitations of this technique, few measurements have been made at pressures exceeding 2.5 GPa. An alternative approach, which has previously been proven up to 2 GPa, is to use high pressure measurements of oxygen self-diffusivity (D) to estimate viscosity (eta) from the Eyring relationship eta = kT/D lambda (where k is the Boltzmann constant, T is absolute temperature, and lambda is the diffusive jump distance). We have performed preliminary experiments on Na2Si4O9 liquid in order to explore the feasibility of this method up to pressures in excess of 10 GPa. Diffusion couples were prepared from glass starting materials with one half of the sample enriched in O-18 and the other half containing the natural abundance (0.2 wt percent). Diffusion experiments were performed at 1600 to 1800 C at pressures in the range of 2.5 to 10 GPa for times up to 6 minutes using a 1200 ton multianvil apparatus. Oxygen diffusivities were derived from the resulting O-18 concentration profiles, which were analyzed using an ion microprobe. At 1800 C, oxygen diffusivities increase continuously from 1 x 10(exp -10) m(exp -2) s(exp -1) at 2.5 Gpa to 5 x 10(exp -10) m(exp -2) s(exp -1) at 10 GPa. From the Eyring relationship, the viscosity of Na2SiO9 liquid is predicted to decrease by approx. 0.7 log units over this pressure range at 1800 C. These trends, which can be related to changes in the structure of the liquid at high pressure, are in agreement with results of molecular dynamics calculations. These results suggest that it should be possible to estimate the viscosities of silicate liquids up to at least 15 GPa and 2200 C by this method.

Rubie, D. C.↗

A Summary of the International Workshops on Space Solar Cell Calibration and Measurement Techniques

The first two of a planned series of international workshops concerning space solar cell calibration and measurement techniques have been held within the past year. The need for these workshops arose from the increasing complexity of space solar cells coupled with the growing international nature of the market for space cells and arrays. The workshops, jointly sponsored by NASDA, ESA and NASA, have the objective of obtaining international agreement on standardized values for the AM0 spectrum and constant, recommendations for laboratory measurement practices and the establishment of a set of protocols for making interlaboratory comparison measurements. The results of the first two workshops, held in Waikiki, Hawaii, USA in 1994 and Madrid, Spain in 1995, are presented.

Brinker, D. J.↗

Synthesis and Photoluminescent Properties of Arylethynyl substituted 9,10-Anthraquinones

A series of arylethynyl substituted anthraquinones were synthesized via Sonogashira coupling reactions of 2,7- dibromo-, 2,6-dibromo- and 2-bromoanthraquinone with para-substituted phenylacetylenes. While the redox properties of those compounds are almost insensitive to substitution, their absorption maxima are linearly related to the Hammett constants for electron donating and electron withdrawing groups separately. All compounds are photoluminescent both in solution (quantum yield of emission approximately 2%) and as solids. X-ray crystallographic characterization of 2,7-bisphenylethynyl anthraquinone indicates a monoclinic p2(l/n) space group and no indication for pi-overlap that would promote self-quenching. The emission maxima are red- shifted by both electron donating and electron withdrawing groups alike. The Stokes shifts of all compounds are significant and are correlated to the electronic properties of the substituents. The reduced forms of these compounds are also photoluminescent and the emission originates from the dihydroanthraquinone core.

Yang, Jin-Hua↗

Effect of Operating Frequency on PDE Driven Ejector Thrust Performance

Results of an on-going study of pulse detonation engine driven ejectors are presented and discussed. The experiments were conducted using a pulse detonation engine (PDE) designed to operate at frequencies up to 50 Hz. The PDE used in these experiments utilizes an equi-molar mixture of oxygen and nitrogen as the oxidizer, and ethylene (C2H4) as the fuel, with the propellant mixture having an equivalence ratio of one. A line of sight laser absorption technique was used to determine the time needed for proper filling of the tube. Thrust measurements were made using an integrated spring damper system coupled with a linear variable displacement transducer. The baseline thrust of the PDE was first measured at each desired frequency and agrees with experimental and modeling results found in the literature. Thrust augmentation measurements were then made for constant diameter ejectors. The ejectors had varying lengths, and two different inlet geometries were tested for each ejector configuration. The parameter space for the study included PDE operation frequency, ejector length, overlap distance and the radius of curvature for the ejector inlets. For the studied experimental matrix, the results showed a maximum thrust augmentation of 106% at an operational frequency of 30 Hz.

Santoro, Robert J.↗

Multimodal Pilot Behavior in Multi-Axis Tracking Tasks with Time-Varying Motion Cueing Gains

In a large number of motion-base simulators, adaptive motion filters are utilized to maximize the use of the available motion envelope of the motion system. However, not much is known about how the time-varying characteristics of such adaptive filters affect pilots when performing manual aircraft control. This paper presents the results of a study investigating the effects of time-varying motion filter gains on pilot control behavior and performance. An experiment was performed in a motion-base simulator where participants performed a simultaneous roll and pitch tracking task, while the roll and/or pitch motion filter gains changed over time. Results indicate that performance increases over time with increasing motion gains. This increase is a result of a time-varying adaptation of pilots' equalization dynamics, characterized by increased visual and motion response gains and decreased visual lead time constants. Opposite trends are found for decreasing motion filter gains. Even though the trends in both controlled axes are found to be largely the same, effects are less significant in roll. In addition, results indicate minor cross-coupling effects between pitch and roll, where a cueing variation in one axis affects the behavior adopted in the other axis.

motion simulators↗

Modeled Impacts of LULC and Climate Change Predictions on Water Quality in Belize

Land use and land cover (LULC) change can have detrimental effects on water quality. In Belize, agricultural expansion creates the risk of increased sediment load and excess nutrients in runoff water, while deforestation removes potential infiltration sites for this outflow. Climate change and evolving precipitation rates can intensify the quantity of runoff, further enabling the flow of sediments and excess nutrients out to the lagoon surrounding the Belize Barrier Reef Reserve System (BBRRS). This study estimated future water quality in Belize by projecting LULC change through 2090 across Belize’s major watersheds based on observed trends from 2008-2018. Those LULC projections were subsequently combined with soil type data, elevation, and precipitation rates into a hydrologic model to produce runoff estimates as a proxy for water quality. The two Coupled Model Intercomparison Project phase 6 (CMIP-6) scenarios employed in the study represented bookend climate change scenarios, and both predicted lower precipitation rates in Belize over the next century due to climate change. When holding LULC change constant, these climate scenarios projected a decrease in runoff, suggesting a positive relationship between precipitation and runoff. The southern watersheds in particular are projected to experience greater decreases in annual rainfall and runoff by 2090. When holding climate constant, runoff increased by approximately 2.8% in the Conservation-focused LULC scenario by 2090, which was 28% lower than the Business as Usual scenario and 42% lower than the Development scenario.

water resources↗

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO) [1,2,3]. PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) [4] CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) [5] program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield [6]. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem [7]. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Thermal Protection Systems↗