Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “coordinate 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,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A scalable planning framework of energy storage systems under frequency dynamics constraints

As the penetration of renewables increases in power systems, the declining system inertia can cause frequency stability issues. Battery energy storage systems (BESSs) respond fast and therefore can relieve the low inertia difficulty but need to be appropriately sized considering the associated cost. This paper presents a novel stochastic optimization model for economically planning BESS capacity while considering the spatial–temporal correlation of wind generation and generator outages under frequency stability constraints, which include the rate-of-change of frequency (RoCoF), frequency nadir (FN), and quasi-steady-state (QSS) frequency. A set of new FN constraints that can be easily linearized is developed. To account for renewable uncertainties, a realistic uncertainty modeling approach, Random Field, is adopted to generate wind generation scenarios by considering both spatial and temporal evolutions of wind speed profiles. The ESS sizing is formulated as a mixed-integer linear programming problem and solved by using a scalable decomposition-and-coordination approach, Surrogate Absolute Value Lagrangian Relaxation (SAVLR). To further improve the scalability and reduce computational burdens, a rolling-horizon-based update is developed and incorporated into SAVLR for providing a practical solution to the long-term planning of very large-scale power systems. Finally, a modified IEEE 118-bus system and the Polish system are used to validate the effectiveness and scalability of the model and solution methodology.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS2017

Current Safety Nets Within the U.S. National Airspace System

There are over 70,000 flights managed per day in the National Airspace System, with approximately 7,000 aircraft in the air over the United States at any given time. Operators of each of these flights would prefer to fly a user-defined “4D” trajectory (4DT), which includes arrival and departure times; preferred gates and runways at the airport; efficient, wind-optimal routes for departure, cruise and arrival phase of flight; and fuel efficient altitude profiles. To demonstrate the magnitude of this achievement a single flight from Los Angeles to Baltimore, accesses over 35 shared or “constrained” resources that are managed by roughly 30 air traffic controllers (at towers, approach control and en route sectors); along with traffic managers at 12 facilities, using over 22 different, independent automation system (including TBFM, ERAM, STARS, ASDE-X, FSM, TSD, GPWS, TCAS, etc.). In addition, dispatchers, ramp controllers and others utilize even more systems to manage each flight’s access to operator-managed resources. Flying an ideal 4DT requires successful coordination of all flight constraints among all flights, facilities, operators, pilots and controllers. Additionally, when conditions in the NAS change, the trajectories of one or more aircraft may need to be revised to avoid loss of flight efficiency, predictability, separation or system throughput. The Aviation Safety Network has released the 2016 airliner accident statistics showing a very low total of 19 fatal airliner accidents, resulting in 325 fatalities1. Despite several high profile accidents, the year 2016 turned out to be a very safe year for commercial aviation, Aviation Safety Network data show. Over the year 2016 the Aviation Safety Network recorded a total of 19 fatal airliner accidents [1], resulting in 325 fatalities. This makes 2016 the second safest year ever, both by number of fatal accidents as well as in terms of fatalities. In 2015 ASN recorded 16 accidents while in 2013 a total of 265 lives were lost. How can we keep it that way and not upset the apple cart by premature insertion of innovative technologies, functions, and procedures? In aviation, safety nets function as the last system defense against incidents and accidents. Current ground-based and airborne safety nets are well established and development to make them more efficient and reliable continues. Additionally, future air traffic control safety nets may emerge from new operational concepts.

Smith, Brian E.↗
NASA NTRS2005

Cognitive Functioning in Space Exploration Missions: A Human Requirement

Solving cognitive issues in the exploration missions will require implementing results from both Human Behavior and Performance, and Space Human Factors Engineering. Operational and research cognitive requirements need to reflect a coordinated management approach with appropriate oversight and guidance from NASA headquarters. First, this paper will discuss one proposed management method that would combine the resources of Space Medicine and Space Human Factors Engineering at JSC, other NASA agencies, the National Space Biomedical Research Institute, Wyle Labs, and other academic or industrial partners. The proposed management is based on a Human Centered Design that advocates full acceptance of the human as a system equal to other systems. Like other systems, the human is a system with many subsystems, each of which has strengths and limitations. Second, this paper will suggest ways to inform exploration policy about what is needed for optimal cognitive functioning of the astronaut crew, as well as requirements to ensure necessary assessment and intervention strategies for the human system if human limitations are reached. Assessment strategies will include clinical evaluation and fitness-to-perform evaluations. Clinical intervention tools and procedures will be available to the astronaut and space flight physician. Cognitive performance will be supported through systematic function allocation, task design, training, and scheduling. Human factors requirements and guidelines will lead to well-designed information displays and retrieval systems that reduce crew time and errors. Means of capturing process, design, and operational requirements to ensure crew performance will be discussed. Third, this paper will describe the current plan of action, and future challenges to be resolved before a lunar or Mars expedition. The presentation will include a proposed management plan for research, involvement of various organizations, and a timetable of deliverables.

Fiedler, Edan↗
NASA NTRS2008

The Whole Heliosphere Interval: Campaign Summaries and Early Results

The Whole Heliosphere Interval (WHI) is an internationally coordinated observing and modeling effort to characterize the 3-dimensional interconnected solar-heliospheric-planetary system - a.k.a. the "heliophysical" system. The heart of the WHI campaign is the study of the interconnected 3-D heliophysical domain, from the interior of the Sun, to the Earth, outer planets, and into interstellar space. WHI observing campaigns began with the 3-0 solar structure from solar Carrington Rotation 2068, which ran from March 20 - April 16, 2008. Observations and models of the outer heliosphere and planetary impacts extended beyond those dates as necessary; for example, the solar wind transit time to outer planets can take months. WHI occurs during solar minimum, which optimizes our ability to characterize the 3-D heliosphere and trace the structure to the outer limits of the heliosphere. A summary of some of the key results from the WHI first workshop in August 2008 will be given.

Thompson, Barbara J.↗
DOE OSTI2026

Group-Additivity–Embedded Multiscale Modeling for Electric Field-Enhanced Nanocatalysis

Elucidating structure-performance relationships remains a central challenge in field-enhanced catalysis, where nanoparticles exhibit nonuniform surface sites with site-dependent responses to electric fields. Low-coordination sites (edges, corners, and tips) are particularly electric field-sensitive (EF), leading to nonuniform charge distribution, adsorption energies, and catalytic activity. Here, using ammonia decomposition on a ruthenium cluster as a model system, we develop a transferable multiscale framework integrating density functional theory, group additivity (GA), Brønsted-Evans-Polanyi scaling, and microkinetic modeling to predict EF-dependent activity across nonuniform cluster sites. Across sites and fields, the nitrogen adsorption energy (E N ) emerges as the governing descriptor, yielding robust volcano relationships whose optimum shifts systematically with field: negative fields strengthen N binding via electron accumulation, while positive fields weaken N binding via charge depletion, moving the optimal E N toward weaker binding. Microkinetic analysis shows that N≡N bond formation remains the key kinetic bottleneck over most conditions; positive fields lower the effective barrier and, critically, increase the fraction of near-optimal active sites, leading to a net enhancement in overall activity relative to zero-field and negative-field cases. By capturing EF- and site-dependent energetics with high accuracy and low computational cost, this GA-embedded multi-scale simulation workflow provides a physically interpretable route to predict and design field-enhanced nanocatalysis.

ammonia decomposition↗
DOE OSTI2023

All-temperature zinc batteries with high-entropy aqueous electrolyte

Electrification of transportation and rising demand for grid energy storage continue to build momentum around batteries across the globe. However, the supply chain of Li-ion batteries is exposed to the increasing challenges of resourcing essential and scarce materials. Therefore, incentives to develop more sustainable battery chemistries are growing. Here, in this paper, we show an aqueous ZnCl 2 electrolyte with introduced LiCl as supporting salt. Once the electrolyte is optimized to Li 2 ZnCl 4 ∙9H 2 O, the assembled Zn–air battery can sustain stable cycling over the course of 800 hours at a current density of 0.4 mA cm -2 between -60 °C and +80 °C, with 100% Coulombic efficiency for Zn stripping/plating. Even at -60 °C, >80% of room-temperature power density can be retained. Advanced characterization and theoretical calculations reveal a high-entropy solvation structure that is responsible for the excellent performance. The strong acidity allows ZnCl 2 to accept donated Cl - ions to form ZnCl 4 2- anions, while water molecules remain within the free solvent network at low salt concentration or coordinate with Li ions. Our work suggests an effective strategy for the rational design of electrolytes that could enable next-generation Zn batteries.

25 ENERGY STORAGE↗
NASA NTRS1997

Application of PML Absorbing Boundary Conditions to the Benchmark Problems of Computational Aeroacoustics

Accurate numerical non-reflecting boundary conditions are important in all the proposed benchmark problems of the Second Workshop. Recently, a new absorbing boundary condition has been developed using Perfectly Matched Layer (PML) equations for the Euler equations. In this approach, a region with a width of a few grid points is introduced adjacent to the non-reflecting boundaries. In the added region, Perfectly Matched Layer equations are constructed and applied so that the out-going waves are absorbed inside the layer with little reflection to the interior domain. It will be demonstrated in the present paper that the proposed absorbing boundary condition is quite general and versatile, applicable to radiation boundaries as well as inflow and outflow boundaries. It is also easy to implement. The emphasis of the paper will be on the application of the PML absorbing boundary condition to problems in Categories 1, 2, and 3. In Category 1, solutions of problems 1 and 2 are presented. Both problems are solved using a multi-domain polar grid system. Perfectly Matched Layer equations for a circular boundary are constructed and their effectiveness assessed. In Category 2, solutions of problem 2 are presented. Here, in addition to the radiation boundary conditions at the far field in the axisymmetric coordinate system, the inflow boundary condition at the duct inlet is also dealt with using the proposed Perfectly Match Layer equations. At the inlet, a PML domain is introduced in which the incident duct mode is simulated while the waves reflected from the open end of the duct are absorbed at the same time. In Category 3, solutions of all three problems are presented. Again, the PML absorbing boundary condition is used at the inflow boundary so that the incoming vorticity wave is simulated while the outgoing acoustic waves are absorbed with very little numerical reflection. All the problems are solved using central difference schemes for spatial discretizations and the optimized Low-Dissipation and Low-Dispersion Runge-Kutta scheme for the time integration. Issues of numerical accuracy and efficiency are also addressed.

Hu, Fang Q.↗
DOE OSTI2021

Addressing Human and Organizational Factors In Nuclear Industry Modernization: A Sociotechnically-Based Strategic Framework

The modernization of nuclear power plants will require an advanced concept of operations, involving an integrated set of tightly coupled systems in which all stakeholders act in a coordinated manner. For this modernization effort to be enabled, we developed a human and organizational factors approach, based on a broad sociotechnical framework. Starting from core human factors principles, we conducted a literature review of the methods and approaches relevant to the modernization problem. These included core disciplines such as cognitive systems engineering, systems theoretic accident modeling and processes (STAMP), human systems integration, resilience engineering, and macroergonomics but also related topics of safety culture and organizational change. From this literature, we developed a conceptual framework centered around the work system with its four interacting components: people, technology, process, and governance. In an effective work system, these four components are jointly optimized according to three systems criteria: efficiency, effectiveness, and safety. System failure may result from excessive emphasis on any one criterion. The actual work of attaining joint optimization in a given work system can be accomplished by utilizing three high level functions: knowledge elicitation, knowledge representation, and cross-functional integration. We illustrated the utility of this approach by applying it to practical problems and case studies. The complete report is available in [1]; this paper is a description of our approach and report contents.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2025

Methods and Tools To Assess Robustness of Nuclear Plant Outages

Refueling outages are one of the most challenging phases in a nuclear power plant (NPP) operating cycle. Refueling outages are extremely costly for an NPP due to the large amount of required resources and because of lost revenue due to plant being off the grid. Outage durations have steadily decreased across the industry over that last few decades primarily due to improved planning and coordination, but there are still many plants that struggle to meet the performance metrics accomplished by other utilities. Schedule resilience is one of the issues. NPP outages require scheduling thousands of activities within 30 days on average. Despite detailed planning, once the outage starts, numerous emergent issues typically appear along with schedule delays requiring continuous replanning and adjustment. When schedule disruption occurs during an outage, plant staff make urgent efforts to recover but are often not able to maintain the planned outage duration. These outage delays can cost a utility several million dollars per day. Tools that could help outage schedulers create a more resilient schedule and allow them to optimally reschedule emergent work could significantly reduce outage delays. One key aspect of creating a resilient schedule is to have accurate estimates for activity duration. Another important outage scheduling capability is the ability to schedule emergent work with minimal disruption. This paper focuses on developing tools and methods to support NPPs with outage schedule optimization and it describes the initial development of tools to support outage management that leverage computational and machine learning methods.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Addressing Human and Organizational Factors in Nuclear Industry Modernization: A Sociotechnically Based Strategic Framework

The modernization of nuclear power plants will require an advanced concept of operations, involving an integrated set of tightly coupled systems in which all stakeholders act in a coordinated manner. For this modernization effort to be enabled, we developed a human and organizational factors approach based on a broad sociotechnical framework. Starting from core human factors principles, we conducted a literature review of the methods and approaches relevant to the modernization problem. These included not only core disciplines such as cognitive systems engineering, systems theoretic accident modeling and processes, human systems integration, resilience engineering, and macroergonomics but also related topics of safety culture and organizational change. From this literature, we developed a conceptual framework centered around the work system with its four interacting components: people, technology, process, and governance. In an effective work system, these four components are jointly optimized according to three systems criteria: efficiency, effectiveness, and safety. System failure may result from excessive emphasis on any one criterion. The actual work of attaining joint optimization in a given work system can be accomplished by utilizing three high-level functions: knowledge elicitation, knowledge representation, and cross-functional integration. Finally, we illustrated the utility of this approach by applying it to practical problems and case studies.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Enabling Modular Autonomous Feedback‐Loops in Materials Science through Hierarchical Experimental Laboratory Automation and Orchestration

Abstract Materials acceleration platforms (MAPs) operate on the paradigm of integrating combinatorial synthesis, high‐throughput characterization, automatic analysis, and machine learning. Within a MAP, one or multiple autonomous feedback loops may aim to optimize materials for certain functional properties or to generate new insights. The scope of a given experiment campaign is defined by the range of experiment and analysis actions that are integrated into the experiment framework. Herein, the authors present a method for integrating many actions within a hierarchical experimental laboratory automation and orchestration (HELAO) framework. They demonstrate the capability of orchestrating distributed research instruments that can incorporate data from experiments, simulations, and databases. HELAO interfaces laboratory hardware and software distributed across several computers and operating systems for executing experiments, data analysis, provenance tracking, and autonomous planning. Parallelization is an effective approach for accelerating knowledge generation provided that multiple instruments can be effectively coordinated, which the authors demonstrate with parallel electrochemistry experiments orchestrated by HELAO. Efficient implementation of autonomous research strategies requires device sharing, asynchronous multithreading, and full integration of data management in experimental orchestration, which to the best of the authors’ knowledge, is demonstrated for the first time herein.

36 MATERIALS SCIENCE↗
NASA NTRS1994

Optimizing the Zeldovich approximation

We have recently learned that the Zeldovich approximation can be successfully used for a far wider range of gravitational instability scenarios than formerly proposed; we study here how to extend this range. In previous work (Coles, Melott and Shandarin 1993, hereafter CMS) we studied the accuracy of several analytic approximations to gravitational clustering in the mildly nonlinear regime. We found that what we called the 'truncated Zeldovich approximation' (TZA) was better than any other (except in one case the ordinary Zeldovich approximation) over a wide range from linear to mildly nonlinear (sigma approximately 3) regimes. TZA was specified by setting Fourier amplitudes equal to zero for all wavenumbers greater than k(sub nl), where k(sub nl) marks the transition to the nonlinear regime. Here, we study the cross correlation of generalized TZA with a group of n-body simulations for three shapes of window function: sharp k-truncation (as in CMS), a tophat in coordinate space, or a Gaussian. We also study the variation in the crosscorrelation as a function of initial truncation scale within each type. We find that k-truncation, which was so much better than other things tried in CMS, is the worst of these three window shapes. We find that a Gaussian window e(exp(-k(exp 2)/2k(exp 2, sub G))) applied to the initial Fourier amplitudes is the best choice. It produces a greatly improved crosscorrelation in those cases which most needed improvement, e.g. those with more small-scale power in the initial conditions. The optimum choice of kG for the Gaussian window is (a somewhat spectrum-dependent) 1 to 1.5 times k(sub nl). Although all three windows produce similar power spectra and density distribution functions after application of the Zeldovich approximation, the agreement of the phases of the Fourier components with the n-body simulation is better for the Gaussian window. We therefore ascribe the success of the best-choice Gaussian window to its superior treatment of phases in the nonlinear regime. We also report on the accuracy of particle positions and velocities produced by TZA.

Melott, Adrian L.↗
DOE OSTI2022

Revisiting metastable immiscibility in SiO 2 –Al 2 O 3 : Structure and phase separation of supercooled liquids and glasses

Understanding and controlling liquid–liquid phase separation in aluminosilicates is crucial for optimizing glass properties. However, the metastable nature of aluminosilicates’ phase separation has made it difficult to study experimentally, and uncertainty persists regarding the compositional and temperature extents of the miscibility gap. Here, we present new experimental evidence that suggests a consolute temperature between 1440 and 1590°C and endmember compositions of 7 and 62 mol.% Al 2 O 3 for the phase-separated glasses. Using containerless melt processing, deeply supercooled liquids over the 0–60 mol.% Al 2 O 3 range are probed with in situ small- and wide-angle X-ray scattering, which simultaneously reveals changes in nanoscale density heterogeneity and atomic structure. Correlations between phase separation and atomic coordination environments are compared for liquids and glasses. Pair distribution function analysis shows mean O–(Si + Al) coordination increases with Al 2 O 3 content and decreases with temperature.

36 MATERIALS SCIENCE↗
NASA NTRS2003

Terrain Portrayal for Head-Down Displays Flight Test

The Synthetic Vision Systems General Aviation (SVS-GA) element of NASA's Aviation Safety Program is developing technology to eliminate low visibility induced General Aviation (GA) accidents through the application of synthetic vision techniques. SVS displays present computer generated 3-dimensional imagery of the surrounding terrain to greatly enhance pilot's situation awareness (SA), reducing or eliminating Controlled Flight into Terrain (CFIT), as well as Low-Visibility Loss of Control (LVLOC) accidents. In addition to substantial safety benefits, SVS displays have many potential operational benefits that can lead to flight in instrument meteorological conditions (IMC) resembling those conducted in visual meteorological conditions (VMC). Potential benefits could include lower landing minimums, more approach options, reduced training time, etc. SVS conducted research will develop display concepts providing the pilot with an unobstructed view of the outside terrain, regardless of weather conditions and time of day. A critical component of SVS displays is the appropriate presentation of terrain to the pilot. The relationship between the realism of the terrain presentation and resulting enhancements of pilot SA and pilot performance has been largely undefined. Comprised of coordinated simulation and flight test efforts, the terrain portrayal for head-down displays (TP-HDD) test series examined the effects of two primary elements of terrain portrayal: variations of digital elevation model (DEM) resolution and terrain texturing. Variations in DEM resolution ranged from sparsely spaced (30 arc-sec/2,953ft) to very closely spaced data (1 arc-sec/98 ft). Variations in texture involved three primary methods: constant color, elevation-based generic, and photo-realistic, along with a secondary depth cue enhancer in the form of a fishnet grid overlay. The TP-HDD test series was designed to provide comprehensive data to enable design trades to optimize all SVS applications, as well as develop requirements and recommendations to facilitate the implementation and certification of SVS displays. The TP-HDD flight experiment utilized the NASA LaRC Cessna 206 Stationaire and evaluated eight terrain portrayal concepts in an effort to confirm and extend results from the previously conducted TP-HDD simulation experiment. A total of 15 evaluation pilots, of various qualifications, accumulated over 75 hours of dedicated research flight time at Newport News (PHF) and Roanoke (ROA), VA, airports from August through October, 2002. This report will present results from the portion of testing conducted at Roanoke, VA.

Hughes, Monica F.↗
DOE OSTI2021

Understanding cation-disordered rocksalt oxyfluoride cathodes

Partial fluorine (F) substitution into the oxygen (O) sublattice has been shown to improve cycling stability of cation-disordered Li-excess rocksalt oxide (DRX) cathodes. Detailed understanding on failure mechanisms and key optimization knowledge of fluorinated-DRX (F-DRX), however, are lacking. In the present study, we incorporate different amounts of F into a baseline DRX system, Li 1.2 Ti 0.4 Mn 0.4 O 2 (LTMO3.0, 3.0 denotes nominal Li/Mn ratio), and synthesize two oxyfluoride compounds, Li-rich Li 1.3 Ti 0.3 Mn 0.4 O 1.7 F 0.3 (LTMOF3.25) and Mn-rich Li 1.2 Ti 0.2 Mn 0.6 O 1.8 F 0.2 (LTMOF2.0) with an increased and reduced Li/Mn ratio of 3.25 and 2.0, respectively. Through careful monitoring of chemical and structural evolution, we show that cycling-induced changes are manifested not only by Mn reduction and degradation of its local coordination environment, but also by F enrichment and formation of LiF-type of domains on the surface. A “concerted-densification” based failure mechanism, involving atomic-level changes in both transition-metal cationic sublattice and oxygen/fluorine anionic sublattice, is proposed for the degradation in F-DRX cathode materials. The study reveals that increasing F content accompanied by reduced Li/Mn ratio mitigates the degradation process, offering key design strategies in achieving balanced cathode capacity and stability.

25 ENERGY STORAGE↗
DOE OSTI2026

4D beam matrix reconstruction in particle accelerators

Transverse beam parameters in particle accelerators are commonly described using the Twiss parameters, which are experimentally accessible yet inherently limited because they neglect correlations between different transverse coordinates. Such correlations frequently arise from uncompensated cathode magnetic fields or misaligned focusing quadrupoles, affecting beam quality and accelerator performance. To address this limitation, we propose and validate a novel diagnostic method for the complete four-dimensional (4D) transverse beam matrix. Our method involves passing the beam through a beamline comprising both conventional and skew quadrupole magnets, followed by downstream measurements of the resulting two-dimensional (2D) beam profiles. These measurements represent distinct 2D projections of the underlying 4D transverse phase–space distribution. By systematically varying quadrupole strengths, multiple independent projections of the beam phase space are obtained. We reconstruct the original 4D beam matrix from these measured projections using an optimization-based least-square fit, providing fast and robust reconstruction regardless of the specific beamline configuration. Through extensive numerical simulations and realistic particle-tracking studies, we demonstrate the diagnostic’s accuracy, robustness, and capability to achieve reconstruction uncertainties smaller than measurement errors, particularly when employing sufficient numbers of quadrupole scans. This method presents a powerful and flexible approach for comprehensive beam characterization and accelerator tuning.

43 PARTICLE ACCELERATORS↗
NASA NTRS1995

Science Rationale for a Micro-Met Mission to Augment InterMarsNet

The 2003 opportunity has the potential to carry out for the first time in Mars exploration history, coordinated measurements from the surface and from orbit that can address fundamental issues associated with the Martian global circulation and climate system. Coordinated measurements are defined here to mean collecting meteorological data from a network of 12-16 globally distributed surface stations simultaneously with an orbiter carrying an atmospheric sounder. With such measurements it is possible to define the horizontally varying (barotropic) and vertically varying (baroclinic) components of the global circulation from which the full 3-dimensional horizontal wind field can be reconstructed. It is also possible to precisely define the CO2 cycle, the main component of the current climate system. InterMarsNet, as currently envisioned, consists of 3-4 landers with a supporting communications orbiter that may carry some instrumentation. The landers are likely to touch down in low latitudes and will probably be configured to optimize seismological objectives. We propose to augment the InterMarsNet meteorological objectives by flying an additional 10-15 "MicroMet" landers equipped to measure surface pressure and nothing else. Surface pressure is the most fundamental meteorological parameter and it is the easiest to measure. The sensors are light, operate with minimal power, and do not require orientation or deployment. Consequently the landers can be very small (< 10 kg) and 10-15 of them can be delivered by a Med-lite launcher. This would enable global network science for meteorology and, when combined wAh the more sophisticated measurements from the InterMarsNet landers and the simultaneous temperature and dust profiling measurements from an atmospheric sounder aboard the orbiter, it would enable us to reconstruct global scale circulation patterns.

Haberle, Robert M.↗
NASA NTRS2010

Adaptive Numerical Algorithms in Space Weather Modeling

Space weather describes the various processes in the Sun-Earth system that present danger to human health and technology. The goal of space weather forecasting is to provide an opportunity to mitigate these negative effects. Physics-based space weather modeling is characterized by disparate temporal and spatial scales as well as by different physics in different domains. A multi-physics system can be modeled by a software framework comprising of several components. Each component corresponds to a physics domain, and each component is represented by one or more numerical models. The publicly available Space Weather Modeling Framework (SWMF) can execute and couple together several components distributed over a parallel machine in a flexible and efficient manner. The framework also allows resolving disparate spatial and temporal scales with independent spatial and temporal discretizations in the various models. Several of the computationally most expensive domains of the framework are modeled by the Block-Adaptive Tree Solar wind Roe Upwind Scheme (BATS-R-US) code that can solve various forms of the magnetohydrodynamics (MHD) equations, including Hall, semi-relativistic, multi-species and multi-fluid MHD, anisotropic pressure, radiative transport and heat conduction. Modeling disparate scales within BATS-R-US is achieved by a block-adaptive mesh both in Cartesian and generalized coordinates. Most recently we have created a new core for BATS-R-US: the Block-Adaptive Tree Library (BATL) that provides a general toolkit for creating, load balancing and message passing in a 1, 2 or 3 dimensional block-adaptive grid. We describe the algorithms of BATL and demonstrate its efficiency and scaling properties for various problems. BATS-R-US uses several time-integration schemes to address multiple time-scales: explicit time stepping with fixed or local time steps, partially steady-state evolution, point-implicit, semi-implicit, explicit/implicit, and fully implicit numerical schemes. Depending on the application, we find that different time stepping methods are optimal. Several of the time integration schemes exploit the block-based granularity of the grid structure. The framework and the adaptive algorithms enable physics based space weather modeling and even forecasting.

Toth, Gabor↗