Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surface complexation modeling”

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 109 records · Page 6

Conservative force model performance for TOPEX/Poseidon precision orbit determination

The TOPEX/Poseidon spacecraft was launched on August 10, 1992 to study the Earth's oceans. To achieve maximum benefit from the altimetric data collected, mission requirements dictate that TOPEX/Poseidon's orbit must be computed at an unprecedented level of accuracy. In order to satisfy these requirements, a model which accounts for the satellite's complex geometry, attitude, and surface properties has been developed. This `box-wing' representation treats the spacecraft as the combination of flat plates arranged in the shape of a box and a connecetd solar array. The nonconservative forces acting on each of the eight surfaces are computed independently, yielding vector accelerations which are summed to compute the total aggregate effect on the satellite center-of-mass. Parameters associated with each flat plate were derived from a finite element analysis of the spacecraft. Certain parameters can be inferred from tracking data and have been adjusted to obtain a better representation of the satellite acceleration history. Changes in the nominal mission profile and the presence of an `anomalistic' force have complicated this tuning process. Model performance, parameter sensitivities, and the `anomalistic' force will be discussed.

Marshall, J. Andrew↗

Evaluation of NASA's High-Resolution Global Composition Simulations: Understanding a Pollution Event in the Chesapeake Bay During the Summer 2017 OWLETS Campaign

Recirculation of pollutants due to a bay breeze effect is a key meteorological mechanism impacting air quality near urban coastal areas, but regional and global chemical transport models have historically struggled to capture this phenomenon. We present a case study of a high ozone (O3) episode observed over the Chesapeake Bay during the NASA Ozone Water-Land Environmental Transition Study (OWLETS) in summer 2017. OWLETS included a complementary suite of ground-based and airborne observations, with which we characterize the meteorological and chemical context of this event and develop a framework to evaluate model performance. Two publicly-available NASA global high-resolution coupled chemistry-meteorology models (CCMMs) are investigated: GEOS-CF and MERRA2-GMI. The GEOS-CF R squared value for comparisons between the NASA Sherpa C-23 aircraft measurements to the GEOS-CF resulted in good agreement (R squared: 0.67) on July 19th and fair agreement (R squared: 0.55) for July 20th. Compared to surface observations, we find the GEOS-CF product with a 25 x 25 km squared grid box, at an hourly (R squared: 0.62 to 0.87) and 15-minute (R squared: 0.64 to 0.87) interval for six regional sites outperforms the hourly nominally 50 x 50 km squared gridded MERRA2-GMI (R squared: 0.53 to 0.76) for four of the six sites, suggesting it is better capable of simulating complex chemical and meteorological features associated with ozone transport within the Chesapeake Bay airshed. When the GEOS-CF product was compared to the TOLNet LiDAR observations at both NASA Langley Research Center (LaRC) and the Chesapeake Bay Bridge Tunnel (CBBT), the median differences at LaRC were -6 to 8% and at CBBT were ± 7% between 400 to 2000 m ASL. This indicates that, for this case study, the GEOS-CF is able to simulate surface level ozone diurnal cycles and vertical ozone profiles at small scales between the surface level and 2000 m ASL. Evaluating global chemical model simulations at sub-regional scales will help air quality scientists understand the complex processes occurring at small spatial and temporal scales within complex surface terrain changes, simulating nighttime chemistry and deposition, and the potential to use global chemical transport simulations in support of regional and sub-regional field campaigns.

NASA Ozone Water-Land Environmental Transition Stu↗

Hector V3.2.0: Functionality and Performance of a Reduced-Complexity Climate Model

Hector is an open-source reduced complexity climate-carbon cycle model that models critical Earth system processes on a global and annual basis. Here we present an updated version of the model, Hector V3.2.0 (hereafter Hector V3) and document its new features, implementation of new science, and performance. Significant new features include permafrost thaw, a reworked energy balance submodel, and updated parameterizations throughout. Hector V3 results are in good general agreement with historical observations of atmospheric CO2 concentrations and global mean surface temperature, and its future temperature projections are consistent with more complex Earth System Model output data from the Sixth Coupled Model Intercomparison Project. We show that Hector V3 is a fully open source, flexible, performant, and robust simulator of global climate changes, note its limitations, and discuss future areas of improvement and research with respect to the model’s scientific, stakeholder, and educational priorities.

Kayln Dorheim↗

Diagnosing the Nature of Land-Atmosphere Coupling During the 2006-7 Dry/Wet Extremes in the U. S. Southern Great Plains

The degree of coupling between the land surface and PBL in NWP models remains largely undiagnosed due to the complex interactions and feedbacks present across a range of scales. In this study, a framework for diagnosing local land-atmosphere coupling (LoCo) is presented using a coupled mesoscale model with observations during the summers of 2006/7 in the U.S. Southern Great Plains. Specifically, the Weather Research and Forecasting (WRF) model has been coupled to NASA's Land Information System (LIS), which enables a suite of PBL and land surface model (LSM) options along provides a flexible and high-resolution representation and initialization of land surface physics and states. This coupling is one component of a larger project to develop a NASA-Unified WRF (NU-WRF) system. A range of diagnostics exploring the feedbacks between soil moisture and precipitation are examined for the dry/wet extremes, along with the sensitivity of PBL-LSM coupling to perturbations in soil moisture.

Santanello, Joseph A., Jr.↗

Regional Impacts of Ice Sheet Surface Representation Improvements in an Earth System Model

Ice sheet surface conditions and their subsequent impact on surface mass balance play an important role in ice sheet dynamics and ice sheet interaction with the overlying atmosphere and surrounding ocean. The NASA Global Modeling and Assimilation Office’s (GMAO) Goddard Earth Observing System (GEOS) ecosystem of models and reanalyses model ice sheet and glacier surface mass balance with a moderately complex snow and ice module, including prognostic surface albedo evolution, fractional snow cover, and snowpack hydrology and meltwater retention. With this configuration, Greenland near surface air temperatures and ice sheet melt are well represented in MERRA-2, GMAO’s current atmospheric reanalysis; however, the similarity to observations in some regions is due, in part, to compensating biases in surface energy budget. Recent improvements to GEOS’ cloud microphysics and longwave radiation parameterizations reduced biases in surface net longwave radiation and exacerbated surface net shortwave radiation biases, effectively increasing near surface temperatures and melting, particularly during the summer months. Here, we demonstrate the Earth system response to an improved representation of the ice sheet surface focused on minimizing the surface energy biases and increasing spatial and temporal representativeness. New changes to the surface albedo parameterization, including the introduction of a prognostic aerosol darkening scheme, help mitigate shortwave biases, while modifications to surface roughness characteristics improve near surface sensible heat fluxes. Focusing on the Arctic and Greenland, we examine the combined impact on these changes on regional weather conditions and ice sheet surface mass balance using a high resolution atmospheric only experiment, and we use a coupled model configuration to assess the impact on oceanic conditions and subseasonal-to-seasonal prediction. Overall, these improvements increase modeled ice sheet surface realism and provide a strong basis for more accurate ice sheet mass balance in GMAO’s future reanalyses and forecasting systems.

Lauren C Andrews↗

Regional Impacts of Ice Sheet Surface Representation Improvements in an Earth System Model

Ice sheet surface conditions and their subsequent impact on surface mass balance play an important role in ice sheet dynamics and ice sheet interaction with the overlying atmosphere and surrounding ocean. The NASA Global Modeling and Assimilation Office’s (GMAO) Goddard Earth Observing System (GEOS) ecosystem of models and reanalyses model ice sheet and glacier surface mass balance with a moderately complex snow and ice module, including prognostic surface albedo evolution, fractional snow cover, and snowpack hydrology and meltwater retention. With this configuration, Greenland near surface air temperatures and ice sheet melt are well represented in MERRA-2, GMAO’s current atmospheric reanalysis; however, the similarity to observations in some regions is due, in part, to compensating biases in surface energy budget. Recent improvements to GEOS’ cloud microphysics and longwave radiation parameterizations reduced biases in surface net longwave radiation and exacerbated surface net shortwave radiation biases, effectively increasing near surface temperatures and melting, particularly during the summer months. Here, we demonstrate the Earth system response to an improved representation of the ice sheet surface focused on minimizing the surface energy biases and increasing spatial and temporal representativeness. New changes to the surface albedo parameterization, including the introduction of a prognostic aerosol darkening scheme, help mitigate shortwave biases, while modifications to surface roughness characteristics improve near surface sensible heat fluxes. Focusing on the Arctic and Greenland, we examine the combined impact on these changes on regional weather conditions and ice sheet surface mass balance using a high resolution atmospheric only experiment, and we use a coupled model configuration to assess the impact on oceanic conditions and subseasonal-to-seasonal prediction. Overall, these improvements increase modeled ice sheet surface realism and provide a strong basis for more accurate ice sheet mass balance in GMAO’s future reanalyses and forecasting systems.

Lauren C. Andrews↗

Using a fast Fourier method to model sound propagation in a stratified atmosphere over a stratified porous-elastic ground

Using a Fast Fourier integration method and a global matrix method for solution of the boundary condition equations at all interfaces simultaneously, a useful tool for predicting acoustic propagation in a stratified fluid over a stratified porous-elastic solid was developed. The model for the solid is a modified Biot-Stoll model incorporating four parameters describing the pore structure corresponding to the Rayleigh-Attenborough rigid-porous structure model. The method is also compared to another Fast Fourier code (CERL-FFP) which models the ground as an impedance surface under a horizontally stratified air. Agreement with the CERL FFP is good. The effects on sound propagation of a combination of ground elasticity, complex ground structure, and atmospheric conditions are demonstrated by theoretical results over a snow layer, and experimental results over a model ground surface.

Tooms, S.↗

A Monte Carlo/response surface strategy for sensitivity analysis: application to a dynamic model of vegetative plant growth

We describe the application of a strategy for conducting a sensitivity analysis for a complex dynamic model. The procedure involves preliminary screening of parameter sensitivities by numerical estimation of linear sensitivity coefficients, followed by generation of a response surface based on Monte Carlo simulation. Application is to a physiological model of the vegetative growth of soybean plants. The analysis provides insights as to the relative importance of certain physiological processes in controlling plant growth. Advantages and disadvantages of the strategy are discussed.

Non-NASA Center↗

Identification and Study of Validation Level Test Cases for Computational Modeling of Non-Charring Ablators

Computational modeling of Thermal Protection System (TPS) materials, used for aerospace applications, provides numerous advantages in preliminary selection and design of a heatshield material and shape for atmospheric entry vehicles. However, to serve as a reliable tool for prediction of material thermal and ablative behavior, the modeling approach needs to be validated against real experimental and flight data, preferably at a range of applied conditions. The validation study is typically very complex as it requires reliable measured data not only for the material thermal response and surface recession, but also well characterized environmental conditions. The validation problem becomes even more complex when the material thermal response is dictated by multi-physics effects such as solid conduction, in-depth thermal decomposition, pyrolysis gas flow and chemical reactions. The multi-physics effects complicate not only the modeling effort, but also the experimental measurement for validation of various aspects of the highly coupled problem. In this study, an attempt is made to identify suitable experimental data that could serve as a source for validation of material thermal response modeling tools. To reduce the computational complexity, this study focuses only on non-charring ablators, where the material thermal response could be modeled with a single governing equation for solid conduction and the ablation is limited only to the surface of the material. With a well characterized and publicly available experimental data being sparse, the study is limited in presenting test cases for only three materials: camphor, graphite and FiberForm® in the sequence of increased modeling complexity. Graphite is a commonly used TPS material for aerospace applications, both for leading edges of high-speed vehicles and internal insulation of solid rocket motors. FiberForm® is a porous carbon pre-form used in preparation of the well known PICA material Tran et al. [1996]. Inclusion of camphor into the list is conditioned with the relative simplicity in modeling the material thermal and chemical response and the low-enthalpy flow environment. In addition, camphor has been used as a simple test material for study of flow transition behavior by Stock and Ginoux [1973] and assessment of a heatshield shape change at flight relevant conditions by Rotondi et al. [2022]. In this work, the identified experimental data was extracted from the public literature and test cases that yet have been published. As it was found from the review, not a single test case contains an exhaustive set of data that would validate every aspect of the material physics. However, in the data collected, various aspects of the material behavior can be still validated, such as surface and in-depth temperature, amount of recession and a shape change. The identified experimental data for each case is accompanied with a characterized flow environment and simulated boundary conditions predicted by a Data-Parallel Line Relaxation (DPLR) code Wright et al. [1998]. In addition, material thermal response numerical simulations in each test case are performed with Kentucky Aerothermodynamics and Thermal Response System (KATS-MR) Zibitsker et al. [2022] providing a comparative study and a sanity check for the proposed validation data. Sample results from the performed numerical study are shown below. Figure 1 shows distribution of surface heat flux and pressure values on a hemi-cylinder model made of FiberForm® and tested in HyMETS arc-jet facility. The results are shown for the high pressure condition among the two tests. Flow simulation was performed with DPLR code on a quarter of original geometry. In the figure, the quarter shape was mirrored across zx and xy planes to show the complete distribution. Figure 2 shows the material response results for the high pressure case (7500 Pa), simulated with KATS-MR and a comparison to the experimental data for the surface temperature and shape shape. The simulation was performed on a 2-D slice, extracted in the xy plane at the middle of the sample. Figure 3 shows the material response simulation for the low pressure case (3500 Pa) and a comparison to the experimental data for surface temperature and shape change.

ablation↗

Unforced decadal fluctuations in a coupled model of the atmosphere and ocean mixed layer

Global average temperature in a 100-year control run of a model used for greenhouse gas response simulations showed low-frequency natural variability comparable in magnitude to that observed over the last 100 years. The model variability was found to be barotropic in the atmosphere, and located in the tropical strip with largest values near the equator in the Pacific. The model variations were traced to complex, low-frequency interactions between the meridional sea surface temperature gradients in the eastern equatorial Pacific, clouds at both high and low levels, and features of the tropical atmospheric circulation. The variations in these and other model parameters appear to oscillate between two limiting climate states. The physical scenario accounting for the oscillations on decadal time scales is almost certainly not found in the real world on shorter time scales due to limited resolution and the omission of key physics (e.g., equatorial ocean dynamics) in the model. The real message is that models with dynamical limitations can still produce significant long-term variability. Only a thorough physical diagnosis of such simulations and comparisons with decadal-length data sets will allow one to decide if faith in the model results is, or is not, warranted.

Barnett, T. P.↗

Ground-Level Digital Terrain Model (DTM) Construction from Tandem-X InSAR Data and Worldview Stereo-Photogrammetric Images

The ground-level digital elevation model (DEM) or digital terrain model (DTM) information are invaluable for environmental modeling, such as water dynamics in forests, canopy height, forest biomass, carbon estimation, etc. We propose to extract the DTM over forested areas from the combination of interferometric complex coherence from single-pass TanDEM-X (TDX) data at HH polarization and Digital Surface Model (DSM) derived from high-resolution WorldView (WV) image pair by means of random volume over ground (RVoG) model. The RVoG model is a widely and successfully used model for polarimetric SAR interferometry (Pol-InSAR) technique for vertical forest structure parameter retrieval [1][2][3][4]. The ground-level DEM have been obtained by complex volume decorrelation in the RVoG model with the DSM using stereo-photogrammetric technique. Finally, the airborne lidar data were used to validate the ground-level DEM and forest canopy height results.

Lee, Seung-Kuk↗

Using SpaceClaim/TD Direct for Modeling Components with Complex Geometries for the Thermal Desktop-Based Advanced Stirling Radioisotope Generator Model

SpaceClaim/TD Direct has been used extensively in the development of the Advanced Stirling Radioisotope Generator (ASRG) thermal model. This paper outlines the workflow for that aspect of the task and includes proposed best practices and lessons learned. The ASRG thermal model was developed to predict component temperatures and power output and to provide insight into the prime contractors thermal modeling efforts. The insulation blocks, heat collectors, and cold side adapter flanges (CSAFs) were modeled with this approach. The model was constructed using mostly TD finite difference (FD) surfaces solids. However, some complex geometry could not be reproduced with TD primitives while maintaining the desired degree of geometric fidelity. Using SpaceClaim permitted the import of original CAD files and enabled the defeaturing repair of those geometries. TD Direct (a SpaceClaim add-on from CRTech) adds features that allowed the mark-up of that geometry. These so-called mark-ups control how finite element (FE) meshes were generated and allowed the tagging of features (e.g. edges, solids, surfaces). These tags represent parameters that include: submodels, material properties, material orienters, optical properties, and radiation analysis groups. TD aliases were used for most tags to allow analysis to be performed with a variety of parameter values. Domain-tags were also attached to individual and groups of surfaces and solids to allow them to be used later within TD to populate objects like, for example, heaters and contactors. These tools allow the user to make changes to the geometry in SpaceClaim and then easily synchronize the mesh in TD without having to redefine these objects each time as one would if using TD Mesher.The use of SpaceClaim/TD Direct has helped simplify the process for importing existing geometries and in the creation of high fidelity FE meshes to represent complex parts. It has also saved time and effort in the subsequent analysis.

Power Systems↗

Mean Flow Augmented Acoustics in Rocket Systems

Oscillatory motion in solid rocket motors and liquid engines has long been a subject of concern. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process and gas dynamics. The customary approach to modeling acoustic waves inside a rocket chamber is to apply the classical inhomogeneous wave equation to the combustion gas. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while the acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The converging section of a rocket nozzle, where gradients in pressure, density, and velocity become large, is a notable region where this approach is not applicable. The expulsion of unsteady energy through the nozzle of a rocket is identified as the predominate source of acoustic damping for most rocket systems. An accurate model of the acoustic behavior within this region where acoustic modes are influenced by the presence of a steady mean flow is required for reliable stability predictions. Recently, an approach to address nozzle damping with mean flow effects was implemented by French [1]. This new approach extends the work originated by Sigman and Zinn [2] by solving the acoustic velocity potential equation (AVPE) formulated by perturbing the Euler equations [3]. The acoustic velocity potential (psi) describing the acoustic wave motion in the presence of an inhomogeneous steady high-speed flow is defined by, (del squared)(psi) − (lambda/c)(exp 2)(psi) − M(dot)[M(dot)(del)(del(psi))] − 2(lambda(M/c) + (M(dot)del(M))(dot)del(psi)−2(lambda)(psi)[M(dot)del(1/c)]=0 (1) with M as the Mach vector, c as the speed of sound, and lambda as the complex eigenvalue. French apply the finite volume method to solve the steady flow field within the combustion chamber and nozzle with inviscid walls. The complex eigenvalues and eigenvector are determined with the use of the ARPACK eigensolver. The present study employs the COMSOL Multphysics framework to solve the coupled eigenvalue problem using the finite element approach. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluated the gas flow inside of a solid rocket motor using St. Robert's law modeling solid propellant burn rate, slip boundary conditions, and the supersonic outflow condition. Results from the HMNF model are used by the coefficient form of the mathematics module to determine the eigenvalues of the AVPE. The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. Pertinent results from these analyses are the complex valued eigenvalue and eigenvectors. Comparisons are made to the French results to evaluate the modeling approach. A comparison of the French results with that of the present analysis is displayed in figures 1 and 2, respectively. The graphic shows the first tangential eigenvector's real (a) and imaginary (b) values.

Fischbach, Sean R.↗

Grand challenge scientific questions in coupled modeling

Most convective field experiments in the past (e.g., SESAME, CCOPE, CINDE) have attempted to resolve only the immediate scales of moist convection using network arrays that spanned two or three atmospheric scales at most. Furthermore, these scales have been defined more on practical considerations (cost, manpower, etc.) than on a clear understanding of their theoretical significance. Unfortunately, this has precluded a description of the entire life cycle of MCS's and their interaction with larger scale systems, the land surface, and trace species. Fortunately, the following factors now make it possible to attempt to simulate scale contraction processes from the synoptic scale down to the cloud scale, as well as interactions between complex meteorological, land surface, precipitation, chemical, and hydrologic processes with coupled, multiscale models: the availability of new technology to sample meteorological fields at high temporal and spatial resolution over a broad region made possible by the weather observing modernization program; increased computer power and improved numerical approaches to run limited area models with nonhydrostatic precipitation physics so as to explicitly resolve MCS (Mesoscale Convective System) processes; and four dimensional assimilation of non-conventional data to provide dynamically consistent datasets for diagnostic analysis of nonlinear scale-interactive dynamics. Several examples of scale-interactive processes which present grand challenges for coupled, multiscale modeling were presented.

Koch, Steven↗

Comparison of NASCAP modelling results with lumped circuit analysis

Engineering design tools that can be used to predict the development of absolute and differential potentials by realistic spacecraft under geomagnetic substorm conditions are described. Two types of analyses are in use: (1) the NASCAP code, which computes quasistatic charging of geometrically complex objects with multiple surface materials in three dimensions; (2) lumped element equivalent circuit models that are used for analyses of particular spacecraft. The equivalent circuit models require very little computation time, however, they cannot account for effects, such as the formation of potential barriers, that are inherently multidimensional. Steady state potentials of structure and insulation are compared with those resulting from the equivalent circuit model.

Stang, D. B.↗

Study of Graphite Ablation at Arc-Jet Conditions using Finite-Rate and Equilibrium Chemistry Models

Moving at hypersonic speeds through the atmosphere generates immense amount of heat that triggers thermo-chemical non-equilibrium effects in the flow approaching the vehicle and leads to surface chemical reactions that cause ablation and consume the material mass. This complex behavior and interaction of the aerothermal environment with a graphite surface is com- putationally modeled in the present work using a fully-coupled simulation between flow solver CHAMPS NBS-Cart and material response solver KATS-MR. In this study, we account for the thermo-chemical non-equilibrium effects occurring in the flow and use an 11-species air-carbon re- action model to simulate the ablation of the material. The reactions at the surface such oxidation, nitridation and recombination are modeled with a newly developed air-carbon reaction model. The coupled simulation is applied to study two experimental cases performed at IHF arc-jet facility at NASA Ames research center. The obtained results are compared against the measured data of surface recession, as well as surface and in-depth temperatures. Additionally, uncoupled sim- ulations, performed with equilibrium chemistry at the wall, point out the substantial differences between the coupled and uncoupled solutions, leading to higher deviation of the predicted temper- atures. Overall, we show excellent agreement between the coupled simulation predictions and the experimental data in all measured quantities.

EDL↗

Simulation gravity modeling to spacecraft-tracking data - Analysis and application

It is proposed that line-of-sight gravity measurements derived from spacecraft-tracking data can be used for quantitative subsurface density modeling by suitable orbit simulation procedures. Such an approach avoids complex dynamic reductions and is analogous to the modeling of conventional surface gravity data. This procedure utilizes the vector calculations of a given gravity model in a simplified trajectory integration program that simulates the line-of-sight gravity. Solutions from an orbit simulation inversion and a dynamic inversion on Doppler observables compare well (within 1% in mass and size), and the error sources in the simulation approximation are shown to be quite small. An application of this technique is made to lunar crater gravity anomalies by simulating the complete Bouguer correction to several large young lunar craters. It is shown that the craters all have negative Bouguer anomalies.

Phillips, R. J.↗

Automated thermal mapping techniques using chromatic image analysis

Thermal imaging techniques are introduced using a chromatic image analysis system and temperature sensitive coatings. These techniques are used for thermal mapping and surface heat transfer measurements on aerothermodynamic test models in hypersonic wind tunnels. Measurements are made on complex vehicle configurations in a timely manner and at minimal expense. The image analysis system uses separate wavelength filtered images to analyze surface spectral intensity data. The system was initially developed for quantitative surface temperature mapping using two-color thermographic phosphors but was found useful in interpreting phase change paint and liquid crystal data as well.

Buck, Gregory M.↗