Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical optimal control”

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 721 records · Page 40

Space Laboratory on a Table Top: A Next Generative ECLSS design and diagnostic tool

This paper describes the development plan for a comprehensive research and diagnostic tool for aspects of advanced life support systems in space-based laboratories. Specifically it aims to build a high fidelity tabletop model that can be used for the purpose of risk mitigation, failure mode analysis, contamination tracking, and testing reliability. We envision a comprehensive approach involving experimental work coupled with numerical simulation to develop this diagnostic tool. It envisions a 10% scale transparent model of a space platform such as the International Space Station that operates with water or a specific matched index of refraction liquid as the working fluid. This allows the scaling of a 10 ft x 10 ft x 10 ft room with air flow to 1 ft x 1 ft x 1 ft tabletop model with water/liquid flow. Dynamic similitude for this length scale dictates model velocities to be 67% of full-scale and thereby the time scale of the model to represent 15% of the full- scale system; meaning identical processes in the model are completed in 15% of the full- scale-time. The use of an index matching fluid (fluid that matches the refractive index of cast acrylic, the model material) allows making the entire model (with complex internal geometry) transparent and hence conducive to non-intrusive optical diagnostics. So using such a system one can test environment control parameters such as core flows (axial flows), cross flows (from registers and diffusers), potential problem areas such as flow short circuits, inadequate oxygen content, build up of other gases beyond desirable levels, test mixing processes within the system at local nodes or compartments and assess the overall system performance. The system allows quantitative measurements of contaminants introduced in the system and allows testing and optimizing the tracking process and removal of contaminants. The envisaged system will be modular and hence flexible for quick configuration change and subsequent testing. The data and inferences from the tests will allow for improvements in the development and design of next generation life support systems and configurations. Preliminary experimental and modeling work in this area will be presented. This involves testing of a single inlet-exit model with detailed 3-D flow visualization and quantitative diagnostics and computational modeling of the system.

Ramachandran, N.↗

A Machine Learning Approach to Improve Air Traffic Management Initiatives

Collaborating closely with commercial air carriers and related organizations, the Federal Aviation Administration(FAA) regulates air traffic and ensures the safety and efficiency of air operations. Air traffic controllers make strategic decisions, such as delaying, rerouting, or canceling flights, partly based on guidance provided by the FAA’s Air TrafficControl System Command Center (ATCSCC). The guidance includes, among other things, control measures known asTraffic Management Initiatives (TMIs) designed to enhance safety and improve operational efficiency. TMIs play a crucial role in managing the demand and capacity within the U.S. National Airspace System (NAS). Two major TMIs that are routinely used (primarily to mitigate the adverse effects of bad weather) are Ground Delay Programs (GDPs) andGround Stops (GSs). In a GDP, flights destined for airports facing thunderstorm activity experience delays at their origin airports. This proactive approach minimizes the risk of routing aircraft through hazardous weather conditions and also replaces (fuel burning) airborne delays with ground delays. In a GS, a temporary restriction is imposed on the departure or arrival of aircraft at a specific airport or within a designated airspace. Although other TMIs (e.g., miles-in-trail) are also implemented as part of (air) traffic flow management in the NAS, the focus of this work is on GDPs and GSs. Since TMIs, by design, lead to flight delays or cancellations, it is crucial to put in place the right set of parameters(e.g., scope and duration of the GDP). For example, when the end time of a GDP extends beyond what is necessary, it imposes unnecessary delays on departing flights. This situation could occur as a result of inaccurate prediction of the(required) duration of the GDP based on the weather forecast. On the other hand, if a GDP ends prematurely before the underlying capacity constraints are resolved at the destination airport, it may result in airborne holding. The delicate balance lies in matching the termination of the GDP precisely with the resolution of capacity constraints, avoiding both the imposition of unnecessary ground delays and the need for airborne holding due to premature program termination.Failing to specify the right parameters for TMIs also leads to flight delays, creating a significant obstacle in managing the increasing traffic volumes causing increased work load for the controllers. To address this issue, we propose the integration of Machine Learning (ML) models in the traffic flow management(TFM) pipeline. In current operations, decisions are made by human experts based on extensive training, historical patterns, available traffic and weather data. Since we have an abundance of data from past events that tell us the likely impact of various TMIs, by ingesting historical data, properly trained ML models can offer valuable insights and aid human decision-making. With the FAA increasingly exploring advanced analytics, ML emerges as a focal point for enhancing TFM within the National Airspace System (NAS). As a first step, this study aims to provide traffic controllers with decision-making support for the issuance and adjustment of TMIs. Data analytics and machine learning have been previously employed to address some of the challenges associated with TMIs. Numerous studies have concentrated on various facets of TMI issuance, exploring factors influencing TMI parameters, including arrival rate, airport capacity, and delay prediction. For example, using weather forecasts, several statistical methods were used to produce probabilistic capacity profiles which in conjunction with deterministic models provided insights into the GDP planning process [1–4]. The downside of using deterministic models is that they rely on fixed inputs and predetermined rules, which lack the ability to account for the inherent uncertainty and variability present in real-world scenarios. In a separate series of studies, researchers aimed to predict the occurrences of GDPs and GSs. The majority of these studies utilized various supervised learning methods, including Decision Trees, Naive Bayes, Support VectorMachines, and Random Forests to analyze the influence of weather conditions and arrival demand on TMI incidents[5–8]. However, these studies primarily focused on predicting the incidence of TMIs without explicitly addressing the scope of TMIs, including their duration and their geographical coverage. Furthermore, the emphasis of these studies was largely on GDPs, given their higher frequency and longer duration when compared to GSs. A limited number of studies focused on predicting the parameters of TMIs, specifically addressing their duration and extent. In one such study focusing on optimizing the TMI parameters at San Francisco International Airport (SFO),the authors utilized a probabilistic forecast of fog [9]. They simulated various capacity scenarios based on the (fog)burn-off forecasts, selecting GDP parameters that minimized airborne and overall ground delays. However, this approach exclusively emphasizes stratus (fog) burn-off as the primary determinant of GDP and GS, neglecting other influential factors like severe weather events, runway closures, lower capacity than traffic demand, and other important variables. Given the complexity of predicting the TMI and determining its scope, we seek a more holistic approach. We aim to consider all significant factors that could impact TMIs and their parameters. What sets this research apart is the fusion of all data sources relevant to the issuance and adjustment of TMIs and it represents the first comprehensive attempt to optimize TMIs in this manner. Since this comprehensive solution involves various aspects, we break down the problem into smaller components and input all parameters into a unified model called the “TMI Adjuster”. Figure 1 shows the overall framework and the list of datasets used in each model. The objective of the TMI Adjuster module is to deliver reliable, consistent and expedited recommendations for the progression, adjustment, and termination of TMIs. The ML solution entails developing a pipeline capable of predicting the necessity of a TMI (e.g., GS or GDP) along with its various parameters. For example, in the case of a GS, this includes the scope of the GS either in terms of distance from the destination airport or based on pre-defined airspace sectors. Here, scope refers to those regions and departing airports that are subject to the GS. In this paper, we concentrate on the issuance of GSs in the three major airports in the New York area — LaGuardia(LGA), John F. Kennedy International (JFK), and Newark Liberty International (EWR). We fuse traffic, weather and other relevant aviation data from years 2017 to 2019 to train and validate the ML models. In particular, we use the following datasets: •Terminal Aerodrome Forecast (TAF): meteorological forecasts specific to each airport, issued four times a day, covering predefined time periods. •TMI data: includes all GSs and GDPs along with their respective parameters. •Aviation System Performance Metrics (ASPM): includes traffic related data such as aircraft delays, arrival, and departure rates. •Notices to Airmen (NOTAMs): utilized to extract runway closure data and manage interdependencies between terminals in close proximity. •Flight cancellation data •Airspace Flow Programs (AFP): includes information on flight airborne holdings caused by TMIs. The data preprocessing entails transforming ASPM, TMI, AFP, NOTAMs, and weather data into an hourly format and consolidating all datasets by merging them based on date and time as the primary key. The TMI Adjuster framework comprises two parallel models: one dedicated to GS and a second model focused on GDP. As previously mentioned, our specific focus is on the GS model as a multi-classification problem. In this framework, each data point of the GS model input summarizes ten hours of data. Specifically, the data loader for the GS model generates the input and output of the model as follows: at a given time step, the input includes the actual traffic, weather, and TMI data from the two-hour window before the time step, alongside the weather forecast and scheduled traffic for the next 8 hours starting from the time step. Based on this information, the output of the GS model for each time interval consists of three dimensions. The first dimension represents a binary decision on whether there should be a GS in place for the next hour or not. The second dimension is related to the scope of the GS in the United States, and the third dimension is related to the scope of the GS in Canada (i.e., to determine if the GS impacts airports in Canada).One of the challenges with TMI modeling is the sparsity of TMI events, particularly regarding its scope. To address this challenge in the scope of the GS model output, we implement grouping. The GS scope for the US region is defined based on a list of centers that should be included when the GS is in place. With 20 centers in the US, we utilized historical data to group them into 4 categories. In particular, we summarized our historical data in a graph format where nodes represent centers, and link weights are defined based on the co-occurrence of centers in the scope parameter ofTMIs. By identified strongly connected components in this graph, we were able to partition the centers into four groups. We consider two model structures for the GS Model. Firstly, a hierarchical classification model [10], where the human decision-making for a GS is of hierarchical nature. The decision-maker first decides whether there is a need fora GS, and if the answer is yes, determines the scope. A hierarchical classification model organizes the problem into a class hierarchy, typically a tree or a Directed Acyclic Graph (DAG) structure, and considers the dependency of the decision in the previous step to the next component [10]. Here, we employ the local classifier per level approach, which involves training one multi-class classifier for each level of the class hierarchy. The second structure is the independent structure. In this setting, as the name suggests, we do not consider the dependency of the decisions in the different dimensions of the output of the model. Instead, for each dimension, we train a multi-class classifier independently. Table 1 summarizes GS model statistics for training, validation and testing. The table documents the effect of limiting data to the time steps when there was actually a TMI in place or when a TMI had just terminated. This resulted in a more balanced distribution of the GS class(GS positive class)versus “No GS”(GS negative class), which might help the training process. While JFK and LGA follow very similar distributions, with 40% and 42% GS positive class respectively, EWR has proportionally fewer GS incidents at 28%. Our subsequent phase involves evaluating the performance of both hierarchical structure and independent structure using different state-of-the-art multi-class classifier models such as Random Forest, Decision Trees, K-nearest Neighbors, and Logistic Regression and forecast the duration and scope of the GSs.

Farzan Masrour Shalmani↗

Using XR for Improving Scientific Discovery With Numerical Weather Models

Earth science (ES) digital twins will help us understand the complex interactions and interrelationships that make up our Earth system and the impacts of earth science phenomena on it. Our work addresses two underdeveloped areas in current ES digital twin work: improving the understanding and interaction with ES model outputs by using Virtual and Mixed Reality (XR) tools and improving the non-intuitive mapping of continuous ES natural phenomena to gridded reference frames in current numerical models. Traditionally, scientists working on ES view and analyze the results of calculated or measured observables with static 1-dimensional (1D), 2D or 3D plots displayed on flat computer screens or paper. Using such limited mediums, it can be very difficult to identify, track and understand the evolution of key features due to poor viewing angles and the nature of flat computer screens. In addition, numerical models, such as the NASA Goddard Earth Observing System (GEOS) ES model, are almost exclusively formulated, visualized and analyzed in an Eulerian reference frame with fixed grid points in space and time. However, ES phenomena such as convective clouds, hurricanes and wildfire smoke plumes are visualized and analyzed in a Lagrangian reference frame: therefore it is often difficult and unnatural to understand these phenomena in relation to each other, visualized either in an Eulerian or Lagrangian context. In 3D visualizations, data generally takes one of three forms: gridded (e.g., voxelized) data, where space is divided into regions; point clouds, where data is represented as a set of points; and meshes, where objects are rendered as surfaces composed of small polygons (usually triangles). A gridded, Eulerian reference frame has been the default representation for the 2D visual analysis of atmospheric data in part because the numerical methods used to generate atmospheric model data in the first place use a gridded approach, with equations defining the relationships between the physical variables in each of a grid's cells across successive timesteps. In our work, we are particularly interested in data from GEOS. Another reason why gridded representations tend to be used for visualizing data from such models is because trajectories are difficult to interpret from representations on 2D surfaces, due to line-of-sight ambiguity. Instead of a fixed grid from GEOS, we embed a trajectory model to simulate particles' movement throughout a GEOS run. We then ingest these particle trajectories as animated point clouds with a NASA open source XR toolkit, the Mixed Reality Exploration Toolkit (MRET), and merge GEOS data with ES phenomena data onto one combined visualization that the user can intuitively interact with. Efficient rendering of arbitrarily large point clouds is an ongoing challenge being addressed by the computer science community, with the GPU-based optimizations and efficient GPU memory utilization a common theme of recent advances, especially for XR, where sustained high frame rate is mandatory to save the user from suffering due to simulation sickness. In this work, we describe and evaluate our progress in choosing and implementing appropriate methods for rendering arbitrarily large point clouds within MRET for XR. While tracking the XR headset enables the immersion of a user within a 3D scene of a data visualization, tracking of XR handheld controllers or user’s hands enables us to implement intuitive user interactions with the visualized datasets. Conventional tools require a user working with an ES visualization to conduct many interactions to commit their intended selections or manipulations with a visualized dataset; for example to specify a set of points in 3D space. Doing so in a 2D flat screen interface has traditionally required specifying a set of points in three distinct 2D coordinate systems (XY, XZ, and YZ), which is cumbersome. In other scientific domains, it has been shown that specifying or selecting a location or volume in XR using handheld controllers or tracked hands allows for greater speed and accuracy. We anticipate the same will hold true for atmospheric data, and we will share initial results of measuring the utility of such an interface. Notably, as the data being visualized is generated by GEOS as a prediction based on initial conditions, an intended application of our tool is to serve as part of an iterative feedback loop. Through XR, a scientist will review and manipulate a GEOS model run, modifying the conditions as needed to do subsequent runs of GEOS. Thereby, XR-based improvements to speed and accuracy of 3D tagging of points minimizes the effort required by both the scientist and the computer cluster conducting the necessary calculations.

Thomas Grubb↗

Stellar Occultation Studies of the Solar System

Earth-based observations of stellar occultations provide extremely high spatial resolution for bodies in the outer solar system, about 10,000 times better than that of traditional imaging observations. Stellar occultation data can be used to establish the structure of atmospheres and rings of solar system bodies at high spatial resolution. Airborne occultation observations are particularly effective, since the controlled mobility of the observing platform allows the observer to fly within the optimum part of the occultation shadow for most events that are visible from Earth. Airborne observations are carried out above any clouds and are nearly free of scintillation noise from the Earth's atmosphere. KAO occultation observations resulted in the first detection of gravity waves in the Martian atmosphere, discovery of the Uranian rings, the first detection of Pluto's atmosphere, the first Earth-based investigations of Triton's atmosphere, and the discovery of narrow jets from Chiron's nucleus. The first SOFIA occultation opportunity will be an investigation of Pluto's atmospheric structure in November, 2002, and will resolve a problem that has lingered since the KAO discovery observation fourteen years earlier. We plan to continue our successful airborne occultation program with the greatly enhanced capability provided by SOFIA. We propose here to replace our KAO occultation photometer with one having twice the throughput, half the noise, a somewhat wider wavelength range, four times the field of view, and ten times the frame rate to optimize its performance and to capitalize on the larger collecting area offered by SOFIA. It will also allow for simultaneous visible and IR occultation observations, greatly enriching the results that we can obtain from occultations. We call this new imaging occultation photometer HOPI (High-speed Occultation Photometer and Imager). HOPI will provide a signal-to-noise ratio two to four times that of our present photometer for a given event, will permit useful observations of the more numerous events involving fainter stars, and will allow higher time resolution observations to be made. HOPI's characteristics also make it an ideal instrument for initial and ongoing evaluations of the SOFIA telescope's performance.

Elliot, James L.↗

Stellar Occultation Studies of the Solar System

Earth-based observations of stellar occultations provide extremely high spatial resolution for bodies in the outer solar system, about 10,000 times better than that of traditional imaging observations. Stellar occultation data can be used to establish the structure of atmospheres and rings of solar system bodies at high spatial resolution. Airborne occultation observations are particularly effective, since the controlled mobility of the observing platform allows the observer to fly within the optimum part of the occultation shadow for most events that are visible from Earth. Airborne observations are carried out above any clouds and are nearly free of scintillation noise from the Earth's atmosphere. KAO occultation observations resulted in the first detection of gravity waves in the Martian atmosphere, discovery of the Uranian rings, the first detection of Pluto's atmosphere, the first Earth-based investigations of Triton's atmosphere, and the discovery of narrow jets from Chiron's nucleus. The first SOFIA occultation opportunity will be an investigation of Pluto's atmospheric structure in November, 2002, and will resolve a problem that has lingered since the KAO discovery observation fourteen years earlier. We plan to continue our successful airborne occultation program with the greatly enhanced capability provided by SOFIA. We propose here to replace our KAO occultation photometer with one having twice the throughput, half the noise, a somewhat wider wavelength range, four times the field of view, and ten times the frame rate to optimize its performance and to capitalize on the larger collecting area offered by SOFIA. It will also allow for simultaneous visible and IR occultation observations, greatly enriching the results that we can obtain from occultations. We call this new imaging occultation photometer HOPI (High-speed Occultation Photometer and Imager). HOPI will provide a signal-to-noise ratio two to four times that of our present photometer for a given event, will permit useful observations of the more numerous events involving fainter stars, and will allow higher time resolution observations to be made. HOPI's characteristics also make it an ideal instrument for initial and ongoing evaluations of the SOFIA telescope's performance.

Elliot, James L.↗

Nonlinear rescaling of control values simplifies fuzzy control

Traditional control theory is well-developed mainly for linear control situations. In non-linear cases there is no general method of generating a good control, so we have to rely on the ability of the experts (operators) to control them. If we want to automate their control, we must acquire their knowledge and translate it into a precise control strategy. The experts' knowledge is usually represented in non-numeric terms, namely, in terms of uncertain statements of the type 'if the obstacle is straight ahead, the distance to it is small, and the velocity of the car is medium, press the brakes hard'. Fuzzy control is a methodology that translates such statements into precise formulas for control. The necessary first step of this strategy consists of assigning membership functions to all the terms that the expert uses in his rules (in our sample phrase these words are 'small', 'medium', and 'hard'). The appropriate choice of a membership function can drastically improve the quality of a fuzzy control. In the simplest cases, we can take the functions whose domains have equally spaced endpoints. Because of that, many software packages for fuzzy control are based on this choice of membership functions. This choice is not very efficient in more complicated cases. Therefore, methods have been developed that use neural networks or generic algorithms to 'tune' membership functions. But this tuning takes lots of time (for example, several thousands iterations are typical for neural networks). In some cases there are evident physical reasons why equally space domains do not work: e.g., if the control variable u is always positive (i.e., if we control temperature in a reactor), then negative values (that are generated by equal spacing) simply make no sense. In this case it sounds reasonable to choose another scale u' = f(u) to represent u, so that equal spacing will work fine for u'. In the present paper we formulate the problem of finding the best rescaling function, solve this problem, and show (on a real-life example) that after an optimal rescaling, the un-tuned fuzzy control can be as good as the best state-of-art traditional non-linear controls.

Vanlangingham, H.↗

Investigation of Condensing Ice Heat Exchangers for MTSA Technology Development

Metabolic heat regenerated Temperature Swing Adsorption (MTSA) technology is being developed for thermal, carbon dioxide (CO2) and humidity control for a Portable Life Support Subsystem (PLSS). Metabolically-produced CO2 present in the ventilation gas of a PLSS is collected using a CO2selective adsorbent via temperature swing adsorption. The temperature swing is initiated through cooling to well below metabolic temperatures. Cooling is achieved with a sublimation heat exchanger using water or liquid carbon dioxide (LCO2) expanded below sublimation temperature when exposed to low pressure or vacuum. Subsequent super heated vapor, as well as additional coolant, is used to further cool the astronaut. The temperature swing on the adsorbent is then completed by warming the adsorbent with a separate condensing ice heat exchanger (CIHX) using metabolic heat from moist ventilation gas. The condensed humidity in the ventilation gas is recycled at the habitat. The water condensation from the ventilation gas is a significant heat transfer mechanism for the warming of the adsorbent bed because it represents as much as half of the energy potential in the moist ventilation gas. Designing a heat exchanger to efficiently transfer this energy to the adsorbent bed and allow the collection of the water is a challenge since the CIHX will operate in a temperature range from 210K to 280K. The ventilation gas moisture will first freeze and then thaw, sometimes existing in three phases simultaneously. A NASA Small Business Innovative Research (SBIR) Phase 1 contract was performed to investigate condensing and icing as applied to MTSA to enable higher fidelity modeling and assess the impact of geometry variables on CIHX performance for future CIHX design optimization. Specifically, a design tool was created using analytical relations to explore the complex, interdependent design space of a condensing ice heat exchanger. Numerous variables were identified as having nontrivial contributions to performance such as hydraulic diameter, heat exchanger effectiveness, ventilation gas mass flow rate and surface roughness. Using this tool, four test articles were designed and manufactured to map to a full MTSA subassembly (the adsorbent bed, the sublimation heat exchanger for cooling and the condensing ice heat exchanger for warming). The design mapping considered impacts due to CIHX geometry as well as subassembly impacts such as thermal mass and thermal resistance through the adsorbent bed. The test articles were tested at simulated PLSS ventilation loop temperature, moisture content and subambient pressure. Ice accumulation and melting were observed. Data and test observations were analyzed to identify drivers of the condensing ice heat exchanger performance. This paper will discuss the analytical models, the test article designs, and testing procedures. Testing issues will be discussed to better describe data and share lessons learned. Data analysis and subsequent conclusions will be presented.

Padilla, Sebastian↗

Active Control of the Forced and Transient Response of a Finite Beam

When studying structural vibrations resulting from a concentrated source, many structures may be modelled as a finite beam excited by a point source. The theoretical limit on cancelling the resulting beam vibrations by utilizing another point source as an active controller is explored. Three different types of excitation are considered, harmonic, random, and transient. In each case, a cost function is defined and minimized for numerous parameter variations. For the case of harmonic excitation, the cost function is obtained by integrating the mean squared displacement over a region of the beam in which control is desired. A controller is then found to minimize this cost function in the control interval. The control interval and controller location are continuously varied for several frequencies of excitation. The results show that control over the entire beam length is possible only when the excitation frequency is near a resonant frequency of the beam, but control over a subregion may be obtained even between resonant frequencies at the cost of increasing the vibration outside of the control region. For random excitation, the cost function is realized by integrating the expected value of the displacement squared over the interval of the beam in which control is desired. This is shown to yield the identical cost function as obtained by integrating the cost function for harmonic excitation over all excitation frequencies. As a result, it is always possible to reduce the cost function for random excitation whether controlling the entire beam or just a subregion, without ever increasing the vibration outside the region in which control is desired. The last type of excitation considered is a single, transient pulse. A cost function representative of the beam vibration is obtained by integrating the transient displacement squared over a region of the beam and over all time. The form of the controller is chosen a priori as either one or two delayed pulses. Delays constrain the controller to be causal. The best possible control is then examined while varying the region of control and the controller location. It is found that control is always possible using either one or two control pulses. The two pulse controller gives better performance than a single pulse controller, but finding the optimal delay time for the additional controllers increases as the square of the number of control pulses.

Post, John Theodore↗

Modified Interior Distance Functions (Theory and Methods)

In this paper we introduced and developed the theory of Modified Interior Distance Functions (MIDF's). The MIDF is a Classical Lagrangian (CL) for a constrained optimization problem which is equivalent to the initial one and can be obtained from the latter by monotone transformation both the objective function and constraints. In contrast to the Interior Distance Functions (IDF's), which played a fundamental role in Interior Point Methods (IPM's), the MIDF's are defined on an extended feasible set and along with center, have two extra tools, which control the computational process: the barrier parameter and the vector of Lagrange multipliers. The extra tools allow to attach to the MEDF's very important properties of Augmented Lagrangeans. One can consider the MIDFs as Interior Augmented Lagrangeans. It makes MIDF's similar in spirit to Modified Barrier Functions (MBF's), although there is a fundamental difference between them both in theory and methods. Based on MIDF's theory, Modified Center Methods (MCM's) have been developed and analyzed. The MCM's find an unconstrained minimizer in primal space and update the Lagrange multipliers, while both the center and the barrier parameter can be fixed or updated at each step. The MCM's convergence was investigated, and their rate of convergence was estimated. The extension of the feasible set and the special role of the Lagrange multipliers allow to develop MCM's, which produce, in case of nondegenerate constrained optimization, a primal and dual sequences that converge to the primal-dual solutions with linear rate, even when both the center and the barrier parameter are fixed. Moreover, every Lagrange multipliers update shrinks the distance to the primal dual solution by a factor 0 less than gamma less than 1 which can be made as small as one wants by choosing a fixed interior point as a 'center' and a fixed but large enough barrier parameter. The numericai realization of MCM leads to the Newton MCM (NMCM). The approximation for the primal minimizer one finds by Newton Method followed by the Lagrange multipliers update. Due to the MCM convergence, when both the center and the barrier parameter are fixed, the condition of the MDF Hessism and the neighborhood of the primal ninimizer where Newton method is 'well' defined remains stable. It contributes to both the complexity and the numerical stability of the NMCM.

Polyak, Roman A.↗

Super-optimal CO2 reduces seed yield but not vegetative growth in wheat

Although terrestrial atmospheric CO2 levels will not reach 1000 micromoles mol-1 (0.1%) for decades, CO2 levels in growth chambers and greenhouses routinely exceed that concentration. CO2 levels in life support systems in space can exceed 10000 micromoles mol-1(1%). Numerous studies have examined CO2 effects up to 1000 micromoles mol-1, but biochemical measurements indicate that the beneficial effects of CO2 can continue beyond this concentration. We studied the effects of near-optimal (approximately 1200 micromoles mol-1) and super-optimal CO2 levels (2400 micromoles mol-1) on yield of two cultivars of hydroponically grown wheat (Triticum aestivum L.) in 12 trials in growth chambers. Increasing CO2 from sub-optimal to near-optimal (350-1200 micromoles mol-1) increased vegetative growth by 25% and seed yield by 15% in both cultivars. Yield increases were primarily the result of an increased number of heads per square meter. Further elevation of CO2 to 2500 micromoles mol-1 reduced seed yield by 22% (P < 0.001) in cv. Veery-10 and by 15% (P < 0.001) in cv. USU-Apogee. Super-optimal CO2 did not decrease the number of heads per square meter, but reduced seeds per head by 10% and mass per seed by 11%. The toxic effect of CO2 was similar over a range of light levels from half to full sunlight. Subsequent trials revealed that super-optimal CO2 during the interval between 2 wk before and after anthesis mimicked the effect of constant super-optimal CO2. Furthermore, near-optimal CO2 during the same interval mimicked the effect of constant near-optimal CO2. Nutrient concentration of leaves and heads was not affected by CO2. These results suggest that super-optimal CO2 inhibits some process that occurs near the time of seed set resulting in decreased seed set, seed mass, and yield.

NASA Discipline Life Support Systems↗

Successful Completion of the JWST OGSE2 Cryogenic Test at JSC Chamber-A While Managing Numerous Challenges

The James Webb Space Telescope (JWST) Optical Telescope Element (OTE) assembly is the largest optically stable infrared-optimized telescope currently being manufactured and assembled, and scheduled for launch in 2018. The JWST OTE, including the primary mirrors, secondary mirror, and the Aft Optics Subsystems (AOS) are designed to be passively cooled and operate at near 45 degrees Kelvin. Due to the size of its large sunshield in relation to existing test facilities, JWST cannot be optically or thermally tested as a complete observatory-level system at flight temperatures. As a result, the telescope portion along with its instrument complement will be tested as a single unit very late in the program, and on the program schedule critical path. To mitigate schedule risks, a set of 'pathfinder' cryogenic tests will be performed to reduce program risks by demonstrating the optical testing capabilities of the facility, characterizing telescope thermal performance, and allowing project personnel to learn valuable testing lessons off-line. This paper describes the 'pathfinder' cryogenic test program, focusing on the recently completed second test in the series called the Optical Ground Support Equipment 2 (OGSE2) test. The JWST OGSE2 was successfully completed within the allocated project schedule while faced with numerous conflicting thermal requirements during cool-down to the final cryogenic operational temperatures, and during warm-up after the cryo-stable optical tests. The challenges include developing a pre-test cool-down and warm-up profiles without a reliable method to predict the thermal behaviors in a rarified helium environment, and managing the test article hardware safety driven by the project Limits and Constraints (L&C's). Furthermore, OGSE2 test included the time critical Aft Optics Subsystem (AOS), a part of the flight Optical Telescope Element that would need to be placed back in the overall telescope assembly integrations. The OGSE2 test requirements included the strict adherence of the project contamination controls due to the presence of the contamination sensitive flight optical elements. The test operations required close coordination of numerous personnel while they being exposed and trained for the 'final' combined OTE and instrument cryo-test in 2017. This paper will also encompass the OGSE2 thermal data look-back review.

JWST OGSE2 Cryogenic Test↗

Instability Wave–Streak Interactions in a High Mach Number Boundary Layer at Flight Conditions

The interaction of stationary streaks undergoing nonmodal growth with modally unstable instability waves in a hypersonic boundary-layer flow is studied using numerical computations. The geometry and flow conditions are selected to match a relevant trajectory location from the ascent phase of the HIFiRE-1 flight experiment; namely, a 7 degree half-angle, circular cone with 2.5 mm nose radius, freestream Mach number equal to 5.30, unit Reynolds number equal to 13.42 m -1 , and wall-to-adiabatic temperature ratio of approximately 0.35 over most of the vehicle. This paper investigates the nonlinear evolution of initially linear optimal disturbances that evolve into finite-amplitude streaks, followed by an analysis of the modal instability characteristics of the perturbed, streaky boundary-layer flow. The investigation is performed with stationary direct numerical simulations (DNS) and plane-marching parabolized stability equations (PSE), in conjunction with partial-differential-equation-based planar eigenvalue analysis. The overall effect of streaks is to reduce the peak amplification factors of instability waves, indicating a possible downstream shift in the onset of laminar-turbulent transition. The present study confirms previous findings that the mean flow distortion of the nonlinear streak perturbation reduces the amplification rates of the Mack-mode instability. More importantly, however, the present results demonstrate that the spanwise varying component of the streak can produce a larger effect on the Mack-mode amplification. The study with selected azimuthal wavenumbers for the stationary streaks reveals that a wavenumber of approximately 1.4 times larger than the optimal wavenumber is more effective in stabilizing the planar Mack-mode instabilities. In the absence of unstable first-mode waves for the present cold-wall condition, transition onset is expected to be delayed until the peak streak amplitude increases to nearly 35 percent of the freestream velocity, when intrinsic instabilities of the boundary-layer streaks begin to dominate the transition process. For streak amplitudes below that limit a significant net stabilization is achieved, yielding a potential transition delay that can exceed 100 percent of the length of the laminar region in the uncontrolled case.

Boundary layer stability↗

Scaling Properties of Particle Density Fields Formed in Simulated Turbulent Flows

Direct numerical simulations (DNS) of particle concentrations in fully developed 3D turbulence were carried out in order to study the nonuniform structure of the particle density field. Three steady-state turbulent fluid fields with Taylor microscale Reynolds numbers (Re(sub lambda)) of 40, 80 and 140 were generated by solving the Navier-Stokes equations with pseudospectral methods. Large scale forcing was used to drive the turbulence and maintain temporal stationarity. The response of the particles to the fluid was parameterized by the particle Stokes number St, defined as the ratio of the particle's stopping time to the mean period of eddies on the Kolmogorov scale (eta). In this paper, we consider only passive particles optimally coupled to these eddies (St approx. = 1) because of their tendency to concentrate more than particles with lesser or greater St values. The trajectories of up to 70 million particles were tracked in the equilibrated turbulent flows until the particle concentration field reached a statistically stationary state. The nonuniform structure of the concentration fields was characterized by the multifractal singularity spectrum, f(alpha), derived from measures obtained after binning particles into cells ranging from 2(eta) to 15(eta) in size. We observed strong systematic variations of f(alpha) across this scale range in all three simulations and conclude that the particle concentration field is not statistically self similar across the scale range explored. However, spectra obtained at the 2(eta), 4(eta), and 8(eta) scales of each flow case were found to be qualitatively similar. This result suggests that the local structure of the particle concentration field may be flow-Independent. The singularity spectra found for 2n-sized cells were used to predict concentration distributions in good agreement with those obtained directly from the particle data. This Singularity spectrum has a shape similar to the analogous spectrum derived for the inertial-range energy dissipation fields of experimental turbulent flows at Re(sub lambda) = 110 and 1100. Based on this agreement, and the expectation that both dissipation and particle concentration are controlled by the same cascade process, we hypothesize that singularity spectra similar to the ones found in this work provide a good characterization of the spatially averaged statistical properties of preferentially concentrated particles in higher Re(sub lambda) turbulent flows.

Hogan, Robert C.↗

Navigation Strategy for the Mars 2001 Lander Mission

The Mars Surveyor Program (MSP) is an ongoing series of missions designed to robotically study, map and search for signs of life on the planet Mars. The MSP 2001 project will advance the effort by sending an orbiter, a lander and a rover to the red planet in the 2001 opportunity. Each vehicle will carry a science payload that will Investigate the Martian environment on both a global and on a local scale. Although this mission will not directly search for signs of life, or cache samples to be returned to Earth, it will demonstrate certain enabling technologies that will be utilized by the future Mars Sample Return missions. One technology that is needed for the Sample Return mission is the capability to place a vehicle on the surface within several kilometers of the targeted landing site. The MSP'01 Lander will take the first major step towards this type of precision landing at Mars. Significant reduction of the landed footprint will be achieved through two technology advances. The first, and most dramatic, is hypersonic aeromaneuvering; the second is improved approach navigation. As a result, the guided entry will produce in a footprint that is only tens of kilometers, which is an order of magnitude improvement over the Pathfinder and Mars Polar Lander ballistic entries. This reduction will significantly enhance scientific return by enabling the potential selection of otherwise unreachable landing sites with unique geologic interest and public appeal. A landed footprint reduction from hundreds to tens of kilometers is also a milestone on the path towards human exploration of Mars, where the desire is to place multiple vehicles within several hundred meters of the planned landing site. Hypersonic aeromaneuvering is an extension of the atmospheric flight goals of the previous landed missions, Pathfinder and Mars Polar Lander (MPL), that utilizes aerodynamic lift and an autonomous guidance algorithm while in the upper atmosphere. The onboard guidance algorithm will control the direction of the lift vector, via bank angle modulation, to keep the vehicle on the desired trajectory. While numerous autonomous guidance algorithms have been developed for use during hypersonic flight at Earth, this will be the first flight of an autonomously directed lifting entry vehicle at Mars. However, without sufficient control and knowledge of the atmospheric entry conditions, the guidance algorithm will not perform effectively. The goal of the interplanetary navigation strategy is to deliver the spacecraft to the desired entry condition with sufficient accuracy and knowledge to enable satisfactory guidance algorithm performance. Specifically, the entry flight path angle must not exceed 0.27 deg. to a 3 sigma confidence level. Entry errors will contribute directly to the size of the landed footprint and the most significant component is entry flight path angle. The size of the entry corridor is limited on the shallow side by integrated heating constraints, and on the steep side by deceleration (g-load) and terminal descent propellant. In order to meet this tight constraint it is necessary to place a targeting maneuver seven hours prior to the time of entry. At this time the trajectory knowledge will be quite accurate, and the effects of maneuver execution errors will be small. The drawback is that entry accuracy is dependent on the success of this final late maneuver. Because propulsive maneuvers are critical events, it is desirable to minimize their occurrence and provide the flight team with as much response time as possible in the event of a spacecraft fault. A mission critical maneuver at Entry - 7 hours does not provide much fault tolerance, and it is desirable to provide a strategy that minimizes reliance on this maneuver. This paper will focus on the Improvements in interplanetary navigation that will decrease entry errors and will reduce the landed footprint, even in the absence of aeromaneuvering. The easiest to take advantage of are Improvements In the knowledge of the Mars ephemeris and gravity field due to the MGS and MSP'98 missions. Improvements In data collection and reduction techniques such as "precislon ranging' and near-simultaneous tracking will also be utilized. In addition to precise trajectory control, a robust strategy for communications and flight operations must also be demonstrated. The result Is a navigation and communications strategy on approach that utilizes optimal maneuver placement to take advantage of trajectory knowledge, minimizes risk for the flight operations team, is responsive to spacecraft hardware limitations, and achieves the entry corridor. The MSP2001 mission Is managed at JPL under the auspices of the Mars Exploration Directorate. The spacecraft flight elements are built and managed by Lockheed-Martin Astronautics in Denver, Colorado.

Mase, Robert A.↗

An Optimization Approach to Support Science Decision Making for Lunar Surface Exploration

Introduction: Scientific exploration is one of the three pillars of NASA’s Moon2Mars architecture, with crew surface extra vehicular activities (EVA) serving a critical enabling function. Development of surface EVA operational planning and execution, specifically integrating science and flight control teams (FCT), is currently being explored through analog scenarios. This integration, exercised, for example, through the Joint EVA and Hu-man Surface Mobility Test Team (JETT), allows for science input on EVA activities in near real-time through a Science Evaluation Room (SER), or Arte-mis science backroom, which integrates with the broader FCT through the Science Officer. The SER works within the FCT to support dynamic EVA planning in response to changes in operational constraints as well as science opportunities and re-prioritization, increasing the mission science return and accelerating the accomplishment of the Moon2Mars science objectives. The SER works within the FCT to provide recommendations to traverse execution in near real-time. One challenge is the requirement to deliver SER inputs to the FCT on operationally relevant timelines. Failure to do so may result in suboptimal execution of science exploration EVAs or even loss of key science objectives. To close this gap, we present a network optimization tool to allow the SER to provide rapid input to the FCT in response to changes in operational constraints or science opportunities. Inputs are predicated on approved science objectives, and clear rationale must be provided to the FCT for any requested change. Accordingly, this tool incorporates the Science Traceability Matrix (STM), SER prioritization scheme, and station characterization and action planning with operational constraints such as duration, traverse speed, and distance to maximize science objectives based on SER priorities, consistent with FCT operational requirements. Method: As a proof of concept, we used an existing linear programing software package used to simulate optimal routes through cellular metabolism. We built a Demonstrative Model with three STM objectives and four stations on a region of the Moon. The objectives were given an arbitrary prioritization and mapped to the stations through four possible crew actions. (Figs. 1 and 2). This station to STM mapping is consistent with the method used by the JETT5 Science Team to develop analog surface EVA science planning. We used a grid system with the landing site at the origin and the four stations placed across the positive x,y quadrant. Actions were assigned to each station and the accomplishment of those actions resulted in a numerical “reward” based on the ability of that action to achieve science objectives. The aggregate reward from each individual STM objective contributes to a global score (Science Yield), weighted by its priority. Operational constraints included a requirement to start and end at the landing site, 5 minutes each for initial station characterization and “clean up,” and variable total EVA time, traverse rate (fixed to 0.5 meters per second in our example), and time to perform each action (10, 5, 7, and 15 min for actions 1, 2, 3, and 4, respectively). Additional constraints and variables will be added in the future (e.g., sample mass, number of stations, traverse route constraints, illumination). Optimization. We converted the connections (arcs) between these stations (nodes) into a mixed integer linear programming optimization problem (arcs = constraints, nodes = variables) with the objective to maximize Science Yield. For any action, the Science Yield is equal to the relevance of that action to an STM objective [3, 2, and 1 point(s) for High, Med., and Low relevance, respectively], multiplied by the STM Objective Priority [3, 2, and 1 point(s) for High, Med., and Low priority, respectively]. This resulted in a model that computes the optimal station and action combination to maximize the Science Yield. These weightings can be adjusted by the SER as desired. Results: We explored three test cases for the Demonstrative Model. First, we set the maximum EVA duration to 120 minutes and computed the optimal route (Fig. 3A). The model suggested per-forming Actions 1 and 2 at Station P01, followed by Actions 1 and 2 at Station P02, and finally Actions 1 and 3 at Station P04 before returning to the Landing Site. Second, we adjusted the STM Objective Priori-ty order and computed the new optimal route (Fig. 3B). Under this situation, the model suggested per-forming all Actions at Station P02 followed by all Actions at Station P03. The previous test cases were relevant to SER planning activities. Next, we explored providing mid-EVA replanning input to the FCT. Scenario: While executing the Route in Fig. 3A the crew finishes at Station P01 and FCT decides that the EVA needs to finish in 45 minutes back at the Landing Site. FCT asks SER to recommend changes to the plan to accommodate this operation-al change. Using the model and incorporating these new constraints (start at Station P01, max. time of 45 min), the model suggested performing Actions 2 and 4 at Station P03 (Fig. 4), requiring 41 minutes to complete and return to the Landing Site. Interestingly, Station 3 was not part of the original route. Using the model, we determined the EVA would need 66 minutes, instead of 45, in order for the original Station P04 to yield a larger Science Yield than Station P03. The parametrization and simulation was per-formed in less than a minute, demonstrating the operational relevance of the approach. Future Efforts: The results from the Demonstrative Model suggest this tool can accelerate SER decision making on operationally relevant timelines. Use in analog activities, such as JETT5 or follow-ons, which have over a dozen stations for a crew to explore and over a dozen actions per station, will provide needed validation of the utility of this tool for planning EVAs, replanning mid-EVA, or planning follow-on EVAs based on previous results. Further integration with FCT execution monitoring tools may provide additional efficiency gains, al-lowing rapid and iterative exploration of operation-al and science decision space by the FCT and SER.

Science Operations↗

Development of a Ground Multi-Mission Low-Cost Optical Terminal (LCOT) for Free-Space Optical Communications

Once confined to the realm of laboratory experiments and theoretical papers, space-based laser communications (lasercomm) are on the verge of achieving mainstream status. Organizations from Facebook to NASA, and missions from cubesats to Orion are employing lasercomm to achieve gigabit communication speeds at mass and power requirements lower than that of traditional radio frequency (RF) methods. Since first demonstrating free-space optical communications services with Lunar Laser Communications Demonstration (LLCD) in 2013, NASA has invested in developing optical communications technologies and capabilities to enhancing its space communications networks. Along with evolving optical space terminals, NASA is also developing lasercomm grounds stations capable of meeting the rapidly increasing data volume demands of upcoming missions, from low-earth to lunar orbits and beyond and integrating these advanced capabilities into its Near-Space and Deep Space Networks. To meet this emerging need, the Low-Cost Optical Terminal (LCOT) project at NASA’s Goddard Space Flight Center (GSFC) is designing, building and validating a prototype for a flexible, multi-mission, and economical optical ground terminal that could be used as a blueprint for a global network of optical ground stations, capable of supporting a wide variety of missions. To date a major impediment to widespread adoption of laser communication has been the lack of an existing ground network infrastructure. A mission that wishes to take advantage of laser communication not only needs to invest in an optical space terminal, but it must also finance the creation of ground terminals to receive the downlink signal. This adds significant additional cost. Missions that do decide to incur the cost of financing a network of ground terminals end up building highly specialized optical receivers that are operable as receivers for that specific mission only. Significant Non-Recurring Engineering (NRE) cost is invested to build highly specialized one-of-a-kind ground terminals that go into storage after that particular mission is over. This is not an economical approach and does nothing to grow the number of optical ground stations available to future missions. In essence each mission that wants to take advantage of the benefits of lasercom has to start from scratch to provide a ground terminal network to support it. As long as this is the case, the cost for using laser communications will be too high for most missions to consider. LCOT intends to close this gap in technology by designing and developing a standard optical ground terminal design that is flexible enough to serve as a receiver for a wide range of future missions – a ground terminal that can be quickly reconfigured to receive downlinks at different wavelengths using different signal formats. Not only does LCOT have the industry-building objectives of utilizing commercial-off-the-shelf (COTS) components to the maximum extent possible, but also spurs the commercial development of other necessary lasercomm components not currently offered by industry. Finally, LCOT will give NASA scientists and engineers a facility where they can gain real-world experience with optical communications. It will give engineers a cost-effective way to try out new concepts and processes by providing the infrastructure for such testing. In this way it is hoped LCOT will serve as a stimulus for innovation in optical communications and speed its widespread adoption by future missions. The LCOT is comprised of five subsystems: Free-Space Optical, Transceiver, Amplifier, Monitor and Control, and Observatory Infrastructure. In August 2021, the LCOT team installed a 70 cm telescope, developed by Planewave Instruments that was optimized for optical communications. Free-Space Optical subsystem comprises of the telescope and its associated hardware, including a transmitter optical assembly, wide field cameras, two optical benches, and an adaptive optics subsystem. The transmit optical assembly, a unique concept design, is a cluster of four functionally independent transmit subassemblies located on the receive telescope. In addition to receiving optical signals and directing the expanded beam with high precision to the space terminal, it also performs tracking functions. The transmit optical assembly will support operations from Low Earth Orbit (LEO) through lunar and will be used as a template for industry manufacturing. The Optical Infrastructure subsystem is responsible for providing environmentally controlled shelters for LCOT equipment and various other systems. To maintain the safety and proper functionality of the telescope, a 16 ft Astrohaven clamshell dome procured which provides all-sky coverage without the need to rotate the dome. Additionally, Atmospheric Monitoring Assembly (AMA) will be part of optical infrastructure subsystem to ensure accurate performance of the LCOT. Like existing optical ground stations, LCOT will measure standard weather station parameters, infrared all sky image of cloud cover, and cloud height. LCOT, however, adds requirements for measuring night time seeing and, in the future, daytime seeing. Unlike other optical ground terminals, the LCOT is transceiver agnostic; user transceivers may be duplex transceivers, standalone receivers, or standalone transmitters with or without acquisition beacon functionality. As such, the LCOT project accommodates testing with external customer transceivers in a flexible manor, further complimenting its intended multi-mission goals. Another unique component of LCOT is the use of a new amplifier technology – the Very Large Mode Area (VLMA) amplifiers. This new technology allows LCOT to avoid the issues faced by previous laser communications ground terminals, gives users more flexibility and modular capability, and is capable of reaching an order of magnitude higher peak power than traditional High Power Optical Amplifiers (HPOA). One drawback of the VLMA HPOA approach is that the amplified light is output into free-space. The solution developed by LCOT is an optics train that couples the output of the VLMA amplifier into a short fiber for transport to the transmit telescopes with high efficiency. Like many of the LCOT components, a set of detailed manufacturing drawings have been created for the optics train to allow any machine shop with a multi-axis Computer Numerical Control (CNC) machine to fabricate the piece parts from commonly available materials. In line with the goals of LCOT, the monitor and control functions are developed as a modular and flexible system with the ability to support future hardware or algorithm changes, minimizing disruptions. A main priority of development in the Monitor and Control Subsystem (MCS) is the safety monitor system.

laser communication↗

Development of a Ground Multi-mission Low Cost Optical Terminal(LCOT) for Free-Space Optical Communication

Once confined to the realm of laboratory experiments and theoretical papers, space-based laser communications (lasercomm) are on the verge of achieving mainstream status. Organizations from Facebook to NASA, and missions from cubesats to Orion are employing lasercomm to achieve gigabit communication speeds at mass and power requirements lower than that of traditional radio frequency (RF) methods. Since first demonstrating free-space optical communications services with Lunar Laser Communications Demonstration (LLCD) in 2013, NASA has invested in developing optical communications technologies and capabilities to enhancing its space communications networks. Along with evolving optical space terminals, NASA is also developing lasercomm grounds stations capable of meeting the rapidly increasing data volume demands of upcoming missions, from low-earth to lunar orbits and beyond and integrating these advanced capabilities into its Near-Space and Deep Space Networks. To meet this emerging need, the Low-Cost Optical Terminal (LCOT) project at NASA’s Goddard Space Flight Center (GSFC) is designing, building and validating a prototype for a flexible, multi-mission, and economical optical ground terminal that could be used as a blueprint for a global network of optical ground stations, capable of supporting a wide variety of missions. To date a major impediment to widespread adoption of laser communication has been the lack of an existing ground network infrastructure. A mission that wishes to take advantage of laser communication not only needs to invest in an optical space terminal, but it must also finance the creation of ground terminals to receive the downlink signal. This adds significant additional cost. Missions that do decide to incur the cost of financing a network of ground terminals end up building highly specialized optical receivers that are operable as receivers for that specific mission only. Significant Non-Recurring Engineering (NRE) cost is invested to build highly specialized one-of-a-kind ground terminals that go into storage after that particular mission is over. This is not an economical approach and does nothing to grow the number of optical ground stations available to future missions. In essence each mission that wants to take advantage of the benefits of lasercom has to start from scratch to provide a ground terminal network to support it. As long as this is the case, the cost for using laser communications will be too high for most missions to consider. LCOT intends to close this gap in technology by designing and developing a standard optical ground terminal design that is flexible enough to serve as a receiver for a wide range of future missions – a ground terminal that can be quickly reconfigured to receive downlinks at different wavelengths using different signal formats. Not only does LCOT have the industry-building objectives of utilizing commercial-off-the-shelf (COTS) components to the maximum extent possible, but also spurs the commercial development of other necessary lasercomm components not currently offered by industry. Finally, LCOT will give NASA scientists and engineers a facility where they can gain real-world experience with optical communications. It will give engineers a cost-effective way to try out new concepts and processes by providing the infrastructure for such testing. In this way it is hoped LCOT will serve as a stimulus for innovation in optical communications and speed its widespread adoption by future missions. The LCOT is comprised of five subsystems: Free-Space Optical, Transceiver, Amplifier, Monitor and Control, and Observatory Infrastructure. In August 2021, the LCOT team installed a 70 cm telescope, developed by Planewave Instruments that was optimized for optical communications. Free-Space Optical subsystem comprises of the telescope and its associated hardware, including a transmitter optical assembly, wide field cameras, two optical benches, and an adaptive optics subsystem. The transmit optical assembly, a unique concept design, is a cluster of four functionally independent transmit subassemblies located on the receive telescope. In addition to receiving optical signals and directing the expanded beam with high precision to the space terminal, it also performs tracking functions. The transmit optical assembly will support operations from Low Earth Orbit (LEO) through lunar and will be used as a template for industry manufacturing. The Optical Infrastructure subsystem is responsible for providing environmentally controlled shelters for LCOT equipment and various other systems. To maintain the safety and proper functionality of the telescope, a 16 ft Astrohaven clamshell dome procured which provides all-sky coverage without the need to rotate the dome. Additionally, Atmospheric Monitoring Assembly (AMA) will be part of optical infrastructure subsystem to ensure accurate performance of the LCOT. Like existing optical ground stations, LCOT will measure standard weather station parameters, infrared all sky image of cloud cover, and cloud height. LCOT, however, adds requirements for measuring night time seeing and, in the future, daytime seeing. Unlike other optical ground terminals, the LCOT is transceiver agnostic; user transceivers may be duplex transceivers, standalone receivers, or standalone transmitters with or without acquisition beacon functionality. As such, the LCOT project accommodates testing with external customer transceivers in a flexible manor, further complimenting its intended multi-mission goals. Another unique component of LCOT is the use of a new amplifier technology – the Very Large Mode Area (VLMA) amplifiers. This new technology allows LCOT to avoid the issues faced by previous laser communications ground terminals, gives users more flexibility and modular capability, and is capable of reaching an order of magnitude higher peak power than traditional High Power Optical Amplifiers (HPOA). One drawback of the VLMA HPOA approach is that the amplified light is output into free-space. The solution developed by LCOT is an optics train that couples the output of the VLMA amplifier into a short fiber for transport to the transmit telescopes with high efficiency. Like many of the LCOT components, a set of detailed manufacturing drawings have been created for the optics train to allow any machine shop with a multi-axis Computer Numerical Control (CNC) machine to fabricate the piece parts from commonly available materials. In line with the goals of LCOT, the monitor and control functions are developed as a modular and flexible system with the ability to support future hardware or algorithm changes, minimizing disruptions. A main priority of development in the Monitor and Control Subsystem (MCS) is the safety monitor system.

Haleh Safavi↗

Analysis and Design of Variable Stiffness Composite Cylinders

An investigation of the possible performance improvements of thin circular cylindrical shells through the use of the variable stiffness concept is presented. The variable stiffness concept implies that the stiffness parameters change spatially throughout the structure. This situation is achieved mainly through the use of curvilinear fibers within a fiber-reinforced composite laminate, though the possibility of thickness variations and discrete stiffening elements is also allowed. These three mechanisms are incorporated into the constitutive laws for thin shells through the use of Classical Lamination Theory. The existence of stiffness variation within the structure warrants a formulation of the static equilibrium equations from the most basic principles. The governing equations include sufficient detail to correctly model several types of nonlinearity, including the formation of a nonlinear shell boundary layer as well as the Brazier effect due to nonlinear bending of long cylinders. Stress analysis and initial buckling estimates are formulated for a general variable stiffness cylinder. Results and comparisons for several simplifications of these highly complex governing equations are presented so that the ensuing numerical solutions are considered reliable and efficient enough for in-depth optimization studies. Four distinct cases of loading and stiffness variation are chosen to investigate possible areas of improvement that the variable stiffness concept may offer over traditional constant stiffness and/or stiffened structures. The initial investigation deals with the simplest solution for cylindrical shells in which all quantities are constant around the circumference of the cylinder. This axisymmetric case includes a stiffness variation exclusively in the axial direction, and the only pertinent loading scenarios include constant loads of axial compression, pressure, and torsion. The results for these cases indicate that little improvement over traditional laminates exists through the use of curvilinear fibers, mainly due to the presence of a weak link area within the stiffness variation that limits the ultimate load that the structure can withstand. Rigorous optimization studies reveal that even though slight increases in the critical loads can be produced for designs with an arbitrary variation of the fiber orientation angle, the improvements are not significant when compared to traditional design techniques that utilize ring stiffeners and frames. The second problem that is studied involves arbitrary loading of a cylinder with a stiffness variation that changes only in the circumferential direction. The end effects of the cylinder are ignored, so that the problem takes the form of an analysis of a cross-section for a short cylinder segment. Various load cases including axial compression, pressure, torsion, bending, and transverse shear forces are investigated. It is found that the most significant improvements in load-carrying capability exist for cases which involve loads that also vary around the circumference of the shell, namely bending and shear forces. The stiffness variation of the optimal designs contribute to the increased performance in two ways: lowering the stresses in the critical areas through redistribution of the stresses; and providing a relatively stiff region that alters the buckling behavior of the structure. These results lead to an in-depth optimization study involving weight optimization of a fuselage structure subjected to typical design constraints. Comparisons of the curvilinear fiber format to traditional stiffened structures constructed of isotropic and composite materials are included. It is found that standard variable stiffness designs are quite comparable in terms of weight and load-carrying capability yet offer the added advantage of tailorability of distinct regions of the structure that experience drastically different loading conditions. The last two problems presented in this work involve the nonlinear phenomenon of long tubes under bending. Though this scenario is not as applicable to fuselage structures as the previous problems, the mechanisms that produce the nonlinear effect are ideally suited to be controlled by the variable stiffness concept. This is due to the fact that the dominating influence for long cylinders under bending is the ovalization of the cross-section, which is governed mainly by the stiffness parameters of the cylindrical shell. Possible improvement of the critical buckling moments for these structures is investigated using either a circumferential or axial stiffness variation. For the circumferential case involving infinite length cylinders, it is found that slight improvements can be observed by designing structures that resist the cross-sectional deformation yet do not detract from the buckling resistance at the critical location. The results also indicate that buckling behavior is extremely dependent on cylinder length. This effect is most easily seen in the solution of finite length cylinders under bending that contain an axial stiffness variation. For these structures, the only mechanism that exhibits improved response are those that effectively shorten the length of the cylinder, thus reducing the cross-sectional deformation due to the forced restraint at the ends. It was found that the use of curvilinear fibers was not able to achieve this effect in sufficient degree to resist the deformation, but that ring stiffeners produced the desired response admirably. Thus, it is shown that the variable stiffness concept is most effective at improving the bending response of long cylinders through the use of a circumferential stiffness variation.

Tatting, Brian F.↗