SEARCH · Search NASA

Results for “particle-in-cell”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

2,856 records · Page 69

All sources NASA NTRS DOE OSTI

NASA NTRS2026

Thermochemically-Closed Sonic-Flow Inversion for Enthalpy and Temperature in Multispecies Arc-Jet Flows

A thermochemically-closed sonic-flow inversion framework (TSIF) is developed to infer bulk enthalpy and total temperature upstream of a choked nozzle in arc-jet flows. The formulation recasts a pressure-rise total enthalpy quantification technique as an inverse problem in characteristic-velocity c * space using measured mass flow rate, upstream total pressure, gas composition, and nozzle throat geometry as inputs. Unlike calorimetric energy-balance approaches or optical diagnostics, the method relies primarily on routinely measured facility quantities combined with explicit thermochemical closure. Thermochemical states are obtained using NASA’s open-source Chemical Equilibrium with Applications (CEA) code, enabling construction of a chemistry-consistent relation between characteristic velocity, total enthalpy, and total temperature under equilibrium or frozen assumptions. A discharge coefficient is self-calibrated using cold-flow (arc-off) operation data and applied to hot-flow (arc-on) measurements, enabling upstream losses to be accounted for without empirical correlations. The framework is applied to air, N 2 , and CO 2 –N 2 arc-jet flows and demonstrates expected trends for the inferred thermochemical states as function of arc power, specific energy input, mass-flow, heater configuration, and test gas. In the air limit, under equilibrium assumptions, the method recovers the classical high-enthalpy asymptotic correlation of Winovich with a mean residual of 4.4%, demonstrating compatibility with established sonic-flow scaling, while extending applicability to arbitrary multi-species mixtures and non-equilibrium chemistry. The framework provides a mixture-flexible methodology for determining bulk thermochemical states in modern arc-jet environments using routine facility pressure, mass-flow, gas-composition, and nozzle-geometry information together with a cold-flow calibration.

inviscid theory↗

DOE OSTI2025

Magnetically and optically active edges in phosphorene nanoribbons

Abstract Nanoribbons, nanometre-wide strips of a two-dimensional material, are a unique system in condensed matter. They combine the exotic electronic structures of low-dimensional materials with an enhanced number of exposed edges, where phenomena including ultralong spin coherence times 1,2 , quantum confinement 3 and topologically protected states 4,5 can emerge. An exciting prospect for this material concept is the potential for both a tunable semiconducting electronic structure and magnetism along the nanoribbon edge, a key property for spin-based electronics such as (low-energy) non-volatile transistors 6 . Here we report the magnetic and semiconducting properties of phosphorene nanoribbons (PNRs). We demonstrate that at room temperature, films of PNRs show macroscopic magnetic properties arising from their edge, with internal fields of roughly 240 to 850 mT. In solution, a giant magnetic anisotropy enables the alignment of PNRs at sub-1-T fields. By leveraging this alignment effect, we discover that on photoexcitation, energy is rapidly funnelled to a state that is localized to the magnetic edge and coupled to a symmetry-forbidden edge phonon mode. Our results establish PNRs as a fascinating system for studying the interplay between magnetism and semiconducting ground states at room temperature and provide a stepping-stone towards using low-dimensional nanomaterials in quantum electronics.

Science & Technology - Other Topics↗

NASA NTRS2026

Spatially Refined Satellite Gravimetry Captures Human Signatures in Global Terrestrial Water Storage Trends

Human activities have directly altered the water cycle through water management, aquifer pumping, agricultural irrigation, and land use change. Although satellite gravimetry has transformed global hydrological research, its coarse resolution limits attribution of freshwater change to human activities at many management-relevant scales. Here we assessed global terrestrial water storage (TWS) trends from April 2002 to November 2025 using “stacked” regression of Level-1B intersatellite ranging data, which leverages temporal information and variability to dramatically improve effective spatial resolution relative to standard approaches. We combined this refined product with rigorous uncertainty analysis, autocorrelation-robust geostatistical methods, and literature assessment to evaluate TWS trend associations with land and water use, climate variability, and glacial mass loss. We identified TWS trend hotspots exhibiting significant spatial associations with anthropogenic 40 drivers, including groundwater and surface-water irrigation, rainfed agriculture, deforestation, and reservoir impoundment. Across these regions, cumulative TWS losses (3,122 Gt) substantially exceeded gains (2,432 Gt). Compared with traditional regression of monthly mascons, our approach yielded regional trend magnitudes that are on average 33% larger, revealing that global freshwater depletion, particularly from groundwater pumping, is considerably more acute than previously estimated. Multivariate regression models show that humans account for a significant share of the spatial variability in TWS trends on every non-polar continent except Australia. We detected localized TWS gains linked to rainfed agriculture, surface water irrigation, and reservoir filling that were unresolved in earlier gravimetric studies. The methodology provides a foundation for future gravity missions to independently track decadal freshwater change with unprecedented spatial fidelity.

groundwater↗

DOE OSTI2026

Transition from Vehicular to Structural Ionic Transport in Electrified Alkali Aqueous Solutions

A molecular understanding of the solvation and dynamics of ions under static electric fields is crucial for modelling a wide range of natural and technological processes. Yet, traditional simulation methods suffer from a trade-off that has to be made between accuracy and statistical convergence. To bridge this gap, herein we extend our recently introduced Perturbed Neural Network Potential Molecular Dynamics (PNNP MD) approach to investigate the solvation structures and ionic transport mechanisms of electrified alkali cationic solutions. We obtain ionic conductivities for Li+, Na+ and Cs+ from the field dependence of the ionic current density in good agreement with experiment. Surprisingly, the migration mechanism is found to be strikingly different for the three ions despite their similar ionic conductivities. While Li+ conducts predominantly through vehicular migration of a stable 4-fold coordinated ion at all field strengths, Cs+ conducts strictly through a structural diffusion mechanism, where 9-12 transient first shell water coordination bonds are broken and renewed. Notably, aqueous Na+ emerges as a “Goldilocks” ion: its ion-water interactions are strong enough to maintain distinct 5-6 fold coordination shells at zero field (unlike Cs+) yet labile enough to be strongly perturbed by electric fields (unlike Li+). As a consequence, we observe an electric field-induced transition from vehicular to structural ionic transport for Na+ that is accompanied by a marked increase in the ionic current density. Our results imply that the conductance mechanism of ions with moderate ion-solvent interactions can be effectively tuned by external electric fields.

Joll, Kit↗

DOE OSTI2025

Internally Catalyzed Hydrogen Atom Transfer (I-CHAT)—A New Class of Reactions in Combustion Chemistry

The current paradigm of low-T combustion and autoignition of hydrocarbons is based on the sequential two-step oxygenation of fuel radicals. The key chain-branching occurs when the second oxygenation adduct (OOQOOH) is isomerized releasing an OH radical and a key ketohydroperoxide (KHP) intermediate. The subsequent homolytic dissociation of relatively weak O–O bonds in KHP generates two more radicals in the oxidation chain leading to ignition. Based on the recently introduced intramolecular “catalytic hydrogen atom transfer” mechanism (J. Phys. Chem. 2024, 128, 2169), abbreviated here as I-CHAT, we have identified a novel unimolecular decomposition channel for KHPs to form their classical isomers—enol hydroperoxides (EHP). The uncertainty in the contribution of enols is typically due to the high computed barriers for conventional (“direct”) keto–enol tautomerization. Remarkably, the I-CHAT dramatically reduces such barriers. The novel mechanism can be regarded as an intramolecular version of the intermolecular relay transfer of H-atoms mediated by an external molecule following the general classification of such processes (Catal. Rev.-Sci. Eng. 2014, 56, 403). Here, we present a detailed mechanistic and kinetic analysis of the I-CHAT-facilitated pathways applied to n-hexane, n-heptane, and n-pentane models as prototype molecules for gasoline, diesel, and hybrid rocket fuels. We particularly examined the formation kinetics and subsequent dissociation of the γ-enol-hydroperoxide isomer of the most abundant pentane-derived isomer γ-C5-KHP observed experimentally. To gain molecular-level insight into the I-CHAT catalysis, we have also explored the role of the internal catalyst moieties using truncated models. All applied models demonstrated a significant reduction in the isomerization barriers, primarily due to the decreased ring strain in transition states. In addition, the longer-range and sequential H-migration processes were also identified and illustrated via a combined double keto–enol conversion of heptane-2,6-diketo-4-hydroperoxide as a potential chain-branching model. To assess the possible impact of the I-CHAT channels on global fuel combustion characteristics, we performed a detailed kinetic analysis of the isomerization and decomposition of γ-C5-KHP comparing I-CHAT with key alternative reactions—direct dissociation and Korcek channels. Calculated rate parameters were implemented into a modified version of the n-pentane kinetic model developed earlier using RMG automated model generation tools (ACS Omega, 2023, 8, 4908). Simulations of ignition delay times revealed the significant effect of the new pathways, suggesting an important role of the I-CHAT pathways in the low-T combustion of large alkanes.

Biochemistry & Molecular Biology↗

DOE OSTI2026

Emergence of low-energy spin waves in superconducting electron-doped cuprates

In order to fully utilize the technological potential of unconventional superconductors, an enhanced understanding of the superconducting mechanism is necessary. In the best performing superconductors, the cuprates, superconductivity is intimately linked with magnetism, although the details of this coupling remain elusive. Here, we address this gap by studying the electron-doped cuprate Nd 1.85 Ce 0.15 CuO 4−δ that has an antiferromagnetic ground state when synthesized and only becomes superconducting after a reductive annealing process. Using neutron spectroscopy, we show that the as-grown crystal exhibits a large spin pseudogap in the magnetic fluctuation spectrum. Annealing removes defects introduced by the commonly employed synthesis method and significantly reduces the spin pseudogap. While the spin pseudogap in the annealed sample likely arises from superconductivity, in the as-grown sample it results from the absence of long-wavelength spin waves. These results reveal a direct connection between defects, magnetism, and superconductivity, offering new insight into the mechanisms underlying high-temperature superconductivity and guiding the design of improved superconducting materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS2003

Columbia Accident Investigation Board Report: Volume Four

Volume IV of the Report contains appendices that were not cited in Volume I. These consist of documents produced by NASA and other organizations, which were provided to the Columbia Accident Investigation Board in support of its inquiry into the February 1, 2003 destruction of the Space Shuttle Columbia. The documents are compiled in this volume in the interest of establishing a complete record, but they do not necessarily represent the views of the Board. Volume I contains the Board's findings, analysis, and recommendations. The documents in Volume III through V are also contained in their original color format on the DVD disc in the back of Volume II. Volumes II – VI contain a number of conclusions and recommendations, several of which were adopted by the Board in Volume I. The other conclusions and recommendations drawn in Volumes II – VI do not necessarily reflect the opinion of the Board, but are included for the record. When there is conflict, Volume I takes precedence.

NASA↗

DOE OSTI2025

Tunable high Néel temperature and large anomalous Hall response in antiferromagnetic Weyl semimetal Mn 3 Sn 1− x Ga x thin films

Antiferromagnetic Weyl semimetals based on Mn 3 X(X = Ge, Sn, Ga) kagome compounds exhibit the same ferromagnetic-like responses, including anomalous Hall, Nernst, and magneto-optical effects, as recently discussed for altermagnets. Driven by the Berry curvature due to Weyl fermions, these materials show a disproportionately large magnitude of electromagnetic effects even in the absence of large magnetization. For applications it is crucial to realize these responses in a wide range of temperatures both below and above 300 K. While stoichiometric Mn 3 X materials do not offer optimal performance, we show that Mn 3 Sn 1−x Ga x sputtered films with a variable composition offers a tunable Néel temperature, T N ≈ 425 ± 6–500 ± 15 K, which is crucial for device applications, together with a large tunable anomalous Hall effect. Our thin film growth method enables continuous and precise control over the film composition between x = 0 and x = 1. Through a detailed magnetization and Hall transport, we establish the magnetic phase diagram for the hexagonal Mn 3 Sn 1−x Ga x . Our results reveal an enhanced T N and antichiral magnetic phase in Ga-doped Mn 3 Sn and an enhanced anomalous Hall magnitude in Sn-doped Mn 3 Ga compared to their stoichiometric undoped forms. Our work demonstrates a route to optimize the technologically relevant antiferromagnets for various applications.

Magnetic properties and materials↗

DOE OSTI2026

Mitigating electrochemical degradation in CsPbBr{sub 3} gamma detectors by organic and inorganic encapsulation.

CsPbBr3 perovskite semiconductors have emerged as a leading candidate for nextgeneration radiation detectors because of their exceptional charge transport properties, defect tolerance, and record-breaking sensitivity and energy resolution. Their long-term stability, however, is hindered by electrode-driven electrochemical decomposition, which is accelerated by moisture- and oxygen-assisted ion migration during operation. Here, we investigated organic and inorganic encapsulation strategies as both environmental barriers and means to suppress interfacial degradation pathways. Atomic layer deposition (ALD) of Al2O3 provided a conformal passivation layer that blocked environmental ingress, suppressed ionic diffusion, reduced leakage current, enhanced energy resolution and expanded the operational electric-field window beyond 5 kV∙cm1 . By contrast, organic encapsulants such as paraffin wax and polystyrene slowed moisture diffusion but did not suppress interfacial reactions, with wax extending stability to over 90 days. These results show that ALD-Al2O3 suppresses dominant interfacial degradation pathways, enabling stable, high-field operation and advancing the practical deployment of CsPbBr3 γ-ray detectors.

Unal, Mustafa↗

DOE OSTI2025

Direct observation of distinct bulk and edge nonequilibrium spin accumulation in ultrathin MoTe 2

Low-symmetry two-dimensional (2D) topological materials such as MoTe 2 host efficient charge-to-spin conversion (CSC) mechanisms that can be harnessed for novel electronic and spintronic devices. However, the nature of the various CSC mechanisms and their correlation with underlying crystal symmetries remain unsettled. In this work, we use local spin-sensitive electrochemical potential measurements to directly probe the spatially dependent nonequilibrium spin accumulation in MoTe 2 flakes down to four atomic layers. We are able to clearly disentangle contributions originating from the spin Hall and Rashba-Edelstein effects and uncover an abundance of unconventional spin polarizations that develop uniquely in the sample bulk and edges with decreasing thickness. Using ab-initio calculations, we construct a unified understanding of all the observed CSC components in relation to the material dimensionality and stacking arrangement. Our findings not only illuminate previous CSC results on MoTe 2 but also have important ramifications for future devices that can exploit the local and layer-dependent spin properties of this 2D topological material.

Science & Technology - Other Topics↗

DOE OSTI2025

Synthesis‐Related Nanoscale Defects in Mo‐Based Janus Monolayers Revealed by Cross‐Correlated AFM and TERS Imaging

2D Janus transition metal dichalcogenides (TMDs) are promising candidates for various applications including non-linear optics, energy harvesting, and catalysis. These materials are usually synthesized via chemical conversion of pristine TMDs. Nanometer-scale characterization of the obtained Janus materials’ morphology and local composition is crucial for both the synthesis optimization and the future device applications. In this work, we present the results of cross-correlated atomic force microscopy (AFM) and tip-enhanced Raman spectroscopy (TERS) study of Janus monolayers synthesized by the hydrogen plasma-assisted chemical conversion of MoSe 2 and MoS 2 . We demonstrate that the choice of both the growth substrate and the starting TMD influences the residual strain, thereby shaping the nanoscale morphology of the resulting Janus material. Furthermore, by employing TERS imaging, we show the presence of nanoscale islands (≈20 nm across) of MoSe 2 - ${\mathrm{Mo}}_{{\mathrm{Se}}}^{\mathrm{S}}$ (MoS 2 -${\mathrm{Mo}}_{\mathrm{S}}^{{\mathrm{Se}}}$) vertical heterostructures originating from the bilayer nanoislands in the precursor monolayer crystals. The understanding of the origins of nanoscale defects in Janus TMDs revealed in this study can help with further optimization of the Janus conversion process towards uniform and wrinkle-/crack-free Janus materials. Moreover, this work shows that cross-correlated AFM and TERS imaging is a powerful and accessible method for studying nanoscale composition and defects in Janus TMD monolayers.

2D materials↗

DOE OSTI2026

Electrically Accelerated Mechanochemical Film Formation by a Phosphonium Phosphate Ionic Liquid: An In Situ Chemical Kinetics Investigation

Powertrains in electric vehicles are exposed to stray currents that accelerate wear and cause failure of mechanical components. These durability issues are further aggravated when using low-viscosity lubricants, which are desired for energy efficiency but create harsher contact conditions at sliding interfaces. This study investigates a phosphonium phosphate ionic liquid as a performance-enhancing additive in a low-viscosity base oil for lubricating electrified sliding interfaces. Ionic liquids can adsorb and react on contact interfaces via stress-assisted chemical reactions, generating nanometric tribofilms that provide protection against wear. However, the effect of electric fields on the mechanochemistry of ionic liquids is poorly understood, hindering their adoption in lubricants for electrified powertrains. This article reports an in situ optical interferometry study of ionic liquid derived tribofilm growth kinetics at stressed sliding/rolling interfaces under direct currents. The application of electric currents accelerated tribofilm formation up to a critical current density (∼1.4 A/mm 2 ), beyond which pitting-induced wear dominated. The tribofilms were composed of iron phosphates, iron oxides, and carbon species, with iron oxides becoming predominant under applied currents. These tribofilms prevented the scuffing failure of steel surfaces under electrified conditions. Based on the results, a kinetic model is proposed that integrates electric current effect into the classical stress-assisted thermal activation framework to allow prediction of tribofilm growth at electrified sliding contacts. This framework provides crucial guidance for designing next-generation lubricants for electrified transportation and power generation systems.

additives↗

DOE OSTI2025

Tree tensor network hierarchical equations of motion based on time-dependent variational principle for efficient open quantum dynamics in structured thermal environments

In this work, we introduce an efficient method, TTN-HEOM, for exactly calculating the open quantum dynamics for driven quantum systems interacting with highly structured bosonic baths by combining the tree tensor network (TTN) decomposition scheme with the bexcitonic generalization of the numerically exact hierarchical equations of motion (HEOM). The method yields a series of quantum master equations for all core tensors in the TTN that efficiently and accurately capture the open quantum dynamics for non-Markovian environments to all orders in the system–bath interaction. These master equations are constructed based on the time-dependent Dirac–Frenkel variational principle, which isolates the optimal dynamics for the core tensors given the TTN ansatz. The dynamics converges to the HEOM when increasing the rank of the core tensors, a limit in which the TTN ansatz becomes exact. We introduce TENSO, tensor equations for non-Markovian structured open systems, as a general-purpose Python code to propagate the TTN-HEOM dynamics. We implement three general propagators for the coupled master equations: two fixed-rank methods that require a constant memory footprint during the dynamics and one adaptive-rank method with a variable memory footprint controlled by the target level of computational error. We exemplify the utility of these methods by simulating a two-level system coupled to a structured bath containing one Drude–Lorentz component and eight Brownian oscillators, which is beyond what can presently be computed using the standard HEOM. Our results show that the TTN-HEOM is capable of simulating both dephasing and relaxation dynamics of driven quantum systems interacting with structured baths, even those of chemical complexity, with an affordable computational cost.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Anticipating Optical Availability in Hybrid RF/FSO Links Using RF Beacons and Deep Learning

Radiofrequency (RF) communications offer reliable but low data rates and energy-inefficient satellite links, while free-space optical (FSO) promises high bandwidth but struggles with disturbances imposed by atmospheric effects. A hybrid RF/FSO architecture aims to achieve optimal reliability along with high data rates for space communications. Accurate prediction of dynamic ground-to-satellite FSO link availability is critical for routing decisions in low-earth orbit constellations. In this paper, we propose a system leveraging ubiquitous RF links to proactively forecast FSO link degradation prior to signal drops below threshold levels. This enables pre-calculation of rerouting to maximally maintain high data rate FSO links throughout the duration of weather effects. We implement a supervised learning model to anticipate FSO attenuation based on the analysis of RF patterns. Through the simulation of a dense lower earth orbit (LEO) satellite constellation, we demonstrate the efficacy of our approach in a simulated satellite network, highlighting the balance between predictive accuracy and prediction duration. An emulated cloud attenuation model is proposed to provide insight into the temporal profiles of RF signals and their correlation to FSO channel dynamics. Our investigation sheds light on the trade-offs between prediction horizon and accuracy arising from RF beacon numbers and proximity.

FSO availability↗

DOE OSTI2026

Uncovering the linear boron environment in Na 3 BP 2 through solid-state 11 B NMR spectroscopy

Boron-based compounds exhibit a wide range of structural diversity, with potential applications spanning organic and inorganic chemistry. Herein, we focus on the characterization of the linear boron-phosphorus unit P═B═P in Na 3 BP 2 using solid-state nuclear magnetic resonance (ssNMR) spectroscopy and density functional theory (DFT) calculations. High-resolution 11 B ssNMR spectra were recorded at two fields, and key parameters such as chemical shift anisotropy (CSA), quadrupolar coupling constants (C Q ), and electric field gradient (EFG) tensors were extracted. The 11 B NMR results revealed a distinct chemical environment for the two-coordinate boron atom, with a CSA span (Ω) of 280 ppm and a C Q of 3.0 MHz. These values were further validated through periodic plane-wave DFT calculations, which showed good agreement with experimental results. The obtained spectral parameters are compared to other linear boron units, such as the BO 2 motif, providing a broader context for understanding boron coordination in inorganic compounds. This work expands the body of NMR knowledge on boron-containing materials, particularly for linear boron motifs. The findings contribute to the growing field of boron chemistry and its potential applications in advanced materials.

Porter, Andrew P. [Ames Laboratory (AMES), Ames, I↗

DOE OSTI2026

Small Co-doping induced magnetic and electrical transitions in single crystal CaF⁢e 0.95 ⁢C⁢o 0.05 ⁢O 3

CaFe⁢O 3 and CaF⁢e 0.95 ⁢C⁢o 0.05 ⁢O 3 single crystals were grown by combing floating-zone method with high oxygen pressure treatment. Both crystals show charge disproportionation of Fe (F⁡e 4+ → F⁡e 3+ + F⁡e 5+ ) accompanied by a metal-to-insulator transition as well as a crystal structural transformation from Pbnm to 𝑃⁢2 1 /𝑛. However, the slight introduction of Co significantly suppresses the critical temperature of charge disproportionation from 290 K in CaFe⁢O 3 to 260 K in CaF⁢e 0.95 ⁢C⁢o 0.05 ⁢O 3 . Different from the single antiferromagnetic phase transition as observed in the polycrystalline CaFe⁢O 3 , two sequential antiferromagnetic transitions are found to occur in these two single crystals with the Néel temperatures around 119 and 112 K for CaFe⁢O 3 crystal and 109 and 98 K for CaF⁢e 0.95 ⁢C⁢o 0.05 ⁢O 3 crystal, respectively. Neutron diffraction illustrates the formation of a spiral antiferromagnetic structure. Moreover, as the temperature further decreases to 65 K, a third magnetic transition accompanied by a second electrical transition is observed in the slightly Co-doped CaF⁢e 0.95 ⁢C⁢o 0.05 ⁢O 3 crystal. In conclusion, first-principles calculations suggest that the introduction of a moderate amount of Co and the peculiar spiral spin texture play an important role in the presence of such new magnetic and electrical transitions.

Xia, Hailiang [Chinese Academy of Sciences (CAS), ↗

NASA NTRS1993

Prediction of prompt NO(x) in hydrocarbon air flames

The gas turbine industry is directing particular attention to very low NOx combustors, whether for aircraft or land based CCGT systems. These low NOx combustors frequently use liquid fuels or natural gas burning under very lean premixed conditions with air or under rich-lean conditions, although only the first case is studied here. In land based systems, diluted steam or nitrogen are sometimes injected into the combustion chambers to reduce flame temperature. The NOx emissions from such systems are the product of three chemical mechanisms which are interrelated: the hydrocarbon prompt NO, the thermal NO (extended Zeldovich mechanism), and the nitrous oxide route to NO. Formation of NO2 from NO also occurs, as well as emission of carbon monoxide and unburnt hydrocarbons. When the fuel-oxidant proportion decreases towards leaner conditions, flame temperatures are lowered, resulting in the total NO being reduced and the thermal-NO contribution greatly diminished to the benefit of the remaining two mechanisms of NO formation. While knowledge of the elementary reactions and their chemical kinetics concerning methane and simple hydrocarbons combustion has existed for a number of years, its use for computer modeling is limited to simple flow dynamics configurations. Nevertheless, understanding of such combusting flows under a wide range of experimental conditions allows for analogies or speculations with more complex actual systems. Such understanding can be achieved by means of one dimensional laminar premixed flame modeling, with a full chemical mechanism which incorporates the three routes of NO formation. Complementary to this understanding is the modeling of the actual combustion system using a full description of the fluid dynamics coupled with a reduced chemical scheme, which is then compared against the first model. The objective of this investigation is to evaluate the relative importance of the three mechanisms of NO formation in lean premixed methane-air combustion with increasing pressure using the one dimensional plug flow package PREMIX, and to test the validity of a three dimensional model with a global chemical mechanism against the one dimensional model in the atmospheric pressure case. Methane is chosen because it is the only mechanism which is reasonably well known and is a good guide to the behavior of other hydrocarbons. The mixture ratio chosen is richer than that in lean gas turbines, but the combustion of this mixture with the low preheat gives realistic gas turbine final flame temperatures. Conditions of NO2 formation are also analyzed in the one dimensional model and results extrapolated to the case of gas turbines.

Valerie Dupont↗

NASA NTRS1993

A SIMULINK Environment for Flight Dynamics and Control Analysis - Application to the DHC-2 'Beaver', Part 1. Implementation of a Model Library in SIMULINK

The design of advanced Automatic Aircraft Control Systems (AACS's) can be improved upon considerably if the designer can access all models and tools required for control system design and analysis through a graphical user-interface, from within one software environment. This MSc-thesis presents the first step in the development of such an environment, which is currently being done at the Section for Stability and Control of Delft University of Technology, Faculty of Aerospace Engineering. The environment is implemented within the commercially available software package MATLAB/SIMULINK. The report consists of two parts. Part I gives a detailed description of the AACS design environment. The heart of this environment is formed by the SIMULINK implementation of a nonlinear aircraft model in block-diagram format. The model has been worked out for the old laboratory aircraft of the Faculty, the De Havilland DHC-2 'Beaver', but due to its modular structure, it can easily be adapted for other aircraft. Part I also describes MATLAB programs which can be applied for finding steady-state trimmed-flight conditions and for linearization of the aircraft model, and it shows how the built-in simulation routines of SIMULINK have been used for open-loop analysis of the aircraft dynamics. Apart from the implementation of the models and tools, a thorough treatment of the theoretical backgrounds is presented. Part II of this report presents a part of an autopilot design process for the 'Beaver' aircraft, which clearly demonstrates the power and flexibility of the AACS design environment from part I. Evaluations of all longitudinal and lateral control laws by means of nonlinear simulations are treated in detail. The AACS design environment from part I proved to be a very useful tool for designing the control laws of the 'Beaver' autopilot within a very tight time-schedule. The autopilot design process itself will be used as a guideline for future AACS research at the Faculty of Aerospace Engineering. Flight tests of the 'Beaver' autopilot, done after evaluating the control laws in the SIMULINK package, proved to be quite successful. In the future, the AACS design package will evolve into a standardized, integrated design environment which can be applied to virtually any type of aircraft. The AACS design cycle will be shortened further by developing tools for automatically porting control laws from the MATLAB/SIMULINK environment to a piloted real-time flight simulator and the Flight Control Computers of the aircraft.

Aircraft Modelling↗