Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed optimization”

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 1,603 records · Page 89

Data traffic reduction schemes for Cholesky factorization on asynchronous multiprocessor systems

Communication requirements of Cholesky factorization of dense and sparse symmetric, positive definite matrices are analyzed. The communication requirement is characterized by the data traffic generated on multiprocessor systems with local and shared memory. Lower bound proofs are given to show that when the load is uniformly distributed the data traffic associated with factoring an n x n dense matrix using n to the alpha power (alpha less than or equal 2) processors is omega(n to the 2 + alpha/2 power). For n x n sparse matrices representing a square root of n x square root of n regular grid graph the data traffic is shown to be omega(n to the 1 + alpha/2 power), alpha less than or equal 1. Partitioning schemes that are variations of block assignment scheme are described and it is shown that the data traffic generated by these schemes are asymptotically optimal. The schemes allow efficient use of up to O(n to the 2nd power) processors in the dense case and up to O(n) processors in the sparse case before the total data traffic reaches the maximum value of O(n to the 3rd power) and O(n to the 3/2 power), respectively. It is shown that the block based partitioning schemes allow a better utilization of the data accessed from shared memory and thus reduce the data traffic than those based on column-wise wrap around assignment schemes.

Naik, Vijay K.↗

Evaluation of Spectral Band Adjustment Factors for Cross-Calibration of Visible Imagers

The NASA CERES EBAF dataset provides TOA SW and LW fluxes for long-term monitoring of the Earth’s energy balance and to validate climate models. The EBAF products, based on the Terra, Aqua, and NOAA20 CERES instrument observed radiances, rely on coincident measurements from the onboard MODIS or VIIRS imagers to determine cloud properties used for angular distribution model scene selection, which is required to convert the CERES observed radiances into flux values. Furthermore, CERES utilizes geostationary imager (GEO) hourly fluxes and clouds to infer the regional averaged daily flux in the SYN1deg product. A seamless transition of fluxes and clouds can only occur if the analogous MODIS, VIIRS, and GEO channels are properly inter-calibrated. The analogous band SRFs differ noticeably and require scene dependent spectral band adjustment factors (SBAF) for proper radiometric scaling between them. Given their disparate but overlapping SRFs, the coincident VIIRS I1 and M5 band TOA reflectance measurements provide the optimal opportunity to validate SBAFs over many surface and cloud conditions. The CERES project maintains SCIAMACHY, GOME-2, and Hyperion scene-stratified hyper-spectral reflectance measurements that can be convolved with sensor pair SRFs to compute the corresponding SBAF. This study highlights the SCIAMACHY, GOME-2, and Hyperion based SBAFs for different MODIS, VIIRS, and GEO inter-calibration targets, including all-sky tropical ocean (ATO), Libya-4, Dome-C, and deep convective clouds. To mitigate the observed sensor radiance fluctuations due to atmospheric absorption and aerosol variations, the Earth target observed radiances are correlated with multiple atmospheric parameters, such as precipitable water and ozone. The remaining M5 and I1 mean spectral band radiance difference should be resolved by the SBAF correction. The formulation of hyper-spectral sensor based SBAFs and validation methods can be verified in future utilizing high-accuracy, SI-traceable hyper-spectral measurements from CLARREO Pathfinder.

David R Doelling↗

Xplo-SA9T: Optimization of Sorbent Technology With Existing and New Amine Variants for Carbon Dioxide Removal and Humidity Control

Scrubbing of the spacesuit air stream to control humidity and regulate carbon dioxide (CO 2 ) relies on amine-loaded resins, colloquially known as sorbents. The efficacy of Rapid Cycle Amine technology relies on the propensity of nitrogens toward adsorption of CO 2 while simultaneously offering regenerative properties under vacuum or positive pressure. The commonly utilized sorbent, SA9T, was recreated and evaluated further to increase stability and serve as a benchmark comparison for newly synthesized amine sorbent technologies. Ongoing storage studies with eight different environments will allow the development of applicable storage protocols and provide insight into possible ammonia generation pathways. Attempts to alleviate the known issues with SA9T included the alteration of the primary amine (tetraethylenepentamine), the choice of resin, and the particle size distribution. In-house CO 2 breakthrough experiments were conducted using a specialized apparatus that can simulate spacesuit conditions, offering high-quality triplicate capacity analysis affording direct comparison to known adsorbents in the literature. Thermo-Desorption Gas Chromatography Mass Spectrometry studies on the synthesized materials revealed any reagent or solvent offgassing while Nuclear Magnetic Resonance spectroscopy tracked the stability of the bound amine. Scanning Electron Microscope images of sieved materials provided a more accurate size interpretation facilitating synthetic design and pressure drop elimination. All efficient sorbents were synthesized to a commercial batch scale subsequent to evaluation including Xplo-SA9T.

John R Tidwell↗

Correction of DIAL Stratospheric Ozone Measurements in the Presence of Pinatubo Aerosols

NASA Langley's airborne lidar system measured aerosol and ozone distributions in the stratosphere from Jan. - Mar. 1992 as part of the Airborne Arctic Stratospheric expedition (AASE-2). The eruption of Mount Pinatubo in Jun. 1991 has increased the aerosol burden of the stratosphere and thereby increased the importance of applying an aerosol correction to the ozone measurements. The correction relies on a Bernoulli solution to derive a backscatter correction to the differential absorption lidar (DIAL) returns at two wavelengths in the ultraviolet spectral region (lambda(sub on) = 301.5 nm, lambda(sub off) = 310.87 nm) as described in earlier works. This paper discusses how the parameters for the correction were optimized for application to the AASE-2 data set.

Fenn, Marta A.↗

Modal characterization of the ASCIE segmented optics testbed: New algorithms and experimental results

New frequency response measurement procedures, on-line modal tuning techniques, and off-line modal identification algorithms are developed and applied to the modal identification of the Advanced Structures/Controls Integrated Experiment (ASCIE), a generic segmented optics telescope test-bed representative of future complex space structures. The frequency response measurement procedure uses all the actuators simultaneously to excite the structure and all the sensors to measure the structural response so that all the transfer functions are measured simultaneously. Structural responses to sinusoidal excitations are measured and analyzed to calculate spectral responses. The spectral responses in turn are analyzed as the spectral data become available and, which is new, the results are used to maintain high quality measurements. Data acquisition, processing, and checking procedures are fully automated. As the acquisition of the frequency response progresses, an on-line algorithm keeps track of the actuator force distribution that maximizes the structural response to automatically tune to a structural mode when approaching a resonant frequency. This tuning is insensitive to delays, ill-conditioning, and nonproportional damping. Experimental results show that is useful for modal surveys even in high modal density regions. For thorough modeling, a constructive procedure is proposed to identify the dynamics of a complex system from its frequency response with the minimization of a least-squares cost function as a desirable objective. This procedure relies on off-line modal separation algorithms to extract modal information and on least-squares parameter subset optimization to combine the modal results and globally fit the modal parameters to the measured data. The modal separation algorithms resolved modal density of 5 modes/Hz in the ASCIE experiment. They promise to be useful in many challenging applications.

Carrier, Alain C.↗

Technical support for digital systems technology development. Task order 1: ISP contention analysis and control

Alternatives for realizing a packet-based network switch for use on a frequency division multiple access/time division multiplexed (FDMA/TDM) geostationary communication satellite were investigated. Each of the eight downlink beams supports eight directed dwells. The design needed to accommodate multicast packets with very low probability of loss due to contention. Three switch architectures were designed and analyzed. An output-queued, shared bus system yielded a functionally simple system, utilizing a first-in, first-out (FIFO) memory per downlink dwell, but at the expense of a large total memory requirement. A shared memory architecture offered the most efficiency in memory requirements, requiring about half the memory of the shared bus design. The processing requirement for the shared-memory system adds system complexity that may offset the benefits of the smaller memory. An alternative design using a shared memory buffer per downlink beam decreases circuit complexity through a distributed design, and requires at most 1000 packets of memory more than the completely shared memory design. Modifications to the basic packet switch designs were proposed to accommodate circuit-switched traffic, which must be served on a periodic basis with minimal delay. Methods for dynamically controlling the downlink dwell lengths were developed and analyzed. These methods adapt quickly to changing traffic demands, and do not add significant complexity or cost to the satellite and ground station designs. Methods for reducing the memory requirement by not requiring the satellite to store full packets were also proposed and analyzed. In addition, optimal packet and dwell lengths were computed as functions of memory size for the three switch architectures.

Stehle, Roy H.↗

The AutoBayes Program Synthesis System: System Description

AUTOBAYES is a fully automatic program synthesis system for the statistical data analysis domain. Its input is a concise description of a data analysis problem in the form of a statistical model; its output is optimized and fully documented C/C++ code which can be linked dynamically into the Matlab and Octave environments. AUTOBAYES synthesizes code by a schema-guided deductive process. Schemas (i.e., code templates with associated semantic constraints) are applied to the original problem and recursively to emerging subproblems. AUTOBAYES complements this approach by symbolic computation to derive closed-form solutions whenever possible. In this paper, we concentrate on the interaction between the symbolic computations and the deductive synthesis process. A statistical model specifies for each problem variable (i.e., data or parameter) its properties and dependencies in the form of a probability distribution, A typical data analysis task is to estimate the best possible parameter values from the given observations or measurements. The following example models normal-distributed data but takes prior information (e.g., from previous experiments) on the data's mean value and variance into account.

Fischer, Bernd↗

Interactive Inverse Design Optimization of Fuselage Shape for Low-Boom Supersonic Concepts

This paper introduces a tool called BOSS (Boom Optimization using Smoothest Shape modifications). BOSS utilizes interactive inverse design optimization to develop a fuselage shape that yields a low-boom aircraft configuration. A fundamental reason for developing BOSS is the need to generate feasible low-boom conceptual designs that are appropriate for further refinement using computational fluid dynamics (CFD) based preliminary design methods. BOSS was not developed to provide a numerical solution to the inverse design problem. Instead, BOSS was intended to help designers find the right configuration among an infinite number of possible configurations that are equally good using any numerical figure of merit. BOSS uses the smoothest shape modification strategy for modifying the fuselage radius distribution at 100 or more longitudinal locations to find a smooth fuselage shape that reduces the discrepancies between the design and target equivalent area distributions over any specified range of effective distance. For any given supersonic concept (with wing, fuselage, nacelles, tails, and/or canards), a designer can examine the differences between the design and target equivalent areas, decide which part of the design equivalent area curve needs to be modified, choose a desirable rate for the reduction of the discrepancies over the specified range, and select a parameter for smoothness control of the fuselage shape. BOSS will then generate a fuselage shape based on the designer's inputs in a matter of seconds. Using BOSS, within a few hours, a designer can either generate a realistic fuselage shape that yields a supersonic configuration with a low-boom ground signature or quickly eliminate any configuration that cannot achieve low-boom characteristics with fuselage shaping alone. A conceptual design case study is documented to demonstrate how BOSS can be used to develop a low-boom supersonic concept from a low-drag supersonic concept. The paper also contains a study on how perturbations in the equivalent area distribution affect the ground signature shape and how new target area distributions for low-boom signatures can be constructed using superposition of equivalent area distributions derived from the Seebass-George-Darden (SGD) theory.

Li, Wu↗

Development of Multi-Sensor Global Cloud and Radiance Composites for DSCOVR EPIC Imager with Subpixel Definition

The Deep Space Climate Observatory (DSCOVR) enables analysis of the daytime Earth radiation budget via the onboard Earth Polychromatic Imaging Camera (EPIC) and National Institute of Standards and Technology Advanced Radiometer (NISTAR). EPIC delivers adequate spatial resolution imagery but only in shortwave bands (317-780 nm), while NISTAR measures the top-of-atmosphere (TOA) whole-disk radiance in shortwave and longwave broadband windows. Accurate calculation of albedo and outgoing longwave flux requires a high-resolution scene identification such as the radiance observations and cloud properties retrievals from low earth orbit (LEO, including NASA Terra and Aqua MODIS, Suomi-NPP VIIRS, and NOAA AVHRR) and geosynchronous (GEO, including GOES east and west, METEOSAT, INSAT-3D, MTSAT-2, and Himawari-8) satellite imagers. The cloud properties are derived using the Clouds and the Earth's Radiant Energy System (CERES) mission Cloud Subsystem group algorithms. These properties have to be co-located with EPIC pixels to provide the scene identification and to select anisotropic directional models (ADMs), which are then used to adjust the NISTAR-measured radiance and subsequently obtain the global daytime shortwave and longwave fluxes. This work presents an algorithm for optimal merging of selected radiance and cloud property parameters derived from multiple satellite imagers to obtain seamless global hourly composites at 5-km resolution. Selection of satellite data for each 5-km pixel is based on an aggregated rating that incorporates five parameters: nominal satellite resolution, pixel time relative to the EPIC time, viewing zenith angle, distance from day/night terminator, and probability of sun glint. To provide a smoother transition in the merged output, in regions where candidate pixel data from two satellite sources have comparable aggregated rating, the selection decision is defined by the cumulative function of the normal distribution so that abrupt changes in the visual appearance of the composite data are avoided. Higher spatial accuracy in the composite product is achieved by using the inverse mapping with gradient search during reprojection and bicubic interpolation for pixel resampling.

Khlopenkov, Konstantin V.↗

A numerical procedure for the parametric optimization of three dimensional scramjet nozzles

A numerical procedure permitting the rapid determination of the internal performance of a class of scramjet nozzle configurations is presented. The approach developed is based on the construction of quasi two dimensional simple wave networks, where lateral expansion effects are incorporated by one dimensional approximations. A numerical procedure following this approach has has been developed and results obtained are highly comparable to those obtained employing a characteristic procedure. The numerical program developed permits the parametric variation of cowl length, turning angles on the cowl and vehicle undersurface and lateral expansion and is subject to fixed constraints such as the vehicle length and nozzle exit height. The program requires uniform initial conditions at the burner exit station and yields the location of all predominant wave zones, accounting for lateral expansion effects. In addition, the program yields the detailed pressure distribution on the cowl and vehicle undersurface and calculates the nozzle thrust, lift and pitching moment.

Dash, S.↗

Buckling and postbuckling behavior of square compression-loaded graphite-epoxy plates with circular cutouts

Results are presented for unidirectional (0, 10)(sub s) and (90,10)(sub s) plates, ((0/90)(sub 5)(sub s)) plates, and for aluminum plates. Results are also presented for ((+/- theta)(sub 6)(sub s)) angle-ply plates for values of theta = 30, 45, and 60 degrees. The results indicate that the change in axial stiffness of a plate at buckling is strongly dependent upon cutout size and plate orthotropy. The presence of a cutout gives rise to an internal load distribution that changes, sometimes dramatically, as a function of cutout size coupled with the plate orthotropy. In the buckled state, the role of orthotropy becomes more significant since bending in addition to membrane orthotropy is present. Most of the plates with cutouts exhibited less postbuckling stiffness than the corresponding plate without a cutout, and the postbuckling stiffness decreased with increasing cutout size. However, some of the highly orthotropic plates with cutouts exhibited more postbuckling stiffness than the corresponding plate without a cutout. These results suggest the possibility of tailoring the cutout size and the stacking sequence of a composite plate to optimize postbuckling stiffness. It was found that plates with large radius cutouts do exhibit some postbuckling strength. The results also indicate that a cutout can influence modal interaction in a plate. Specifically, results are presented that show a plate with a relatively small cutout buckling at a higher load than the corresponding plate without a cutout, due to modal interaction. Other results are presented that indicate the presence of nonlinear prebuckling deformations, due to material nonlinearity, in the angle-ply plates with theta = 45 and 60 degrees. The nonlinear prebuckling deformations are more pronounced in the plates with theta = 45 degrees and become even more pronounced as the cutout size increases. Results are also presented that show how load-path eccentricity due to improper machining of the test specimens affects the buckling behavior. Some of the plates with cutouts and eccentricity exhibited a snap-through type of buckling behavior.

Nemeth, Michael P.↗

Ground Based, Millimeter Wave Measurement of Ozone in the Middle Atmosphere

There is a need for highly reliable measurements of stratospheric ozone. Policy makers worldwide concerned with public health rely oil a clear consensus from the scientific community as a basis for ozone-related environmental policy that has a significant impact oil national economies. Tile latest Such consensus was presented in WMO, and used in a 1999 meeting of the parties considering amendments to the Montreal Protocol oil Substances that Deplete the Ozone Layer. Tile scientific community, in turn, needs highly precise and accurate measurements of ozone levels, and small time derivatives of these levels, both in continued - development of its understanding of the physical and chemical processes involved and as clear evidence that these processes are occurring as stated. Over most of the world, changes in ozone levels are small. For example, over the heavily populated northern midlatitudes, the linearized rate of ozone decline is between 0.2% per year and 0.7% per year, depending on altitude. These values are small enough to make measurement requirements technically challenging. Data quality may suffer from imperfections in individual instruments. In one instance, early results from a satellite-borne ozone sensor were later found to be invalid because of calibration drift. Even in the absence of drift, tile absolute calibration of a new sensor may differ slightly from that of its predecessor in service. Most ozone remote sensing instruments operate at ultraviolet or infrared wavelengths where scattering from dust and aerosols must be taken into account; results from these systems may be or are affected following a major volcanic eruption, such as tile one at Mt. Pinatubo in 1991. Given these difficulties, a consensus of measurements from several independent systems is required to insure a reliable understanding of stratospheric ozone levels. Because of the above-described need for highly precise and accurate ozone measurements using several independent techniques, there was interest in developing several techniques which were known but not highly developed in the 1980's into systems capable of being used in ail operational manner to make measurements with the level of quality needed to-detect small trends in ozone levels. A ground-based microwave instrument capable of remotely sensing stratospheric ozone had been designed by tile Principal Investigator of the present project. This instrument was built at tile Millitech Corporation in South Deerfield, Massachusetts before tile present work began. (Funding for design and construction of the instrument came from sources other than the present grant.) Tile instrument measured the spectrum of one of the many emission lines produced by purely rotational transitions of ozone, one at a frequency of 110.8 GHz. The altitude distribution of ozone can. in principle, be retrieved from the details of the pressure-broadened spectrum of the ozone transition. However, the level of contamination of the spectral measurement by instrumentally induced artifacts must be very low in order to retrieve a ozone profile of useful quality from it. The Millitech instrument demonstrated spectral purity at ail adequate level, and there were promising ideas for instrumental improvements and for further development of the technique. The initial objectives of the present project, then, were to develop techniques for calibrating the Millitech instrument, to minimize artifacts in tile spectra it produces, to optimally retrieve ozone profiles from tile spectra, to test tile quality of the microwave profiles by comparing them with profiles obtained using several other, independent techniques over both short and periods of time, and to perform research using the ozone data gathered with the instrument.

Parrish, Alan↗

Local Estimators for Spacecraft Formation Flying

A formation estimation architecture for formation flying builds upon the local information exchange among multiple local estimators. Spacecraft formation flying involves the coordination of states among multiple spacecraft through relative sensing, inter-spacecraft communication, and control. Most existing formation flying estimation algorithms can only be supported via highly centralized, all-to-all, static relative sensing. New algorithms are needed that are scalable, modular, and robust to variations in the topology and link characteristics of the formation exchange network. These distributed algorithms should rely on a local information-exchange network, relaxing the assumptions on existing algorithms. In this research, it was shown that only local observability is required to design a formation estimator and control law. The approach relies on breaking up the overall information-exchange network into sequence of local subnetworks, and invoking an agreement-type filter to reach consensus among local estimators within each local network. State estimates were obtained by a set of local measurements that were passed through a set of communicating Kalman filters to reach an overall state estimation for the formation. An optimization approach was also presented by means of which diffused estimates over the network can be incorporated in the local estimates obtained by each estimator via local measurements. This approach compares favorably with that obtained by a centralized Kalman filter, which requires complete knowledge of the raw measurement available to each estimator.

Fathpour, Nanaz↗

Initial Study of Friction Pull Plug Welding

Pull plug friction welding is a new process being developed to conveniently eliminate defects from welded plate tank structures. The general idea is to drill a hole of precise, optimized dimensions and weld a plug into it, filling the hole perfectly. A conically-shaped plug is rotated at high angular velocity as it is brought into contact with the plate material in the hole. As the plug is pulled into the hole, friction rapidly raises the temperature to the point at which the plate material flows plastically. After a brief heating phase, the plug rotation is terminated. The plug is then pulled upon with a forging force, solidly welding the plug into the hole in the plate. Three aspects of this process were addressed in this study. The transient temperature distribution was analyzed based on slightly idealized boundary conditions for different plug geometries. Variations in hole geometry and ram speed were considered, and a program was created to calculate volumes of displaced material and empty space, as well as many other relevant dimensions. The relation between the axially applied forging force and the actual forging pressure between the plate and plug surfaces was determined for various configurations.

Rich, Brian S.↗

Simulation and Analysis of Technology and Operational Procedures to Reduce the Combined Effects of Emissions and Contrails

The development and evaluation of concepts and technology to support future air traffic management systems require a hierarchy of models ranging from real-time simulations to extensive field evaluations. Air traffic simulation models such as Airspace Concept Evaluation System, Center Tracon Automation System, Future Air traffic management Concept Evaluation Tool and others are used to design air traffic systems balancing the conflicting objectives of maximizing safety, meeting future demands for airports and airspace and increase efficiency of traffic flows in the presence of uncertain weather. The impact of aviation emissions and contrails on climate imposes another constraint on the design of aircraft and aviation operations. The understanding of the complex interaction between physical climate system, the carbon and other greenhouse gas emissions and aviation activity can be improved by the development of integrated assessment models that include emission and climate models together with air traffic simulations. The modeling of aircraft emissions and their interaction with each other to change the concentration levels of different gasses in the atmosphere and the resulting impact of the radiative forcing on the equilibrium of the Earth's atmosphere is complex and requires the use of coupled atmosphere-ocean general circulation models together with three-dimensional models of carbon cycle and chemistry of other non-CO2 greenhouse gases. These models are computationally intensive and unsuitable for studies involving the generation of multiple scenarios. Simple emission and climate models, based on the input-output relations of linear systems, capture the fundamental emission to climate impact behavior by careful selection of key variables and their dynamics. The impact of various greenhouse gases depends on the total concentration, effect per unit change in atmospheric concentration and the spatial distribution of the gas. All these quantities are influenced by the lifetime of the gas. The impact of a greenhouse gas depends on the interval of assessment, which may vary from a few decades to a few centuries. Climate metrics are aimed at providing a common scale to compare different greenhouse gases. If the metrics are to be used as a tool in developing and evaluating aviation operations, they should be transparent and easy to apply. Global Warming Potential and Aggregate Global Temperature Potential are some of the commonly used metrics. This paper integrates a national-level air traffic simulation and optimization capability with simple climate models and carbon cycle models, and climate metrics to assess the impact of aviation on climate. The capability brings together metrics, which are useful in aviation operations together with metrics used in climate studies. The capability can be used to make trade-offs between extra fuel cost and reduction in climate impact. There is considerable uncertainty in our understanding of the radiative forcing associated with emissions and contrails. The parameters in the simulation can be used to evaluate the effect of various uncertainties in emission models and contrails. It can also be used to evaluate the impact of different decision horizons. Alternatively, the optimization results from the simulation can be used as inputs to other tools that monetize global climate impacts like the FAA's Aviation Environmental Portfolio Management Tool for Impacts.

impact of aviation on the environment↗

Accelerated Simulation of Air Pollution Using NVIDIA RAPIDS

Atmospheric chemistry models are a central tool to study and forecast the impact of air pollution on the environment, vegetation, and human health. However, the numerical simulation of chemical kinetics is computationally expensive due to the stiffness of the system of ordinary differential equations that describes atmospheric chemistry. Here we present an alternative approach to the computation of atmospheric chemistry based on machine learning. Our training data set is produced using the NASA Goddard Earth Observing System (GEOS) model with GEOS-Chem chemistry, run on the NASA Center for Climate Simulation (NCCS) Discover supercomputing cluster on 384 Intel Xeon Haswell cores. This model spends more than 50% of total run time on solving atmospheric chemistry. The data set contains as input features the air pollution concentrations before solving the differential equations, together with some key physical parameters such as temperature and sun intensity. As target variables we define the air pollution concentrations after solving the differential equations. Using Dask-cuDF and Dask-XGBoost on the NVIDIA RAPIDS platform on 8 Tesla V100 GPUs, we generate from this training set gradient boosted decision tree models that can reproduce the simulation of chemical kinetics. We do this on the NCCS Advanced Data Analytics Platform (ADAPT) science cloud environment. Our application takes full advantage of recent advances in Dask-XGBoost, such as multi-node and multi-GPU scaling for distributed training with large data sets. The increase in training data size enabled by this is critical to capture the full range of chemical environments encountered across the globe and all annual seasons.The boosted tree models offer good predictability and show many of the features of the full chemistry reference simulation. Further improvements can be achieved through mass balance considerations and by accounting for error correlations. We incorporate the boosted tree models into the GEOS reference model using XGBoost's C API. This enables a seamless integration of the GPU trained models into GEOS-Chem, which is written in Fortran and optimized for use in a massively parallel CPU environment. We show the benefits of this approach and discuss the potential speedup of this machine learning accelerated atmospheric chemistry model.

Keller, Christoph A.↗

Constraining Global Coronal Models with Multiple Independent Observables

Global coronal models seek to produce an accurate physical representation of the Sun's atmosphere that can be used, for example, to drive space-weather models. Assessing their accuracy is a complex task, and there are multiple observational pathways to provide constraints and tune model parameters. Here, we combine several such independent constraints, defining a model-agnostic framework for standardized comparison. We require models to predict the distribution of coronal holes at the photosphere, and neutral line topology at the model's outer boundary. We compare these predictions to extreme-ultraviolet (EUV) observations of coronal hole locations, white-light Carrington maps of the streamer belt, and the magnetic sector structure measured in situ by Parker Solar Probe and 1 au spacecraft. We study these metrics for potential field source surface (PFSS) models as a function of source surface height and magnetogram choice, as well as comparing to the more physical Wang–Sheeley–Arge (WSA) and the Magnetohydrodynamic Algorithm outside a Sphere (MAS) models. We find that simultaneous optimization of PFSS models to all three metrics is not currently possible, implying a trade-off between the quality of representation of coronal holes and streamer belt topology. WSA and MAS results show the additional physics that they include address this by flattening the streamer belt while maintaining coronal hole sizes, with MAS also improving coronal hole representation relative to WSA. We conclude that this framework is highly useful for inter- and intra-model comparisons. Integral to the framework is the standardization of observables required of each model, evaluating different model aspects.

S. T. Badman↗

Statistical Wind-Tunnel Experimentation Advancements for eVTOL Aircraft Aero-Propulsive Model Development

Compared to many conventional aircraft, electric vertical takeoff and landing (eVTOL) vehicle designs present greater aero-propulsive complexity and many interacting factors requiring development of novel testing and model development strategies. This paper describes several advances in statistically-rigorous experiment design methods for wind-tunnel testing of eVTOL aircraft enabling accurate aero-propulsive model development. The techniques are applied to a subscale tilt-rotor eVTOL aircraft configuration built at NASA Langley Research Center that has 24 independent control effectors. The general testing approach, leveraging design of experiments and response surface methodology, is outlined and followed by a detailed description of the experiment design and modeling approach employed for different phases of the wind-tunnel testing. A new approach for gravitational tare modeling is developed and validated, which reduces the required wind-tunnel test time by nearly 50%. A mostly automated trim envelope determination approach is formulated to rapidly determine the aircraft flight envelope to inform transition testing for aero-propulsive model development. Efficient powered-airframe characterization tests are conducted using a nested I-optimal design approach formulated in previous simulation-based testing. The collective wind-tunnel testing approach facilitates rapid characterization of eVTOL aircraft configurations and produces accurate aero-propulsive models that can be used in flight dynamics simulations.

design of experiments↗