Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “global 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,459 records · Page 81

A processing centre for the CNES CE-GPS experimentation

CNES is involved in a GPS (Global Positioning System) geostationary overlay experimentation. The purpose of this experimentation is to test various new techniques in order to select the optimal station synchronization method, as well as the geostationary spacecraft orbitography method. These new techniques are needed to develop the Ranging GPS Integrity Channel services. The CNES experimentation includes three transmitting/receiving ground stations (manufactured by IN-SNEC), one INMARSAT 2 C/L band transponder and a processing center named STE (Station de Traitements de l'Experimentation). Not all the techniques to be tested are implemented, but the experimental system has to include several functions; part of the future system simulation functions, such as a servo-loop function, and in particular a data collection function providing for rapid monitoring of system operation, analysis of existing ground station processes, and several weeks of data coverage for other scientific studies. This paper discusses system architecture and some criteria used in its design, as well as the monitoring function, the approach used to develop a low-cost and short-life processing center in collaboration with a CNES sub-contractor (ATTDATAID), and some results.

Suard, Norbert↗

Flight data results of estimate fusion for spacecraft rendezvous navigation from shuttle mission STS-69

A recently developed rendezvous navigation fusion filter that optimally exploits existing distributed filters for rendezvous and GPS navigation to achieve the relative and inertial state accuracies of both in a global solution is utilized here to process actual flight data. Space Shuttle Mission STS-69 was the first mission to date which gathered data from both the rendezvous and Global Positioning System filters allowing, for the first time, a test of the fusion algorithm with real flight data. Furthermore, a precise best estimate of trajectory is available for portions of STS-69, making possible a check on the performance of the fusion filter. In order to successfully carry out this experiment with flight data, two extensions to the existing scheme were necessary: a fusion edit test based on differences between the filter state vectors, and an underweighting scheme to accommodate the suboptimal perfect target assumption made by the Shuttle rendezvous filter. With these innovations, the flight data was successfully fused from playbacks of downlinked and/or recorded measurement data through ground analysis versions of the Shuttle rendezvous filter and a GPS filter developed for another experiment. The fusion results agree with the best estimate of trajectory at approximately the levels of uncertainty expected from the fusion filter's covariance matrix.

Carpenter, J. Russell↗

Performance Analysis and Portability of the PLUM Load Balancing System

The ability to dynamically adapt an unstructured mesh is a powerful tool for solving computational problems with evolving physical features; however, an efficient parallel implementation is rather difficult. To address this problem, we have developed PLUM, an automatic portable framework for performing adaptive numerical computations in a message-passing environment. PLUM requires that all data be globally redistributed after each mesh adaption to achieve load balance. We present an algorithm for minimizing this remapping overhead by guaranteeing an optimal processor reassignment. We also show that the data redistribution cost can be significantly reduced by applying our heuristic processor reassignment algorithm to the default mapping of the parallel partitioner. Portability is examined by comparing performance on a SP2, an Origin2000, and a T3E. Results show that PLUM can be successfully ported to different platforms without any code modifications.

Oliker, Leonid↗

ViDI: Virtual Diagnostics Interface: The Future of Wind Tunnel Testing - Volume 1

The quality of data acquired in a given test facility ultimately resides within the fidelity and implementation of the instrumentation systems. Over the last decade, the emergence of robust optical techniques has vastly expanded the envelope of measurement possibilities. At the same time the capabilities for data processing, data archiving and data visualization required to extract the highest level of knowledge from these global, on and off body measurement techniques have equally expanded. Yet today, while the instrumentation has matured to the production stage, an optimized solution for gaining knowledge from the gigabytes of data acquired per test (or even per test point) is lacking. A technological void has to be filled in order to possess a mechanism for near-real time knowledge extraction during wind tunnel experiments. Under these auspices, the Virtual Diagnostics Interface, or ViDI, was developed.

Fleming, Gary A.↗

Industry Consolidation and Future Airline Network Structures in Europe

In the current downturn in demand for air travel, major airlines are revising and rationalizing their networks in an attempt to improve financial performance and strengthen their defences against both new entrants and traditional rivals. Expansion of commercial agreements or alliances with other airlines has become a key reaction to the increasingly competitive marketplace. In the absence, for regulatory reasons, of cross-border mergers these are the principal means by which the industry can consolidate internationally. This paper analyzes the developments which have been taking place and attempts to itentify the implications for airline network structures and the function of different hub airports. The range of services available to passengers in long-haul markets to/from Europe is evaluated before and after recent industry reorganization. Hubs are crucial to interlink the route networks of parmers in an alliance. However, duplication between nearby hub airports that find themselves within the same airline alliance can lead to loss of service at the weaker locations. The extent to which the alliance hubs in Europe duplicate or complement each other in terms of network coverage is assessed and this methodology also enables the optimal partnerships for "unattached" airlines to be identified. The future role of the various European hubs is considered under different scenarios of global alliance development. The paper concludes by considering possible longer-term developments. In an environment where the low-cost carriers will provide a major element of customer choice, it is suggested that the traditional airlines will retrench around their hubs, surrendering many secondary cities to the low-cost sector. Further reduction in the number of alliances could threaten more of the European hubs. For both regulatory and commercial reasons, the end result may be just one airline alliance - so recreating in the deregulated market the historic rule of IATA.

Dennis, Nigel↗

Equuleus Launch Window Analysis and Mission Design

This paper presents the trajectory design process for EQUULEUS, a 6U CubeSat developed by JAXA and the University of Tokyo that is scheduled to launch as a piggyback of NASA’s Artemis 1. After separation from the upper stage of the Space Launch System, EQUULEUS will maneuver along a low-energy transfer to an Earth–Moon quasi-rectilinear halo orbit in 1-to-4 resonance with the lunar synodic period. As a secondary payload, the trajectory of EQUULEUS needs to be compatible with the requirements of the primary mission, but also robust against disturbances and potential changes in the deployment state. Realistic initial conditions spanning two years of potential launch windows are processed and the solution structure for optimal lunar transfers is analyzed. A host of candidate solutions is presented, compatibly with the fuel and power limitations of EQUULEUS. The global understanding of the solution space is shown to be insightful for the design of robust trajectories for limited control-authority spacecraft.

Kawakatsu, Yasuhiro↗

Dynamically-Spaced Geo-Grid Segmentation for Weighted Point Sampling on a Polygon Map Layer

Geo-grid algorithms divide a large polygon area into several smaller polygons, which are important for studying or executing a set of operations on underlying topological features of a map. The current geo-grid algorithms divide a large polygon in to a set of smaller but equal size polygons only (e.g. is ArcMaps Fishnet). The time to create a geo-grid is typically proportional to number of smaller polygons created. This raises two problems - (i) They cannot skip unwanted areas (such as water bodies, given about 71% percent of the Earth's surface is water-covered); (ii) They are incognizant to any underlying feature set that requires more deliberation. In this work, we propose a novel dynamically spaced geo-grid segmentation algorithm that overcomes these challenges and provides a computationally optimal output for borderline cases of an uneven polygon. Our method uses an underlying topological feature of population distributions, from the LandScan Global 2016 dataset, for creating grids as a function of these weighted features. We benchmark our results against available algorithms and found our approach improves geo-grid creation. Later on, we demonstrate the proposed approach is more effective in harvesting Points of Interest data from a crowd-sourced platform.

Sims, Kelly↗

WE-H-201-00: Opportunities for Physicists to Support Low and Mid-Income Countries

The desperate need for radiotherapy in low and mid-income countries (LMICs) has been well documented. Roughly 60 % of the worldwide incidence of cancer occurs in these resource-limited settings and the international community alongside governmental and non-profit agencies have begun publishing reports and seeking help from qualified volunteers. However, the focus of several reports has been on how dire the situation is and the magnitude of the problem, leaving most to feel overwhelmed and unsure as to how to help and why to get involved. This session will help to explain the specific ways that Medical Physicists can uniquely assist in this grand effort to help bring radiotherapy to grossly-underserved areas. Not only can these experts fulfill an important purpose, they also can benefit professionally, academically, emotionally and socially from the endeavor. By assisting others worldwide with their skillset, Medical Physicists can end up helping themselves. Learning Objectives: Understand the need for radiotherapy in LMICs. Understand which agencies are seeking Medical Physicists for help in LMICs. Understand the potential research funding mechanisms are available to establish academic collaborations with LMIC researchers/physicians. Understand the potential social and emotional benefits for both the physicist and the LMIC partners when collaborations are made. Understand the potential for collaboration with other high-income scientists that can develop as the physicist partners with other large institutions to assist LMICs. Wil Ngwa - A recent United Nations Study reports that in developing countries more people have access to cell phones than toilets. In Africa, only 63% of the population has access to piped water, yet, 93% of Africans have cell phone service. Today, these cell phones, Skype, WhatsApp and other information and communication technologies (ICTs) connect us in unprecedented ways and are increasingly recognized as powerful, indispensable to global health. Thanks to ICTs, there are growing opportunities for Medical Physicists to reach out beyond the bunker and impact the world far beyond, without even having to travel. These growing opportunities in global health for Medical Physicists, powered by ICTs, will be highlighted in this presentation, illustrated by high impact examples/models across the globe that are improving patient safety and healthcare outcomes, saving lives. Learning Objectives: Published definitions of global health and the emerging field of global radiation oncology Learn about the transformative potential of ICTs in global radiation oncology care, research and education with focus on Medical Physics Learn about high impact examples of ICT-powered global radiation oncology and the increasing opportunities for participation by Medical Physicists. Yakov Pipman - The number and scope of volunteer Medical Physics activities in support of low-to-middle income countries has been increasing gradually. This happens through a variety of formal channels and to some extent through less formal but personal initiatives. A good deal of effort is dedicated by many, but many more could be recruited through a structured framework to volunteer. We will look into typical volunteer activities and how they fit with organizations already involved in advancing Medical Physics in LMIC. We will identify the range of these organizational activities and their scope to reveal areas of further need. We will point to a few key features of MPWB ( http://www.mpwb.org ) as a volunteering and collaborating structure and how members can get involved and contribute to these efforts at the grass roots level. The goal is that scarce resources can thus be channeled to complement rather than compete with those already in place. Learning Objectives: Understand the strengths and limitations of various organizations that support Medical Physics efforts in LMIC. Learn about ways to volunteer and contribute to Global Health through a grass roots organization focused on Medical Physics in LMIC. Perry Sprawls - With the growing capability and complexity of medical imaging methods in all countries of the world, the expanding role of medical physicists includes optimizing imaging procedures with respect to image quality, radiation dose, and other conflicting factors. With access to appropriate educational resources local medical physicists in all countries can provide direct clinical support and educational for other medical professionals. This is being supported through the process of Collaborative Teaching that combines the capabilities and experience of medical physicists in countries spanning the spectrum of economic, technological, and clinical development. The supporting resources are on the web at: http://www.sprawls.org/resources . Learning Objectives: Identify the medical physics educational needs to support effective and optimized medical imaging procedures. Use collaborative teaching resources to enhance the role of medical physicists in all countries of the world.

60 APPLIED LIFE SCIENCES↗

Petascale solvers for anisotropic PDEs in atmospheric modelling on GPU clusters

Memory bound applications such as solvers for large sparse systems of equations remain a challenge for GPUs. Fast solvers should be based on numerically efficient algorithms and implemented such that global memory access is minimised. To solve systems with trillions ( O ( 10 12 ) ) unknowns the code has to make efficient use of several million individual processor cores on large GPU clusters. We describe the multi-GPU implementation of two algorithmically optimal iterative solvers for anisotropic PDEs which are encountered in (semi-) implicit time stepping procedures in atmospheric modelling. In this application the condition number is large but independent of the grid resolution and both methods are asymptotically optimal, albeit with different absolute performance. In particular, an important constant in the discretisation is the CFL number; only the multigrid solver is robust to changes in this constant. We parallelise the solvers and adapt them to the specific features of GPU architectures, paying particular attention to efficient global memory access. We achieve a performance of up to 0.78 PFLOPs when solving an equation with 0.55 · 10 12 unknowns on 16384 GPUs; this corresponds to about 3% of the theoretical peak performance of the machine and we use more than 40% of the peak memory bandwidth with a Conjugate Gradient (CG) solver. Although the other solver, a geometric multigrid algorithm, has a slightly worse performance in terms of FLOPs per second, overall it is faster as it needs less iterations to converge; the multigrid algorithm can solve a linear PDE with half a trillion unknowns in about one second.

Computer Science↗

The Great Geospace Observatory and Simultaneous Missions of Opportunity

A predictive understanding of the sun to geospace environment is one of the main goals of ILWS. This can only be achieved through a "system-level" approach, meaning long-term, simultaneous, continuous observations across the relevant scales of the magnetosphere and ionosphere/thermosphere (IT). To date such an approach, which must involve simultaneous, multi-scale, global imaging of different geospace regions, has not been carried out for a complete geomagnetic storm. Such imagery, now routine for the Solar community, is of critical scientific importance and captures public imagination. Its absence in geospace studies has limited the growth and impact of geospace science. In this presentation, we discuss a concept called the Great Geospace Observatory, which would involve coordinated geospace imaging through an international effort of multiple, simultaneous Missions of Opportunity. In this way, the cost would be spread among different agencies as well as putting remote sensors in vantage points optimized for each type of imaging. 24/7 auroral imaging from weather satellites on Molniya (or similar) orbits, EUV imaging of the plasmasphere from high-inclination orbits, continuous and global ENA imaging from geosynchronous commercial satellites, and continuous X-ray imaging of the cusp and magnetosheath from a high-altitude dedicated probe would quantitatively track system-level dynamics at through substorms, sawtooth events, steady magnetospheric convection, and storms; studying energy and mass coupling between the solar wind, magnetosphere, and the upper atmosphere. In our minds, The Great Geospace Observatory represents the next strategic step for ILWS and needs to be seriously considered.

Donovan, Eric↗

The Great Geospace Observatory and Simultaneous Missions of Opportunity

A predictive understanding of the sun to geospace environment is one of the main goals of ILWS. This can only be achieved through a "system-level" approach, meaning long-term, simultaneous, continuous observations across the relevant scales of the magnetosphere and ionosphere/thermosphere (IT). To date such an approach, which must involve simultaneous, multi-scale, global imaging of different geospace regions, has not been carried out for a complete geomagnetic storm. Such imagery, now routine for the Solar community, is of critical scientific importance and captures public imagination. Its absence in geospace studies has limited the growth and impact of geospace science. In this presentation, we discuss a concept called the Great Geospace Observatory, which would involve coordinated geospace imaging through an international effort of multiple, simultaneous Missions of Opportunity. In this way, the cost would be spread among different agencies as well as putting remote sensors in vantage points optimized for each type of imaging. 24/7 auroral imaging from weather satellites on Molniya (or similar) orbits, EUV imaging of the plasmasphere from high-inclination orbits, continuous and global ENA imaging from geosynchronous commercial satellites, and continuous X-ray imaging of the cusp and magnetosheath from a high-altitude dedicated probe would quantitatively track system-level dynamics at through substorms, sawtooth events, steady magnetospheric convection, and storms; studying energy and mass coupling between the solar wind, magnetosphere, and the upper atmosphere. In our minds, The Great Geospace Observatory represents the next strategic step for ILWS and needs to be seriously considered.

Donovan, Eric↗

Mars Mission Optimization Based on Collocation of Resources

This paper presents a powerful approach for analyzing Martian data and for optimizing mission site selection based on resource collocation. This approach is implemented in a program called PROMT (Planetary Resource Optimization and Mapping Tool), which provides a wide range of analysis and display functions that can be applied to raw data or imagery. Thresholds, contours, custom algorithms, and graphical editing are some of the various methods that can be used to process data. Output maps can be created to identify surface regions on Mars that meet any specific criteria. The use of this tool for analyzing data, generating maps, and collocating features is demonstrated using data from the Mars Global Surveyor and the Odyssey spacecraft. The overall mission design objective is to maximize a combination of scientific return and self-sufficiency based on utilization of local materials. Landing site optimization involves maximizing accessibility to collocated science and resource features within a given mission radius. Mission types are categorized according to duration, energy resources, and in-situ resource utilization. Optimization results are shown for a number of mission scenarios.

G E Chamitoff↗

WE-H-207A-03: The Universality of the Lognormal Behavior of [F-18]FLT PET SUV Measurements

Purpose: Log transforming [F-18]FDG PET standardized uptake values (SUVs) has been shown to lead to normal SUV distributions, which allows utilization of powerful parametric statistical models. This study identified the optimal transformation leading to normally distributed [F-18]FLT PET SUVs from solid tumors and offers an example of how normal distributions permits analysis of non-independent/correlated measurements. Methods: Forty patients with various metastatic diseases underwent up to six FLT PET/CT scans during treatment. Tumors were identified by nuclear medicine physician and manually segmented. Average uptake was extracted for each patient giving a global SUVmean (gSUVmean) for each scan. The Shapiro-Wilk test was used to test distribution normality. One parameter Box-Cox transformations were applied to each of the six gSUVmean distributions and the optimal transformation was found by selecting the parameter that maximized the Shapiro-Wilk test statistic. The relationship between gSUVmean and a serum biomarker (VEGF) collected at imaging timepoints was determined using a linear mixed effects model (LMEM), which accounted for correlated/non-independent measurements from the same individual. Results: Untransformed gSUVmean distributions were found to be significantly non-normal (p<0.05). The optimal transformation parameter had a value of 0.3 (95%CI: −0.4 to 1.6). Given the optimal parameter was close to zero (which corresponds to log transformation), the data were subsequently log transformed. All log transformed gSUVmean distributions were normally distributed (p>0.10 for all timepoints). Log transformed data were incorporated into the LMEM. VEGF serum levels significantly correlated with gSUVmean (p<0.001), revealing log-linear relationship between SUVs and underlying biology. Conclusion: Failure to account for correlated/non-independent measurements can lead to invalid conclusions and motivated transformation to normally distributed SUVs. The log transformation was found to be close to optimal and sufficient for obtaining normally distributed FLT PET SUVs. These transformations allow utilization of powerful LMEMs when analyzing quantitative imaging metrics.

60 APPLIED LIFE SCIENCES↗

Bayesian Optimization of a Free-Electron Laser

The Linac Coherent Light Source X-ray free-electron laser is a complex scienti c apparatus which changes con gurations multiple times per day, necessitating fast tuning strategies to reduce setup time for successive experiments. To this end, we employ a Bayesian approach to maximizing X- ray laser pulse energy by controlling groups of quadrupole magnets. A Gaussian process model provides probabilistic predictions for the machine response with respect to control parameters, enabling a balance of exploration and exploitation in the search for the global optimum. We show that the model parameters can be learned from archived scans, and correlations between devices can be extracted from the beam transport. Furthermore, the result is a sample-e cient optimization routine, combining both historical data and knowledge of accelerator physics to signi cantly outperform existing optimizers.

47 OTHER INSTRUMENTATION↗

Optimized Extraction Method To Remove Humic Acid Interferences from Soil Samples Prior to Microbial Proteome Measurements

The microbial composition and their activities in soil environments play a critical role in organic matter transformation and nutrient cycling, perhaps most specifically with respect to impact on plant growth but also more broadly to global impact on carbon and nitrogen-cycling. Liquid chromatography coupled to high performance mass spectrometry provides a powerful approach to characterize soil microbiomes; however, the limited microbial biomass and the presence of abundant interferences in soil samples present major challenges to soil proteome extraction and subsequent MS measurement. To address some of the major issues, we have designed and optimized an experimental method to enhance microbial proteome extraction concomitant with minimizing the soil-borne humic substances co-extraction from soils. Among the range of interferences, humic substances are often the worst in terms of adversely impacting proteome extraction and mass spectrometry measurement. Our approach employs an in-situ detergent-based microbial lysis / TCA precipitation coupled with an additional acidification precipitation step at the peptide level which efficiently removes humic acids. By combing filtration and pH adjustment of the final peptide solution, the remaining humic acids can be differentially precipitated and removed with a membrane filter, thereby leaving much cleaner proteolytic peptide samples for MS measurement. As a result, this modified method is a reliable and straight-forward protein extraction method that efficiently removes soil-borne humic substances without inducing proteome sample loss or reducing or biasing protein identification in mass spectrometry.

59 BASIC BIOLOGICAL SCIENCES↗

10-th order of accuracy for numerical solution of 3-D elasticity equations for heterogeneous materials on unfitted Cartesian meshes

We have developed the Optimal Local Truncation Error Method (OLTEM) with 10-th order of accuracy on unfitted Cartesian meshes for a system of 3-D elasticity equations with smooth irregular interfaces. 5 x 5 x 5 = 125-point stencils (similar to those for quadratic finite elements) for elastic heterogeneous materials are used for OLTEM. There are no unknowns at the interface points between different materials; the structure of the global discrete equations is the same for homogeneous and heterogeneous materials. The calculation of unknown stencil coefficients is based on the minimization of the local truncation error of the stencil equations and yields the optimal 10-th order of accuracy for OLTEM on unfitted Cartesian meshes, i.e., the increase by 7 orders in accuracy compared to quadratic finite elements on conformal meshes. A new post-processing procedure provides the 9-th order of accuracy for stresses in the 3-D case. Similar to basic computations it uses OLTEM with the 125-point stencils, the interface conditions and the elasticity equations. It was shown that the use of the elasticity equations for post-processing improves the accuracy of 0.1% stresses by 6 orders compared to post-processing without the use of PDEs. At an accuracy of for stresses, OLTEM with the new post-processing procedure reduces the number of degrees of freedom by 360 - 8000 times compared to quadratic finite elements with similar stencils. OLTEM with the 125-point stencils yields even more accurate results than high-order finite elements with much wider stencils. OLTEM provides accurate numerical results for compressible and nearly incompressible materials.

elasticity equations↗

Differences Between OCO‐2 and GOME‐2 SIF Products From a Model‐Data Fusion Perspective

Space-borne retrievals of solar-induced chlorophyll fluorescence (SIF) over land surfaces have recently become a resource for studying and quantifying the broad scale dynamics of gross carbon uptake (gross primary productivity—GPP) across ecosystems. To prepare for the assimilation of SIF data in terrestrial biosphere models, we examine how differences between SIF products (due to differences in acquisition characteristics and processing chain) may affect the optimization of model parameters and the resultant GPP estimate. We compare recent daily mean SIF products (one from the Orbiting Carbon Observatory-2 [OCO-2] and two from the Global Ozone Monitoring Experiment–2 [GOME-2], GlobFluo [GF] and NASA-v28 [N28], missions), averaged at 0.5° × 0.5° spatial resolution and 16-day temporal resolution, at the biome level. Phase differences between these products are relatively small. A first-order correction of the difference in spectral sampling between the two instruments shows that OCO-2 and N28 are consistent in terms of magnitude and amplitude, while GF is twice as large as the others. Using a bias-blind toy data assimilation framework, we analyze how biases between SIF products, and between model and products, can be partially alleviated by optimizing the slope and intercept parameters of a linear GPP-SIF operator. As observation biases can transfer to biases in other optimized process-based parameters and to modeled carbon fluxes— thereby resulting in unidentified inaccurate parameter values—we argue that potential SIF biases should be treated cautiously in real-world experiments in order to achieve realistic and reliable future simulations.

Gross primary production↗

First Use of Synoptic Vector Magnetograms for Global Nonlinear, Force-Free Coronal Magnetic Field Models

Context. The magnetic field permeating the solar atmosphere is generally thought to provide the energy for much of the activity seen in the solar corona, such as flares, coronal mass ejections (CMEs), etc. To overcome the unavailability of coronal magnetic field measurements, photospheric magnetic field vector data can be used to reconstruct the coronal field. Currently, there are several modelling techniques being used to calculate three-dimensional field lines into the solar atmosphere. Aims. For the first time, synoptic maps of a photospheric-vector magnetic field synthesized from the vector spectromagnetograph (VSM) on Synoptic Optical Long-term Investigations of the Sun (SOLIS) are used to model the coronal magnetic field and estimate free magnetic energy in the global scale. The free energy (i.e., the energy in excess of the potential field energy) is one of the main indicators used in space weather forecasts to predict the eruptivity of active regions. Methods. We solve the nonlinear force-free field equations using an optimization principle in spherical geometry. The resulting threedimensional magnetic fields are used to estimate the magnetic free energy content E(sub free) = E(sub nlfff) − E(sub pot), which is the difference of the magnetic energies between the nonpotential field and the potential field in the global solar corona. For comparison, we overlay the extrapolated magnetic field lines with the extreme ultraviolet (EUV) observations by the atmospheric imaging assembly (AIA) on board the Solar Dynamics Observatory (SDO). Results. For a single Carrington rotation 2121, we find that the global nonlinear force-free field (NLFFF) magnetic energy density is 10.3% higher than the potential one. Most of this free energy is located in active regions.

Sun: corona↗