Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Response surface model”

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 181 records · Page 10

Thorium concentrations in the lunar surface. II - Deconvolution modeling and its application to the regions of Aristarchus and Mare Smythii

The broad angular response which characterized the Apollo gamma ray spectrometer resulted in a loss of spatial resolution and some of the contrast in determining surface concentrations within lunar regions small compared to the field of view. A deconvolution technique has been developed which removes much of this instrumental effect, thereby improving both spatial resolution and accuracy at the cost of a loss in precision. Geometric models of regional thorium distribution are convoluted through the response function of the instrument to yield a predicted distribution which is compared with the observed data field for quality of fit. Application to areas which include Aristarchus and Mare Smythii confirm some geological relationships and fail to support others.

Haines, E. L.↗

Assessment of Pioneer Venus Entry Heating with Coupled Radiation and Ablation

The objective of this work is to reconstruct thermocouple temperature measurements taken at the stagnation point of the Pioneer Venus large probe using uncoupled and coupled ablation calculations. The time history of the large probe forebody entry heating was obtained from a trajectory-based analysis that employed a material response solver to model equilibrium and nonequilibrium ablation chemistry at the surface. In the uncoupled approach, equilibrium ablation chemistry was assumed in the material response calculations, but the material response outputs were not returned to the flowfield solver. Therefore, the effects of char and pyrolysis gas injection into the boundary layer due to ablation (e.g., convective blockage and radiative energy absorption) did not influence the flowfield solutions. In the coupled approach, nonequilibrium ablation chemistry was performed and the material response outputs were returned to the flowfield solver. Therefore, the effects of char and pyrolysis gas injection into the boundary layer were accounted for in the flowfield solutions. Both uncoupled and coupled approaches produced reconstructed temperatures that were within 10% of the flight measurements through the entry heat pulse, with the coupled results comparing more favorably with the measurements. Overall, the reconstructed thermocouple temperatures in this work significantly improved on the results of previous reconstructions for the large probe.

Jannuel V V Cabrera↗

Decoupling of Iron and Phosphate in the Global Ocean

Iron is an essential micronutrient for marine phytoplankton, limiting their growth in high nutrient, low chlorophyll regions of the ocean. I use a hierarchy of ocean circulation and biogeochemistry models to understand controls on global iron distribution. I formulate a mechanistic model of iron cycling which includes scavenging onto sinking particles and complexation with an organic ligand. The iron cycle is coupled to a phosphorus cycling model. Iron's aeolian source is prescribed. In the context of a highly idealized multi-box model scheme, the model can be brought into consistency with the relatively sparse ocean observations of iron in the oceans. This biogeochemical scheme is also implemented in a coarse resolution ocean general circulation model. This model also successfully reproduces the broad regional patterns of iron and phosphorus. In particular, the high macronutrient concentrations of the Southern Ocean result from iron limitation in the model. Due to the potential ability of iron to change the efficiency of the carbon pump in the remote Southern Ocean, I study Southern Ocean surface phosphate response to increased aeolian dust flux. My box model and GCM results suggest that a global ten fold increase in dust flux can support a phosphate drawdown of 0.25-0.5 micromolar.

GLOBAL OCEAN MODEL↗

Mesoscale Convective Systems During SCSMEX: Simulations with a Regional Climate Model and a Cloud-Resolving Model

The South China Sea Monsoon Experiment (SCSMEX) was conducted in May-June 1998. One of its major objectives is to better understand the key physical processes for the onset and evolution of the summer monsoon over Southeast Asia and southern China. Multiple observation platforms (e.g., upper-air soundings, Doppler radar, ships, wind profilers, radiometers, etc.) during SCSMEX provided a first attempt at investigating the detailed characteristics of convection and circulation changes associated with monsoons over the South China Sea region. SCSMEX also provided precipitation derived from atmospheric budgets and comparison to those obtained from the Tropical Rainfall Measuring Mission (TRMM). In this paper, a regional scale model (with grid size of 20 km) and Goddard Cumulus Ensemble (GCE) model (with 1 km grid size) are used to perform multi-day integration to understand the precipitation processes associated with the summer monsoon over Southeast Asia and southern China. The regional climate model is used to understand the soil-precipitation interaction and feedback associated with a flood event that occurred in and around China's Yantz River during SCSMEX Sensitivity tests on various land surface models, sea surface temperature (SST) variations, and cloud processes are performed to understand the precipitation processes associated with the onset of the monsoon over the S. China Sea during SCSMEX. These tests have indicated that the land surface model has a major impact on the circulation over the S. China Sea. Cloud processes can effect the precipitation pattern while SST variation can effect the precipitation amounts over both land and ocean. The exact location (region) of the flooding can be effected by the soil-rainfall feedback. The GCE-model results captured many observed precipitation characteristics because it used a fine grid size. For example, the model simulated rainfall temporal variation compared quite well to the sounding-estimated rainfall. The results show there are more latent heat fluxes prior to the onset of the monsoon. However, more rainfall was simulated after the onset of the monsoon. This modeling study indicates the latent heat fluxes (or evaporation) have more of an impact on precipitation processes and rainfall in the regional climate model simulations than in the cloud-resolving model simulations. Research is underway to determine if the difference in the grid sizes or the moist processes used in these two models is responsible for the differing influence of surface fluxes an precipitation processes.

Tao, W.-K.↗

Inferred Sea Level Prediction in the NASA GMAO Seasonal Forecasting System

Reliable predictions of sea level anomalies on seasonal timescales with lead times of 1 to 9 months may have relevance to stakeholders – for example, in the advance deployment of resources for coastal flood mitigation. Routine prediction and analysis may also highlight physical processes associated with sea level change and modeling capabilities on seasonal and other timescales. These forecasts may represent interannual changes in the seasonal slope of the ocean surface, teleconnection effects such as the El Niño/Southern Oscillation phenomenon, and variations in seasonal hydrology including precipitation and coastal runoff. Coupled atmosphere/ocean models are routinely used in the seasonal prediction of temperature anomalies, precipitation anomalies, sea ice cover, and climate indices such as the Niño3.4 predictions under the North American Multi-Model Ensemble (NMME) protocol. Within the limits of their configuration, these complex Earth-system models have a potential for depicting regional changes in oceanic column properties, including the sea surface height. Seasonal prediction models generally have no representation of long-term mass contributions from melting land ice, or changes in vertical land motion; their output may be more specifically characterized as predictions of the ocean dynamic sea level. In practice however, the sea surface height prognostic variable is substantially compromised by the forecast model response to initial conditions. Imbalances between the initial, observed hydrologic cycle and the forecast model state produce abrupt adjustments in the model sea surface height. As a result, most seasonal prediction systems employ a constraint on the globally-averaged sea surface height that is applied at each time step. This essentially renders the prognostic sea surface height variable as unserviceable. Several approaches have previously been used to retrieve sea level information from seasonal forecasts beyond the use of the sea surface height variable. Here, we extend a method of relating other prognostic values, including ocean circulation and climate indices, to observed sea level variations. We use the merged altimetry record of the NASA MEaSUREs Gridded Sea Surface Height Anomalies data set and monthly revised local reference gauge observations from the National Oceanography Centre Permanent Service for Mean Sea Level (PSMSL) to evaluate derived prognostic variables from the NASA Global Modeling and Assimilation Office subseasonal-to-seasonal system version 2.1 (GMAO S2S v2.1). We focus on results for the midlatitudes with particular emphasis on US gauge locations. As shown in previous studies, prognostic ENSO-related indices in boreal winter are well correlated with gauge observations for the US west coast, but also for other locations in the southeastern US. Other forecast climate indices such as the North Atlantic Oscillation have relations to sea level that are limited both seasonally and spatially. As expected, surface atmospheric pressure (e.g., inverse barometer effect) is found to be particularly well correlated with observed sea level. We provide a characterization of forecast skill for seasonal sea level with this method.

Richard I Cullather↗

Heat Transfer Analysis in Wire Bundles for Aerospace Vehicles

Design of wiring for aerospace vehicles relies on an understanding of "ampacity" which refers to the current carrying capacity of wires, either, individually or in wire bundles. Designers rely on standards to derate allowable current flow to prevent exceedance of wire temperature limits due to resistive heat dissipation within the wires or wire bundles. These standards often add considerable margin and are based on empirical data. Commercial providers are taking an aggressive approach to wire sizing which challenges the conventional wisdom of the established standards. Thermal modelling of wire bundles may offer significant mass reduction in a system if the technique can be generalized to produce reliable temperature predictions for arbitrary bundle configurations. Thermal analysis has been applied to the problem of wire bundles wherein any or all of the wires within the bundle may carry current. Wire bundles present analytical challenges because the heat transfer path from conductors internal to the bundle is tortuous, relying on internal radiation and thermal interface conductance to move the heat from within the bundle to the external jacket where it can be carried away by convective and radiative heat transfer. The problem is further complicated by the dependence of wire electrical resistivity on temperature. Reduced heat transfer out of the bundle leads to higher conductor temperatures and, hence, increased resistive heat dissipation. Development of a generalized wire bundle thermal model is presented and compared with test data. The steady state heat balance for a single wire is derived and extended to the bundle configuration. The generalized model includes the effects of temperature varying resistance, internal radiation and thermal interface conductance, external radiation and temperature varying convective relief from the free surface. The sensitivity of the response to uncertainties in key model parameters is explored using Monte Carlo analysis.

S L Rickman↗

Satellite Soil Moisture Data Assimilation Impacts on Modeling Weather Variables and Ozone in the Southeastern US – Part 2: Sensitivity to Dry-Deposition Parameterizations

Ozone (O3) dry deposition is a major O3 sink. As a follow-up study of Huang et al. (2021), we quantify the impact of satellite soil moisture (SM) on model representations of this process when different dry-deposition parameterizations are implemented, based on which the implications for interpreting O3 air pollution levels and assessing the O3 impacts on human and ecosystem health are provided. The SM data from NASA's Soil Moisture Active Passive mission are assimilated into the Noah-Multiparameterization (Noah-MP) land surface model within the NASA Land Information System framework, semicoupled with Weather Research and Forecasting model with online Chemistry (WRF-Chem) regional-scale simulations covering the southeastern US. Major changes in the modeling system used include enabling the dynamic vegetation option, adding the irrigation process, and updating the scheme for the surface exchange coefficient. Two dry-deposition schemes are implemented, i.e., the Wesely scheme and a “dynamic” scheme, in the latter of which dry-deposition parameterization is coupled with photosynthesis and vegetation dynamics. It is demonstrated that, when the dynamic scheme is applied, the simulated O3 dry-deposition velocities vd and their stomatal and cuticular portions, as well as the total O3 fluxes Ft, are larger overall; vd and Ft are 2–3 times more sensitive to the SM changes due to the data assimilation (DA). Further, through case studies at two forested sites with different soil types and hydrological regimes, we highlight that, applying the Community Land Model type of SM factor controlling stomatal resistance (i.e., β factor) scheme in replacement of the Noah-type β factor scheme reduced the v(d) sensitivity to SM changes by ∼75 % at one site, while it doubled this sensitivity at the other site. Referring to multiple evaluation datasets, which may be associated with variable extents of uncertainty, the model performance of vegetation, surface fluxes, weather, and surface O3 concentrations shows mixed responses to the DA, some of which display land cover dependency. Finally, using model-derived concentration- and flux-based policy-relevant O3 metrics as well as their matching exposure–response functions, the relative biomass/crop yield losses for several types of vegetation/crops are estimated to be within a wide range of 1 %–17 %. Their sensitivities to the model's dry-deposition scheme and the implementation of SM DA are discussed.

Min Huang↗

Coupled Ablation, Heat Conduction, Pyrolysis, Shape Change and Spallation of the Galileo Probe

The Galileo probe enters the atmosphere of Jupiter in December 1995. This paper presents numerical methodology and detailed results of our final pre-impact calculations for the heat shield response. The calculations are performed using a highly modified version of a viscous shock layer code with massive radiation coupled with a surface thermochemical ablation and spallation model and with the transient in-depth thermal response of the charring and ablating heat shield. The flowfield is quasi-steady along the trajectory, but the heat shield thermal response is dynamic. Each surface node of the VSL grid is coupled with a one-dimensional thermal response calculation. The thermal solver includes heat conduction, pyrolysis, and grid movement owing to surface recession. Initial conditions for the heat shield temperature and density were obtained from the high altitude rarefied-flow calculations of Haas and Milos. Galileo probe surface temperature, shape, mass flux, and element flux are all determined as functions of time along the trajectory with spallation varied parametrically. The calculations also estimate the in-depth density and temperature profiles for the heat shield. All this information is required to determine the time-dependent vehicle mass and drag coefficient which are necessary inputs for the atmospheric reconstruction experiment on board the probe.

Milos, Frank S.↗

Shortwave feedbacks and El Nino-Southern Oscillation: Forced ocean coupled ocean-atmosphere experiments

Changes in tropical sea surface temperature (SST) can produce changes in cloudiness that modify incoming solar shortwave (SW) radiation, which in turn affects SST. The effects of this negative feedback on Pacific interannual variability are examined in forced ocean model and hybrid coupled ocean-atmosphere model simulations. Two empirical schemes are used to model the large-scale, low-frequency response of surface SW to SST anomalies. The first scheme attempts to account for the nonlocal nature of the atmospheric response to SST based patterns of covariability analyzed through singular value decomposition. In the observations the primary coupled mode of variability is composed of a SW anomaly in the central Pacific that covaries with anomalous SST in the eastern Pacific. This is applied in the model as a nonlocal feedback. The second scheme examines the effects of a purely local feedback with a spatially varying coefficient of magnitude chosen similar to the first scheme. In almost all cases the second scheme behaved similarly to the first, presumably because the correlation scale of SST is large enough for El Nino-Southern Oscillation (ENSO) dynamics that there is little sensitivity to the local approximation in the SW feedback. In simulations forced by time series of observed wind stress the SW feedback induced very minor SST damping. Results for a simpified heat budget analysis showed that while the SW feedback increased the local heat flux damping on SST, it also induced a mean shallowing of the mixed layer. The resulting changes in both the local mean vertical temperature gradient and the zonal velocity response to the wind stress acted to oppose the local heat flux damping effects. When the observed SW anomalies were applied to forced simulations, the simulated SST anomalies were modified as expected, and agreement with observed SST improved. In coupled simulations the SW feedbacks had greater impact than in the case of specified stress. The main effects were to decrease the magnitude of the warm and cold SST anomalies in the central Pacific, while leaving the pattern and evolution of ENSO anomalies essentially unchanged elsewhere. The SW feedbacks thus produce a modest improvement of the model ENSO SST pattern compared with observations, although they tended to shorten the period of the model ENSO cycle. Overall the results suggest that large-scale SW feedbacks are of quantitative importance to simulating some aspects of the ENSO cycle but are not critical to the overall occurrence of the phenomenon.

Waliser, Duane E.↗

Analysis of an electrohydraulic aircraft control surface servo and comparison with test results

An analysis of an electrohydraulic aircraft control-surface system is made in which the system is modeled as a lumped, two-mass, spring-coupled system controlled by a servo valve. Both linear and nonlinear models are developed, and the effects of hinge-moment loading are included. Transfer functions of the system and approximate literal factors of the transfer functions for several cases are presented. The damping action of dynamic pressure feedback is analyzed. Comparisons of the model responses with results from tests made on a highly resonant rudder control-surface servo indicate the adequacy of the model. The effects of variations in hinge-moment loading are illustrated.

Edwards, J. W.↗

Development of an Effective Finite-rate Oxidation Model for NuSil-coated Charred Carbon Preform Ablators

Recently, a detailed effective finite-rate surface chemistry model was developed for the oxidation of FiberForm [1] using the molecular beam experimental data of Poovathingal et al., [2]. FiberForm is the major building block of the thermal protection system (TPS) material Phenolic Impregnated Carbon Ablator (PICA), commonly used by NASA. The surface chemistry model consists of detailed surface reaction mechanisms such as adsorption, desorption, and several types of Langmuir-Hinshelwood (LH) reactions to characterize the oxygen-carbon interactions at the surface. This model provides excellent agreement with the experimental data for oxidation product compositions and corresponding translational energy distributions. Further, an effective oxidation model was constructed that captures the equivalent interaction of oxygen inside the microstructure (including multiple surface collisions) via a modified reactivity for a smooth wall boundary condition. This enables the use of this model directly into Computational Fluid Dynamics (CFD) codes and Material Response (MR) codes to accurately simulate the gas–surface interactions within FiberForm without using the detailed micro-structure. However, this effective model is valid only for the virgin FiberForm. As the TPS material undergoes ablation, the preform carbon burns and turns into char. In addition, PICA is also coated with a protective silicon coating called NuSil for the purpose of mitigating the spread of phenolic dust, and limit contamination during clean room operations. Fig. 1 shows a X-ray microtomography image of a charred NuSil-coated PICA. The NuSil layer introduces species containing silicon into the product mixture. Further, the reactivity of the carbon within the char layer is different from the virgin FiberForm as shown in Fig. 2. Recently a new set of molecular beam experiments were performed on this NuSil-coated TPS material [3] with the same type of oxygen beam used in the previous experiments. Using the latest experimental data, the previously developed effective model will be extended to account for the charred carbon as well as the NuSil coating. Finally, this new effective model with three phases – preform carbon, char, and NuSil; will be compared and validated against the experimental product compositions.

K Swaminathan Gopalan↗

Variability of Phenology and Fluxes of Water and Carbon with Observed and Simulated Soil Moisture in the Ent Terrestrial Biosphere Model (Ent TBM Version 1.0.1.0.0)

The Ent Terrestrial Biosphere Model (Ent TBM) is a mixed-canopy dynamic global vegetation model developed specifically for coupling with land surface hydrology and general circulation models (GCMs). This study describes the leaf phenology submodel implemented in the Ent TBM version 1.0.1.0.0 coupled to the carbon allocation scheme of the Ecosystem Demography (ED) model. The phenology submodel adopts a combination of responses to temperature (growing degree days and frost hardening), soil moisture (linearity of stress with relative saturation) and radiation (light length). Growth of leaves, sapwood, fine roots, stem wood and coarse roots is updated on a daily basis. We evaluate the performance in reproducing observed leaf seasonal growth as well as water and carbon fluxes for four plant functional types at five Fluxnet sites, with both observed and prognostic hydrology, and observed and prognostic seasonal leaf area index. The phenology submodel is able to capture the timing and magnitude of leaf-out and senescence for temperate broadleaf deciduous forest (Harvard Forest and Morgan- Monroe State Forest, US), C3 annual grassland (Vaira Ranch, US) and California oak savanna (Tonzi Ranch, US). For evergreen needleleaf forest (Hyytiäla, Finland), the phenology submodel captures the effect of frost hardening of photosynthetic capacity on seasonal fluxes and leaf area. We address the importance of customizing parameter sets of vegetation soil moisture stress response to the particular land surface hydrology scheme. We identify model deficiencies that reveal important dynamics and parameter needs.

Soil moisture↗

The GLAS physical inversion method for analysis of TIROS N data

The Goddard Laboratory for Atmospheric Sciences approach to temperature sounding from satellite observations is fundamentally different from that of NESS in that heavy reliance is placed on the ability to model accurately the instrumental response to atmospheric and surface conditions, while no use is made of statistical relationships between satellite observations and atmospheric temperature profiles. The method involves starting with a guess set of atmospheric and surface conditions from which expected brightness temperatures for the satellite observations are computed. Then, iterative relaxation of atmospheric and surface conditions is performed according to the difference between observed and computed brightness temperatures until sufficient agreement is reached. In the absence of sufficient agreement, no retrieval is produced for that location.

Susskind, J.↗

Mechanisms of interannual ocean-atmosphere interactions

Results of an investigation to study the interaction between ocean and atmosphere on the annual to decadal time scale are reported. Separate studies of the ocean response to atmospheric forcing, and of atmospheric response to ocean forcing were also conducted. The main findings are the dynamics of sea surface temperature anomalies, the role of short time scale weather fluctuations in the seasonal cycle of the upper ocean variability, and the planetary wave response to sea surface temperature anomalies. A numerical model of the ocean atmosphere continent system and a two layer quasi-geostropic ocean model is discussed.

Frankignoul, C.↗

On transient rheology and glacial isostasy

The effect of transient creep on the inference of long-term mantle viscosity is investigated using theoretical predictions from self-gravitating, layered earth models with Maxwell, Burgers' body, and standard linear solid rheologies. The interaction between transient and steady-state rheologies is studied. The responses of the standard linear solid and Burgers' body models to transient creep in the entire mantle, and of the Burgers' body and Maxwell models to creep in the lower mantle are described. The models' responses are examined in terms of the surface displacement, free air gravity anomaly, wander of the rotation pole, and the secular variation of the degree 2 zonal coefficient of the earth's gravitational potential field. The data reveal that transient creep cannot operate throughout the entire mantle.

Yuen, David A.↗

MOM3D method of moments code theory manual

MOM3D is a FORTRAN algorithm that solves Maxwell's equations as expressed via the electric field integral equation for the electromagnetic response of open or closed three dimensional surfaces modeled with triangle patches. Two joined triangles (couples) form the vector current unknowns for the surface. Boundary conditions are for perfectly conducting or resistive surfaces. The impedance matrix represents the fundamental electromagnetic interaction of the body with itself. A variety of electromagnetic analysis options are possible once the impedance matrix is computed including backscatter radar cross section (RCS), bistatic RCS, antenna pattern prediction for user specified body voltage excitation ports, RCS image projection showing RCS scattering center locations, surface currents excited on the body as induced by specified plane wave excitation, and near field computation for the electric field on or near the body.

Shaeffer, John F.↗

Variable-Complexity Multidisciplinary Optimization on Parallel Computers

This report covers work conducted under grant NAG1-1562 for the NASA High Performance Computing and Communications Program (HPCCP) from December 7, 1993, to December 31, 1997. The objective of the research was to develop new multidisciplinary design optimization (MDO) techniques which exploit parallel computing to reduce the computational burden of aircraft MDO. The design of the High-Speed Civil Transport (HSCT) air-craft was selected as a test case to demonstrate the utility of our MDO methods. The three major tasks of this research grant included: development of parallel multipoint approximation methods for the aerodynamic design of the HSCT, use of parallel multipoint approximation methods for structural optimization of the HSCT, mathematical and algorithmic development including support in the integration of parallel computation for items (1) and (2). These tasks have been accomplished with the development of a response surface methodology that incorporates multi-fidelity models. For the aerodynamic design we were able to optimize with up to 20 design variables using hundreds of expensive Euler analyses together with thousands of inexpensive linear theory simulations. We have thereby demonstrated the application of CFD to a large aerodynamic design problem. For the predicting structural weight we were able to combine hundreds of structural optimizations of refined finite element models with thousands of optimizations based on coarse models. Computations have been carried out on the Intel Paragon with up to 128 nodes. The parallel computation allowed us to perform combined aerodynamic-structural optimization using state of the art models of a complex aircraft configurations.

Grossman, Bernard↗

Parametric Studies of Flow Separation using Air Injection

Boundary Layer separation causes the airfoil to stall and therefore imposes dramatic performance degradation on the airfoil. In recent years, flow separation control has been one of the active research areas in the field of aerodynamics due to its promising performance improvements on the lifting device. These active flow separation control techniques include steady and unsteady air injection as well as suction on the airfoil surface etc. This paper will be focusing on the steady and unsteady air injection on the airfoil. Although wind tunnel experiments revealed that the performance improvements on the airfoil using injection techniques, the details of how the key variables such as air injection slot geometry and air injection angle etc impact the effectiveness of flow separation control via air injection has not been studied. A parametric study of both steady and unsteady air injection active flow control will be the main objective for this summer. For steady injection, the key variables include the slot geometry, orientation, spacing, air injection velocity as well as the injection angle. For unsteady injection, the injection frequency will also be investigated. Key metrics such as lift coefficient, drag coefficient, total pressure loss and total injection mass will be used to measure the effectiveness of the control technique. A design of experiments using the Box-Behnken Design is set up in order to determine how each of the variables affects each of the key metrics. Design of experiment is used so that the number of experimental runs will be at minimum and still be able to predict which variables are the key contributors to the responses. The experiments will then be conducted in the 1ft by 1ft wind tunnel according to the design of experiment settings. The data obtained from the experiments will be imported into JMP, statistical software, to generate sets of response surface equations which represent the statistical empirical model for each of the metrics as a function of the key variables. Next, the variables such as the slot geometry can be optimized using the build-in optimizer within JMP. Finally, a wind tunnel testing will be conducted using the optimized slot geometry and other key variables to verify the empirical statistical model. The long term goal for this effort is to assess the impacts of active flow control using air injection at system level as one of the task plan included in the NASAs URETI program with Georgia Institute of Technology.

Zhang, Wei↗