Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “G&C”

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 37 records · Page 2

Expanding the genetic toolkit: adenine and cytosine base editors for gene disruption in Aspergillus niger

Despite revolutionizing fungal genetic engineering, conventional CRISPR/Cas9-mediated knockouts rely on DNA double-strand breaks (DSBs), which can cause unwanted insertions and deletions, chromosomal abnormalities, and cytotoxicity. Base editors such as adenine base editors (ABEs), which convert A‧T to G‧C, and cytosine base editors (CBEs), which convert C‧G to T‧A, offer a safer alternative by enabling predictable, target-specific single-nucleotide changes without introducing DSBs. To overcome the limitations of traditional genome editing in filamentous fungi, we developed efficient base-editing systems in Aspergillus niger . For the first time, we constructed an ABE in A. niger , achieving up to 80% editing efficiency and inducing predictable A-to-G mutations at the intended intron sites, disrupting gene function through mRNA mis-splicing. We also developed a highly efficient CBE system, capable of introducing premature stop codons with 50–100% efficiency. To broaden the editing scope, we implemented a Cas9-NG variant recognizing a relaxed PAM sequence requiring only a single guanine (G), enabling editing at start codons and splice sites. Leveraging this expanded scope, we established gene disruption approaches by targeting start codons via ABE-mediated A-to-G conversions (ATG-to-GTG and ATG-to-ACG) and CBE-mediated C-to-T conversion (ATG-to-ATA). Additionally, our base-editing systems enable multiplex gRNA delivery and marker-free editing of multiple genes. Collectively, the scope-expanding strategies increase the number of genes targetable for disruption by base-editing in A. niger by 26.3% and enable near-complete coverage of 96% of the coding genes. Overall, this work demonstrates the potential of ABE and CBE systems as versatile, efficient, and safer alternatives to DSBs-based gene disruption in filamentous fungi.

Aspergillus↗

Watson-Crick Base-Pairing Requirements for ssDNA Recognition and Processing in Replication-Initiating HUH Endonucleases

Replication-initiating HUH endonucleases (Reps) are sequence-specific nucleases that cleave and rejoin single-stranded DNA (ssDNA) during rolling-circle replication. These functions are mediated by covalent linkage of the Rep to its substrate post cleavage. Here, we describe the structures of the endonuclease domain from the Muscovy duck circovirus Rep in complex with its cognate ssDNA 10-mer with and without manganese in the active site. Structural and functional analyses demonstrate that divalent cations play both catalytic and structural roles in Reps by polarizing and positioning their substrate. Further structural comparisons highlight the importance of an intramolecular substrate Watson-Crick (WC) base pairing between the -4 and +1 positions. Subsequent kinetic and functional analyses demonstrate a functional dependency on WC base pairing between these positions regardless of the pair’s identity (i.e., A·T, T·A, G·C, or C·G), highlighting a structural specificity for substrate interaction. Finally, considering how well WC swaps were tolerated in vitro, we sought to determine to what extent the canonical -4T·+1A pairing is conserved in circular Rep-encoding single-stranded DNA viruses and found evidence of noncanonical pairings in a minority of these genomes. Altogether, our data suggest that substrate intramolecular WC base pairing is a universal requirement for separation and reunion of ssDNA in Reps.

59 BASIC BIOLOGICAL SCIENCES↗

Different tertiary interactions create the same important 3D features in a distinct flavivirus xrRNA

During infection by a flavivirus (FV), cells accumulate noncoding subgenomic flavivirus RNAs (sfRNAs) that interfere with several antiviral pathways. These sfRNAs are formed by structured RNA elements in the 3' untranslated region (UTR) of the viral genomic RNA, which block the progression of host cell exoribonucleases that have targeted the viral RNA. Previous work on these exoribonuclease-resistant RNAs (xrRNAs) from mosquito-borne FVs revealed a specific three-dimensional fold with a unique topology in which a ring-like structure protectively encircles the 5' end of the xrRNA. Conserved nucleotides make specific tertiary interactions that support this fold. Examination of more divergent FVs reveals differences in their 3' UTR sequences, raising the question of whether they contain xrRNAs and if so, how they fold. To answer this, we demonstrated the presence of an authentic xrRNA in the 3' UTR of the Tamana bat virus (TABV) and solved its structure by X-ray crystallography. The structure reveals conserved features from previously characterized xrRNAs, but in the TABV version these features are created through a novel set of tertiary interactions not previously seen in xrRNAs. This includes two important A–C interactions, four distinct backbone kinks, several ordered Mg 2+ ions, and a C + –G–C base triple. The discovery that the same overall architecture can be achieved by very different sequences and interactions in distantly related flaviviruses provides insight into the diversity of this type of RNA and will inform searches for undiscovered xrRNAs in viruses and beyond.

59 BASIC BIOLOGICAL SCIENCES↗

Green’s function coupled cluster simulation of the near-valence ionizations of DNA-fragments

Accurate description of the ionization process in DNA is crucial to the understanding of the DNA damage under exposure to ionizing radiation and the exploration of the potential application of DNA strands in nanoelectronics. In this work, by employing our recently developed Green’s function coupled-cluster library on supercomputing facilities, we have studied the spectral functions of several guanine–cytosine (G–C) base pair structures ([G–C]n, n = 1–3) for the first time in a relatively broad near-valence regime ([−25.0, −5.0] eV) in the coupled-cluster with singles and doubles level. Our focus is to give a preliminary many-body coupled-cluster understanding and guideline of the vertical ionization energy (VIE), spectral profile, and ionization feature changes of these systems as the system size expands in this near-valence regime. The results show that, as the system size expands, even though the lowest VIEs keep decreasing, the changes of spectral function profile and the relative peak positions get unexpectedly smaller. Further analysis of the ionized states associated with the most intensive peak in the spectral functions reveals non-negligible |2h, 1p⟩’s in the ionized wave functions of the considered G–C base pair systems. The leading |2h, 1p⟩’s associated with the main ionizations from the cytosine part of the G–C base pairs feature a transition from the intra-base-pair cytosine π → π* excitation to the inter-base-pair electron excitation as the size of G–C base pairs expands, which also indicates the minimum quantum region in the many-body calculations of DNA systems.

Peng, Bo (ORCID:0000000242267294)↗

Three-Dimensional Analysis of Nuclear Size, Shape and Displacement in Clover Root Cap Statocytes from Space and a Clinostat

Under normal (l-g) conditions the statocytes of root caps have a characteristic polarity with the nucleus in tight association with the proximal cell wall; but, in altered gravity environments including microgravity (mu-g) and the clinostat (c-g) movement of the nucleus away from the proximal cell wall is not uncommon. To further understand the cause of gravity-dependent nuclear displacement in statocytes, three-dimensional cell reconstruction techniques were used to precisely measure the volumes, shapes, and positions of nuclei in white clover (Trifolium repens) flown in space and rotated on a clinostat. Seeds were germinated and grown for 72 hours aboard the Space Shuttle (STS-63) in the Fluid Processing Apparatus (BioServe Space Technologies, Univ. of Colorado, Boulder). Clinorotation experiments were performed on a two-axis clinostat (BioServe). Computer reconstruction of selected groups of statocytes were made from serial sections (0.5 microns thick) using the ROSS (Reconstruction Of Serial Sections) software package (Biocomputation Center, NASA Ames Research Center). Nuclei were significantly displaced from the tops of cells in mu-g (4.2 +/- 1.0 microns) and c-g (4.9 +/- 1.4 microns) when compared to l-g controls (3.4 +/- 0.8 gm); but, nuclear volume (113 +/- 36 cu microns, 127 +/- 32 cu microns and 125 +/- 28 cu microns for l-g, mu-g and c-g respectively) and the ratio of nuclear volume to cell volume (4.310.7%, 4.211.0% and 4.911.4% respectively) were not significantly dependent on gravity treatment (ANOVA; alpha = 0.05). Three-dimensional analysis of nuclear shape and proximity to the cell wall, however, showed that nuclei from l-g controls appeared ellipsoidal while those from space and the clinostat were more spherically shaped. This change in nuclear shape may be responsible for its displacement under altered gravity conditions. Since the cytoskeleton is known to affect nuclear polarity in root cap statocytes, those same cytoskeletal elements could also control nuclear shape. This alteration in nuclear shape and position in mu-g and c-g when compared to l-g may lead to functional differences in the gravity signaling systems of plants subjected to altered gravity environments.

Smith, J.D.↗

Conundrum for the free energy of a holonomous gluonic plasma at cubic order

We compute the term ~ g 3 in the free energy for a S U ( N ) gauge theory with nonzero holonomy at nonzero temperature. If the holonomy is generated kinematically by the introduction of gauge invariant sources coupled to Polyakov loops, the contribution of charged (off-diagonal) gluons to the free energy at order g 3 , F ( 3 : c . g . ) , is singular: F ( 3 : c . g . ) ≠ 0 without holonomy, but F ( 3 : c . g . ) = 0 when the holonomy is nonzero, even infinitesimally. We show that the absence of the charged gluon contribution is required by gauge invariance alone and is therefore a universal feature.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Modelling Temporal Variability in the Carbon Balance of a Spruce/Moss Boreal Forest

A model of the daily carbon balance of a black spruce/feathermoss boreal forest ecosystem was developed and results compared to preliminary data from the 1994 BOREAS field campaign in northern Manitoba, Canada. The model, driven by daily weather conditions, simulated daily soil climate status (temperature and moisture profiles), spruce photosynthesis and respiration, moss photosynthesis and respiration, and litter decomposition. Model agreement with preliminary field data was good for net ecosystem exchange (NEE), capturing both the asymmetrical seasonality and short-term variability. During the growing season simulated daily NEE ranged from -4 g C m(exp -2) d(exp -1) (carbon uptake by ecosystem) to + 2 g C m(exp -2) d(exp -1) (carbon flux to atmosphere), with fluctuations from day to day. In the early winter simulated NEE values were + 0.5 g C m(exp -2) d(exp -1), dropping to + 0.2 g C m(exp -2) d(exp -1) in mid-winter. Simulated soil respiration during the growing season (+ 1 to + 5 g C m(exp -2) d(exp -1)) was dominated by metabolic respiration of the live moss, with litter decomposition usually contributing less than 30% and live spruce root respiration less than 10% of the total. Both spruce and moss net primary productivity (NPP) rates were higher in early summer than late summer. Simulated annual NEE for 1994 was -51 g C m(exp -2) y(exp -1), with 83% going into tree growth and 17% into the soil carbon accumulation. Moss NPP (58 g C m(exp -2) d(exp -1)) was considered to be litter (i.e. soil carbon input; no net increase in live moss biomass). Ecosystem respiration during the snow-covered season (84 g Cm(exp -2)) was 58% of the growing season net carbon uptake. A simulation of the same site for 1968-1989 showed about 10-20% year-to-year variability in heterotrophic respiration (mean of + 113 g C m-2 y@1). Moss NPP ranged from 19 to 114 g C m(exp -2) y(exp -1); spruce NPP from 81 to 150 g C nt-2 y,@l; spruce growth (NPP minus litterfall) from 34 to 103 g C m(exp -2) y(exp -1); NEE ranged from +37 to -142 g C m(exp -2) y(exp -1). Values for these carbon balance terms in 1994 were slightly smaller than the 1969 - 89 means. Higher ecosystem productivity years (more negative NEE) generally had early springs and relatively wet summers; lower productivity years had late springs and relatively dry summers.

Frolking, S.↗

Integrated G and C Implementation within IDOS: A Simulink Based Reusable Launch Vehicle Simulation

The implementation of multiple Integrated Guidance and Control (IG&C) algorithms per flight phase within a vehicle simulation poses a daunting task to coordinate algorithm interactions with the other G&C components and with vehicle subsystems. Currently being developed by Universal Space Lines LLC (USL) under contract from NASA, the Integrated Development and Operations System (IDOS) contains a high fidelity Simulink vehicle simulation, which provides a means to test cutting edge G&C technologies. Combining the modularity of this vehicle simulation and Simulink s built-in primitive blocks provide a quick way to implement algorithms. To add discrete-event functionality to the unfinished IDOS simulation, Vehicle Event Manager (VEM) and Integrated Vehicle Health Monitoring (IVHM) subsystems were created to provide discrete-event and pseudo-health monitoring processing capabilities. Matlab's Stateflow is used to create the IVHM and Event Manager subsystems and to implement a supervisory logic controller referred to as the Auto-commander as part of the IG&C to coordinate the control system adaptation and reconfiguration and to select the control and guidance algorithms for a given flight phase. Manual creation of the Stateflow charts for all of these subsystems is a tedious and time-consuming process. The Stateflow Auto-builder was developed as a Matlab based software tool for the automatic generation of a Stateflow chart from information contained in a database. This paper describes the IG&C, VEM and IVHM implementations in IDOS. In addition, this paper describes the Stateflow Auto-builder.

Fisher, Joseph E.↗

Observations of Stably Stratified Flow through a Microscale Gap

This paper reports the findings of a comprehensive field investigation on flow through a mountain gap subject to a range of stably stratified environmental conditions. This study was embedded within the Perdigão field campaign, which was conducted in a region of parallel double-ridge topography with ridge-normal wind climatology. One of the ridges has a well-defined gap (col) at the top, and an array of in situ and remote sensors, including a novel triple Doppler lidar system, was deployed around it. The experimental design was mostly guided by previous numerical and theoretical studies conducted with an idealized configuration where a flow (with characteristic velocity U 0 and buoyancy frequency N ) approaches normal to a mountain of height h with a gap at its crest, for which the governing parameters are the dimensionless mountain height G= (Nh)/U 0 and various gap aspect ratios. Modified forms of G were proposed to account for real-world atmospheric variability, and the results are presented in terms of a gap-averaged value G c . The nature of gap flow was highly dependent on G c , wherein a nearly neutral flow regime (G c <1), a transitional mountain wave regime (G c ~O(1)), and a gap-jetting regime ( G c >O(1)) were identified. The measurements were in broad agreement with previous numerical and theoretical studies on a single ridge with a gap or double ridge topography, although details vary. This is the first ever detailed field study reported on microscale (O(100 m)) gap flows, and provides useful data and insights for future theoretical and numerical studies.

Microscale, Gap, environmental fluid dynamics, win↗

Targeted synthesis and reaction mechanism discussion of Mo 2 C based insertion-type electrodes for advanced pseudocapacitors

Mo 2 C is one of the few compounds that possess good electronic conductivity. Meanwhile, it possesses a natural 1D zigzag tunnel structure that is ideally suited for fast ion diffusion. Here, an effective approach is demonstrated for fabrication of structurally stable N-doped Mo 2 C/C nanobelts. They demonstrate high and fast energy storage ability with initial capacitances of 1139 C g -1 at 1 mV s-1, 151 C g -1 at an extremely high scan rate of 2000 mV s -1 and 208 C g -1 at a discharge current density of 200 A g -1 . After electrochemical activation of cycling, the capacity is continuously enhanced and much higher capacitances of 2523 C g -1 at 1 mV s -1 and 1403 C g -1 at 50 A g -1 are achieved after 15,000 cycles at 50 mV s -1 . Using the power law, it is evaluated that a surface-controlled capacitive process makes the main contribution to the capacity, which is over 90% when the scan rates are higher than 10 mV s -1 and still high as 73% at 1 mV s -1 . From in situ synchrotron XRD, it is found that there is a negligible change in the crystal structure and volume during charging/discharging, reflecting an insertion-type charge storage mechanism.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Perturbative study of the one-dimensional quantum clock model

In this report we calculate the ground-state energy density ε(g) for the one-dimensional N-state quantum clock model up to order 18, where g is the coupling and N = 3,4,5,...,10, 20. Using methods based on the Padé approximation, we extract the singular structure of ε"(g) or ε(g). They correspond to the specific heat and free energy of the classical two-dimensional (2D) clock model. We find that, for N = 3, 4, there is a single critical point at g c = 1. The heat capacity exponent of the corresponding 2D classical model is α = 0.34 ± 0.01 for N = 3, and α = - 0.01 ± 0.01 for N = 4. For N > 4, there are two exponential singularities related by g c1 = 1/g c2 , and ε(g) behaves as Ae - $\frac{c}{|g_c -g|^a}$ + analytic terms near g c . The exponent σ gradually grows from 0.2 to 0.5 as N increases from 5 to 9, and it stabilizes at 0.5 when N > 9. The phase transitions exhibited in these models should be generalizations of the Kosterlitz-Thouless transition, which has σ = 0.5 .

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Relative Controls of Vapor Pressure Deficit and Soil Water Stress on Canopy Conductance in Global Simulations by an Earth System Model

Climate change impacts drought and vapor pressure deficit (VPD) that critically influence land-atmosphere water vapor and CO 2 exchange. There is a lack of model-based global evaluation of the relative roles of the two types of water stresses in limiting canopy conductance (G c ). Using the Energy Exascale Earth System Model (E3SM), we conducted four land-atmosphere coupled simulations for the historical period with and without a plant hydraulics scheme (PHS) and under present-day and quadrupling of CO 2 concentrations. By fitting an empirical model of G c to the E3SM simulations, we evaluated the sensitivity of G c to water stresses caused by VPD and soil water (SW) in the growing season. The empirical model is based on a multiplicative algorithm that adjusts G c according to VPD, soil water saturation, and their interaction. Our results show positive trends of VPD over land during the historical period, with large differences between simulations with and without PHS. Despite the differences, a majority of land pixels that experience both SW and VPD stresses are in semi-arid and tropical savanna ecosystems. G c is more sensitive to VPD than SW in less than 20% of land pixels where the empirical model provides a good fit for the simulated G c . Ignoring SW and VPD interaction may underestimate the relative role of VPD on G c by as much as 28%. Including PHS accentuates SW stress limitation on G c by simulating wetter soil and higher ET that lower VPD. Elevated CO 2 reduces G c by increasing water use efficiency and ameliorates SW stress on G c .

54 ENVIRONMENTAL SCIENCES↗

Pterodactyl: Trade Study for an Integrated Control System Design of a Mechanically Deployed Entry Vehicle

This paper presents a trade study method used to evaluate and down-select from a set of guidance and control (G&C) system designs for a mechanically deployable entry vehicle (DEV). The Pterodactyl project, funded by NASA's Space Technology Mission Directorate (STMD), was prompted by the challenge to develop an effective G&C system for a vehicle without a backshell, which is the case for DEVs. For the DEV, the project assumed a specific aeroshell geometry pertaining to an Adaptable, Deployable, Entry Placement Technology (ADEPT) vehicle, which was successfully developed by STMD prior to this study. The Pterodactyl project designed three different G&C systems for the vehicle's precise entry, which this paper briefly discusses. This paper details the Figures of Merit (FOMs) and metrics used during the course of the project's G&C system assessment. Each G&C configuration was traded against the three FOMs categories: G&C system performance, affordability and life cycle costs, and safety and mission success. The relative importance of the FOMs was determined from the Analytical Hierarchy Process (AHP), which was used to develop weights that were combined with quantitative design metrics and engineering judgement to rank the G&C systems against one another. This systematic method takes into consideration the project's input while simultaneously reducing unintentional judgement bias and ultimately was used to select a single G&C design for the project to continue pursuing in the next prototyping and testing phase.

Alunni, Antonella I.↗

Ponderosa pine hydraulic stress predicts more extreme wildfire behavior under future conditions in Bandelier National Monument, New Mexico

Background Live fuel moisture contributes to wildfire spread and reflects plant stress and physiological traits. The anticipated change in live fuel moisture under future conditions is likely non-linear, owing to physiological plant thresholds in water hydraulics. We constructed a mechanistic model of live fuel moisture’s response to water stress to understand the impact of future climate on live fuel moisture. We first gathered data on plant physiology and live fuel moisture for Pinus ponderosa at Bandelier National Monument, NM, USA, and modeled their relationship. We then parameterized a mechanistic plant hydrodynamics model (FATES-HYDRO) to simulate changes in plant stress and a statistical model to simulate the resulting impact on live fuel moisture. We then simulated FATES-HYDRO under future climate anomalies (SSP2-4.5 and SSP5-8.5: 2080–2100) to understand the change in plant stress and estimate its impact on live fuel moisture. Results We found that the number of days below crucial thresholds of live fuel moisture (100% and 79%) increased from contemporary levels (< 100%: 72 days, < 79%: 1.4 days) under SSP2-4.5 (< 100%: 185 days, < 79%: 10.2 day) and increased exponentially under SSP5-8.5 (< 100%: 215 days, < 79%: 65 days). We found that gross primary productivity decreased under both future climate scenarios (contemporary: 336 g C m −2 , SSP2-4.5: 203 g C m −2 , SSP5-8.5: 243 g C m −2 ); however, spring productivity increased under SSP5-8.5, potentially altering fuel loading. We additionally see a potentially lethal loss of conductivity in hydraulic P. ponderosa under SSP5-8.5. Conclusions The overall increase in plant water stress (as represented by loss of hydraulic conductivity and leaf water potential) and lower live fuel moisture appear to be driven by reduced precipitation during late summer monsoons typical of the region, extending the fire season later in the year. We predict increasing variability in the P. ponderosa wildfire regime under both climate projections driven by changing productivity, rising mortality, and an overall decrease in live fuel moisture.

54 ENVIRONMENTAL SCIENCES↗

Carbon budget of the Harvard Forest Long-Term Ecological Research site: pattern, process, and response to global change

How, where, and why carbon (C) moves into and out of an ecosystem through time are long-standing questions in biogeochemistry. Here, we bring together hundreds of thousands of C-cycle observations at the Harvard Forest in central Massachusetts, USA, a mid-latitude landscape dominated by 80–120-yr-old closed-canopy forests. These data answered four questions: (1) where and how much C is presently stored in dominant forest types; (2) what are current rates of C accrual and loss; (3) what biotic and abiotic factors contribute to variability in these rates; and (4) how has climate change affected the forest’s C cycle? Harvard Forest is an active C sink resulting from forest regrowth following land abandonment. Soil and tree biomass comprise nearly equal portions of existing C stocks. Net primary production (NPP) averaged 680–750 g C·m –2 ·yr –1 ; belowground NPP contributed 38–47% of the total, but with large uncertainty. Mineral soil C measured in the same inventory plots in 1992 and 2013 was too heterogeneous to detect change in soil-C pools; however, radiocarbon data suggest a small but persistent sink of 10–30 g C·m –2 ·yr –1 . Net ecosystem production (NEP) in hardwood stands averaged ~300 g C·m –2 ·yr –1 . NEP in hemlock-dominated forests averaged ~450 g C·m –2 ·yr –1 until infestation by the hemlock woolly adelgid turned these stands into a net C source. Since 2000, NPP has increased by 26%. For the period 1992–2015, NEP increased 93%. The increase in mean annual temperature and growing season length alone accounted for ~30% of the increase in productivity. Interannual variations in GPP and NEP were also correlated with increases in red oak biomass, forest leaf area, and canopy-scale light-use efficiency. Compared to long-term global change experiments at the Harvard Forest, the C sink in regrowing biomass equaled or exceeded C cycle modifications imposed by soil warming, N saturation, and hemlock removal. Furthermore, the results of this synthesis and comparison to simulation models suggest that forests across the region are likely to accrue C for decades to come but may be disrupted if the frequency or severity of biotic and abiotic disturbances increases.

54 ENVIRONMENTAL SCIENCES↗

Small Body GN&C Research Report: A G-Guidance Enhancement to Increase Mission Functionality

G-Guidance is a robust G&C (guidance and control) algorithm developed under the small-body GN&C task. The G-Guidance scheme utilizes a model predictive control approach, along with a convexification of the governing dynamics, control constraints, and trajectory/state constraints. This report details an enhancement to the FY2005 G-Guidance algorithm; the addition of a fire-second approach helps to nullify velocity errors and hit desired final velocities much more precisely than the original Fire-first scheme developed in FY2005. The enhancement preserves the primary benefit of the algorithm, which is to ensure required thruster silent times during trajectory maneuvers. The fire-second scheme increases the versatility of the algorithm for missions employing G-Guidance. For instance, a landing sequence could employ the fire-second scheme to ensure a null final velocity at landing, whereas an ascent sequence could utilize a fire-first scheme to ensure there is no drift of the spacecraft toward the ground at the beginning of the maneuver. Examples are provided within to demonstrate a fire-first versus fire-second guidance scheme. As in the existing G-Guidance algorithm, the examples and schemes incorporate gravity models and thruster firing times into discrete dynamics that are solved as a optimal control problem to minimize fuel consumption or thruster energy expenditure.

velocity errors↗