Search NASA⌕ Search

SEARCH · Search NASA

Results for “iterative analysis”

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 307 records · Page 17

Anderson acceleration stability in NDA-accelerated k-eigenvalue problems

Anderson acceleration (AA) has been used to improve the stability and convergence rate of multiphysics iterative methods for reactor analysis. Most applications studied assume a tightly converged solution for the different physics problems, and AA is usually applied to state variables like temperature, density, and heat generation rate. In this paper, we study the theoretical performance of AA in NDA-accelerated k-eigenvalue problems. The problems and algorithms studied are simplified from the coupled iteration scheme adopted by MPACT and many other high-fidelity whole-core reactor codes. Compared to previous analyses of AA for these iteration schemes, we study the case with a partially converged neutronics solution and possibly partially converged nonlinear diffusion acceleration (NDA)/coarse mesh finite difference (CMFD) solutions. We observe that the performance of the iteration scheme with AA is very sensitive to the initial guess and is affected by the partially converged CMFD solutions. When the NDA solution is fully converged, using AA cannot achieve the optimal convergence rate in large-sized problems. Conversely, if the NDA solution is partially converged, the iteration scheme with AA can diverge or converge extremely slowly. It is found that the loss of robustness for AA is due to the fact that it is applied to the iterative subspace of state variables rather than the fundamental unknowns of the governing equations. To improve the robustness, the scalar flux should also be considered in the implementation of AA. After considering the residuals of flux, we observe that the stability is regardless of the partial convergence of NDA solutions. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Integrated numerical methods for hypersonic aircraft cooling systems analysis

Numerical methods have been developed for the analysis of hypersonic aircraft cooling systems. A general purpose finite difference thermal analysis code is used to determine areas which must be cooled. Complex cooling networks of series and parallel flow can be analyzed using a finite difference computer program. Both internal fluid flow and heat transfer are analyzed, because increased heat flow causes a decrease in the flow of the coolant. The steady state solution is a successive point iterative method. The transient analysis uses implicit forward-backward differencing. Several examples of the use of the program in studies of hypersonic aircraft and rockets are provided.

Petley, Dennis H.↗

Image segmentation by iterative parallel region growing with application to data compression and image analysis

Image segmentation can be a key step in data compression and image analysis. However, the segmentation results produced by most previous approaches to region growing are suspect because they depend on the order in which portions of the image are processed. An iterative parallel segmentation algorithm avoids this problem by performing globally best merges first. Such a segmentation approach, and two implementations of the approach on NASA's Massively Parallel Processor (MPP) are described. Application of the segmentation approach to data compression and image analysis is then described, and results of such application are given for a LANDSAT Thematic Mapper image.

Tilton, James C.↗

Comprehensive analysis of disruption mitigation methods using gas and pellet-like injections in ITER-like Tokamaks

Abstract Inert-gas shielding could be an effective mechanism for protection of plasma facing surfaces (PFS) against plasma particles impact and photon radiation heat loads during transient events in fusion devices. Neutral gas injection is one promising way to mitigate erosion of tokamak components and contamination. The objective of this work is to study and optimize mitigation methods using neutral gas and pellet-like injections to decrease the heat load to the divertor surfaces and to prevent vaporization of the various internal surfaces due to transient events in ITER-like devices. The integrated self-consistent models implemented in the HEIGHTS package was used for detailed analysis of the potential secondary plasma generation from the injected inert gas, its radiative characteristics, and shielding effectiveness. We varied the density, size, and location of an argon gas cloud to minimize the disruption energy deposited into the divertor components. We also investigated innovative ways for minor changes in ITER-like internal design to mitigate disruptions. We found the optimum parameters to fully protect ITER tokamak surfaces from erosion and vaporization during plasma instabilities. This preliminary analysis showed that using Ar gas injection methods could lead to enhancement in components lifetime in ITER-like and future DEMO devices with minor design changes.

Science & Technology - Other Topics↗

Valuation-Based Framework for Considering Distributed Generation Photovoltaic Tariff Design: Preprint

While an export tariff is only one element of a larger regulatory framework for distributed generation, we choose to focus on tariff design because of the significant impact this program design component has on the various flows of value among power sector stakeholders. In that context, this paper is organized into a series of steps that can be taken during the design of a DGPV export tariff design. To that end this paper outlines a holistic, high-level approach to the complex undertaking of DGPV tariff design, the crux of which is an iterative cost-benefit analysis process. We propose a multi-step progression that aims to promote transparent, focused, and informed dialogue on CBA study methodologies and assumptions. When studies are completed, the long-run marginal avoided cost of the DGPV program should be compared against the costs imposed on utilities and non-participating customers, recognizing that these can be defined differently depending on program objectives. The results of this comparison can then be weighed against other program objectives to formulate tariff options. Potential changes to tariff structures can be iteratively fed back into established analytical tools to inform further discussions.

14 SOLAR ENERGY↗

Regression Analysis and Calibration Recommendations for the Characterization of Balance Temperature Effects

Analysis and use of temperature-dependent wind tunnel strain-gage balance calibration data are discussed in the paper. First, three different methods are presented and compared that may be used to process temperature-dependent strain-gage balance data. The first method uses an extended set of independent variables in order to process the data and predict balance loads. The second method applies an extended load iteration equation during the analysis of balance calibration data. The third method uses temperature-dependent sensitivities for the data analysis. Physical interpretations of the most important temperature-dependent regression model terms are provided that relate temperature compensation imperfections and the temperature-dependent nature of the gage factor to sets of regression model terms. Finally, balance calibration recommendations are listed so that temperature-dependent calibration data can be obtained and successfully processed using the reviewed analysis methods.

wind tunnel test↗

RSRM Chamber Pressure Oscillations: Transit Time Models and Unsteady CFD

Space Shuttle solid rocket motor low frequency internal pressure oscillations have been observed since early testing. The same type of oscillations also are present in the redesigned solid rocket motor (RSRM). The oscillations, which occur during RSRM burn, are predominantly at the first three motor cavity longitudinal acoustic mode frequencies. Broadband flow and combustion noise provide the energy to excite these modes at low levels throughout motor burn, however, at certain times during burn the fluctuating pressure amplitude increases significantly. The increased fluctuations at these times suggests an additional excitation mechanism. The RSRM has inhibitors on the propellant forward facing surface of each motor segment. The inhibitors are in a slot at the segment field joints to prevent burning at that surface. The aft facing segment surface at a field joint slot burns and forms a cavity of time varying size. Initially the inhibitor is recessed in the field joint cavity. As propellant burns away the inhibitor begins to protrude into the bore flow. Two mechanisms (transit time models) that are considered potential pressure oscillation excitations are cavity-edge tones, and inhibitor hole-tones. Estimates of frequency variation with time of longitudinal acoustic modes, cavity edge-tones, and hole-tones compare favorably with frequencies measured during motor hot firing. It is believed that the highest oscillation amplitudes occur when vortex shedding frequencies coincide with motor longitudinal acoustic modes. A time accurate computational fluid dynamic (CFD) analysis was made to replicate the observations from motor firings and to observe the transit time mechanisms in detail. FDNS is the flow solver used to detail the time varying aspects of the flow. The fluid is approximated as a single-phase ideal gas. The CFD model was an axisymmetric representation of the RSRM at 80 seconds into burn.Deformation of the inhibitors by the internal flow was determined through an iterative structural and CFD analysis. The analysis domain ended just upstream of the nozzle throat. This is an acoustic boundary condition that caused the motor to behave as a closed-open organ pipe. This differs from the RSRM which behaves like a closed-closed organ pipe. The unsteady CFD solution shows RSRM chamber pressure oscillations predominately at the longitudinal acoustic mode frequencies of a closed-open organ pipe. Vortex shedding in the joint cavities and at the inhibitors contribute disturbances to the flow at the second longitudinal acoustic mode frequency. Further studies are planned using an analysis domain that extends downstream of the nozzle throat.

Nesman, Tom↗

TRANSP integrated modeling code for interpretive and predictive analysis of tokamak plasmas

This paper provides a comprehensive review of the TRANSP code, a sophisticated tool for interpretive and predictive analysis of tokamak plasmas, detailing its major capabilities and features. It describes the equations for particle, power, and momentum balance analysis, as well as the poloidal field diffusion equations. The paper outlines the spatial and time grids used in TRANSP and details the equilibrium assumptions and solvers. Various models for heating and current drive and radiation, including updates to the NUBEAM model, are discussed. The handling of large-scale events such as sawtooth crashes and pellet injections is examined, along with the predictive capabilities for advancing plasma profiles. The integration of TRANSP with the ITER Integrated Modeling and Analysis Suite (IMAS) is highlighted, demonstrating enhanced data access and analysis capabilities. Additionally, the paper discusses best practices and continuous integration techniques to enhance TRANSP's robustness. The suite of TRANSP tools, designed for efficient data analysis and simulation, further supports the optimization of tokamak operations and coupling with other tokamak codes. Continuous development and support ensure that TRANSP remains a major code for the analysis of experimental data for controlled thermonuclear fusion, maintaining its critical role in supporting the optimization of tokamak operations and advancing fusion research.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Description of the prototype diagnostic residual gas analyzer for ITER

The diagnostic residual gas analyzer (DRGA) system to be used during ITER tokamak operation is being designed at Oak Ridge National Laboratory to measure fuel ratios (deuterium and tritium), fusion ash (helium), and impurities in the plasma. The eventual purpose of this instrument is for machine protection, basic control, and physics on ITER. Prototyping is ongoing to optimize the hardware setup and measurement capabilities. The DRGA prototype is comprised of a vacuum system and measurement technologies that will overlap to meet ITER measurement requirements. Three technologies included in this diagnostic are a quadrupole mass spectrometer, an ion trap mass spectrometer, and an optical penning gauge that are designed to document relative and absolute gas concentrations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Microchannel-based Membrane-less Extraction of Li from Unconventional Lithium Sources & the Separation of REE

This final report provides an overview of the Project's entire duration, covering July 1, 2021 to December 31, 2023. It primarily focuses on the achievements, technological developments, and unique challenges the team faced while working on separating and extracting Lithium from produced waters. The project's primary aim was to create an integrated, high-throughput, membrane-less, and modular microfluidic platform that could extract Lithium from unconventional sources. We have successfully met all goals and milestones envisioned in the SOPO document. The most critical primary milestones, including the Go-No-Go milestone (refer to the Gantt chart in the Appendices), were successfully accomplished. We demonstrated phase separation (>90%) and extraction (>85%) performance in the MPSE using synthetic, and representative produced water composition feed at 50 ml/min total flow through MPSE 36. We have also performed a parametric study of the MPSE operations, beyond the scope of SOPO, exploring operating conditions of current and broader interest. The extended investigation of operational parameters is concurrent with our efforts to seek further development of the MPSE technology beyond the scope of the Project. Along these lines of development, we have made efforts to be responsive to DOE calls for technological developments of other types of resources (beyond PW) for the recovery of Critical Materials and higher TRL development (beyond TRL 4). During the work on this Project, we developed and implemented three innovative technical approaches that emerged from our efforts to successfully meet the Project milestones. The innovative & original technical approaches developed and implemented in this Project are now the contributions to process engineering that could be clearly credited to the Project. First, Convergent Design Approach is a comprehensive feedforward & feedback loop of four design phases: i) design for functionality, ii) design for manufacturing, iii) design for sustainability, and iv) design for market. Next was Process Intensification. A major aim of this Project was to create an innovative phase separation & extraction microscale-based technology for Li separation – thus the words microchannel-based in the Project title. A microscale-based technology is intrinsically in the center of the Process Intensification domain as defined by its unique principles. Therefore, Process Intensification was implicitly envisioned in the Project’s SOPO. Lastly, Time Scale Analysis is a novel tool for discovering the needs and directions of Process Intensification implementations in any process technology. This Project is fully credited for developing and implementing the three novel technical approaches mentioned above. These are general contributions to process engineering that emerged from this Project. Beyond the original SOPO scope, the OSU-U.Pitt research group utilized a Convergent Design methodology, integrating first-principles mathematical modeling with experimental validation on the Minimum Development Vehicle. By creating these Digital Twins, the team rapidly assessed manufacturing iterations to support TEA analysis. This framework further enabled the development of advanced Surface Modification Techniques, where hydrophobic and oleophobic coating strategies were optimized via Digital Twin tools and validated through rigorous 100-hour longevity testing. TEA Analysis: The closing efforts of this Project were focused on the TEA analysis. TEA analysis had two primary functions: i) enabling critical assessments of design variations withing 10 the Concurrent Design Approach, thus enabling evolution of the MPSE design to reach faster- better-cheaper alternatives; and ii) to create a bridge between the accomplishments of this Project and future projects of higher TRL, beyond TRL 6 level. It is important to note that the TEA model created in the Project stirred the technological solutions for the recovery of critical materials toward a vision of a very profitable modular plant that has unique zero-waste water discharge signature. More importantly, thanks to our experimental performance data and conservative assumptions, the TEA model predicts minimal technological and investment risks. Low cost of a modular unit of a nominal capacity of [1000 tons of Li 2 CO 3 /year] positions the MPSE based technology within the reach of community investors, thus offering a paradigm shift in the development of critical technologies. The project successfully navigated two primary challenges: solvent selection and manufacturing adaptation. Restricted by the SOPO to existing literature for lithium recovery, the team identified a critical need for a "material excellence program" to develop next-generation solvents, eventually concluding with a preliminary investigation into promising Ionic Liquids (ILs). Simultaneously, COVID-19 supply chain disruptions forced a pivot from traditional manufacturing to advanced additive methods at ATAMI-OSU. By transitioning from stainless steel to 3D-printed polymer substrates, the team achieved a transformative three-order-of- magnitude reduction in manufacturing costs and compressed prototyping timelines from several months to just two days. The MPSE technology offers significant energy, environmental, and economic advantages by overcoming the traditional bottlenecks of phase-separation hardware and contactor size. Unlike conventional mixer-settlers or membrane-based systems, MPSE operates without moving parts or fouling-prone membranes, achieving robust performance even with challenging, viscous, or particulate-heavy feeds. Key performance metrics include an energy intensity reduction of 5–50x (3–40 kJ/m 3 ) compared to incumbent technologies and a dramatic reduction of processing time to under 60 seconds, which drastically reduces the physical plant footprint. These technical efficiencies translate into superior economic outcomes; for a 100 t/year Li 2 CO 3 facility, implementing MPSE is projected to nearly halve contactor CAPEX (from $\$$6.08M to $\$$3.01M) and significantly increase the project's Net Present Value (NPV), derisking new investment and enabling distributed critical-mineral processing configurations. The commercialization of MPSE technology is being spearheaded by Vigsur Dynamics Inc., which has adopted a structured, parallel approach to technical and business development since its formation in January 2026. Following extensive customer discovery and engagement with the Oregon State University accelerator, Vigsur Dynamics is working to establish a business model that transitions from pilot demonstrations to modular hardware sales, ultimately aiming for a "build-own-operate" service strategy. Current technical milestones—including 100 hours of continuous operation, superior energy efficiency, and successful 6-unit modular scale-up— provide a foundation for this transition. Backed by ongoing IP licensing and a growing network of industrial and venture advisors, the company is actively de-risking the platform to replace conventional mixer-settler systems in the critical minerals market.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Statistical analysis of the Nb 3 Sn strand production for the ITER toroidal field coils

The ITER toroidal field (TF) strand procurement initiated the largest Nb 3 Sn superconducting strand production hitherto. The industrial-scale production started in Japan in 2008 and finished in summer 2015. Six ITER partners (so-called Domestic Agencies, or DAs) are in charge of the procurement and involved eight different strand suppliers all over the world, of which four are using the bronze route (BR) process and four the internal-tin (IT) process. In total more than 500 tons have been produced including excess material covering losses during the conductor manufacturing process, in particular the cabling. The procurement is based on a functional specification where the main strand requirements like critical current, hysteresis losses, Cu ratio and residual resistance ratio are specified but not the strand production process or layout. We introduce the analysis on the data acquired during the quality control (QC) process that was carried out to ensure the same conductor performance requirements are met by the different strand suppliers regardless of strand design. The strand QC is based on 100% billet testing and on applying statistical process control (SPC) limits. Throughout the production, samples adjacent to the strand pieces tested by the suppliers are cross-checked ('verified') by their respective DAs reference labs. The level of verification was lowered from 100% at the beginning of the procurement progressively to approximately 25% during the final phase of production. Based on the complete dataset of the TF strand production, an analysis of the SPC limits of the critical strand parameters is made and the related process capability indices are calculated. In view of the large-scale production and costs, key manufacturing parameters such as billet yield, number of breakages and piece-length distribution are also discussed. Furthermore, the results are compared among all the strand suppliers, focusing on the difference between BR and IT processes. Following the completion of the largest Nb 3 Sn strand production, our experience gained from monitoring the execution of the QC activities and from auditing the results from the measurements is summarised for future superconducting strand material procurement activities.

Vostner, A.↗

Graphical and Numerical Description of Strain-Gage Balance Interactions

A new approach for the graphical and numerical description of balance interactions is presented. The approach uses data from single-component loads as input. This choice has the two advantages. First, the number of applied loads is at the minimum needed for interactions to be observed. In addition, the loads used for the description of interactions can easily be repeated at different sites. Output differences relative to the outputs of the zero load point of a load series are used for the description of interactions. Similarly, load differences relative to the loads of the zero load point of a load series are used for the description of the loads. Interactions are plotted versus the load differences for each load component while omitting the outputs of the primary gage of the chosen load component. The resulting plots have a characteristic star pattern as interactions are zero at zero load. Numerical estimates of the slopes of the interactions can be reverse-engineered from the load prediction equations of a balance if calibration data is examined. The slopes are the off-diagonal coefficients of the inverse of the matrix that has the coefficients of the linear terms of the fitted loads if the Non-Iterative Method is used for the analysis. Similarly, the slopes are the off-diagonal coefficients of the inverse of the matrix that is the non-iterative part of the primary load iteration equation if the Iterative Method is used for the analysis. Data from the calibration of a force balance is processed to illustrate the proposed graphical and numerical description of interactions.

wind tunnel test↗

Graphical and Numerical Description of Strain-Gage Balance Interactions

A new approach for the graphical and numerical description of strain-gage balance interactions is presented. The approach uses data from single-component loads as input. This choice has the two advantages. First, the number of applied loads is at the minimum needed for interactions to be observed. In addition, the loads used for the description of interactions can easily be repeated at different sites. Output differences relative to the outputs of the zero load point of a load series are used for the description of interactions. Similarly, load differences relative to the loads of the zero load point of a load series are used for the description of the loads. Interactions are plotted versus the load differences for each load component while omitting outputs of the primary gage of the chosen load component. The resulting plots have a star pattern as all interactions are zero at zero load. Estimates of the slopes of the interactions can be obtained from the load prediction equations of a balance if calibration data is examined. The slopes are the off-diagonal coefficients of the inverse of the matrix that has the coefficients of the linear terms of the fitted loads if the Non-Iterative Method is used for the analysis. Similarly, the slopes are the off-diagonal coefficients of the inverse of the matrix that is the non-iterative part of the primary load iteration equation if the Iterative Method is used for the analysis. Data sets from a manual and a machine calibration of a force balance are processed to illustrate the proposed description of interactions.

strain-gage balance↗

Elevated temperature mechanical properties of Inconel 617 surface oxide using nanoindentation

Inconel 617 is a principal candidate material for helium gas cooled very-high-temperature reactors with outlet temperatures of 700–950 °C. Recent findings showed that this alloy develops unique surface oxide layers especially at high temperature (HT) helium environment with distinctive wear, friction and contact properties. This study investigates the elevated temperature mechanical properties of Inconel 617 top surface layers aged in HT helium environment. Nanoindentation technique is used to obtain load-displacement graphs of the alloy top surface oxide in temperatures ranging from 25 °C up to 600 °C. In addition, using finite element analysis along with an iterative regression technique, a semi-numerical method is developed to further measure and quantify the material parameters and, in particular, time-independent and creep characteristics of the oxide. Although Young's modulus of the oxide is found to be relatively close to the bulk for the tested temperatures, the yield strength and hardness, in comparison to the bulk material, increase significantly as the material is oxidized after aging. The oxide exhibits significant softening as the temperature increases to 600 °C. Unlike the bulk material, diffusion through the grains is found to be the dominant creep mechanism for the oxide. Considerable difference between the mechanical properties of the oxide and the bulk material shows the need for accurate measurements of near surface mechanical properties, if reliable predictive contact and tribological models are sought at elevated temperatures.

36 MATERIALS SCIENCE↗

Experimental study of ECH pre-ionization on J-TEXT

An experimental study on electron cyclotron heating (ECH) pre-ionization has been conducted on J-TEXT in support of the joint experiment research for ITER plasma initiation. In this experiment, ECH power was injected to the vessel before the application of loop voltage ionize the neutral gas and form the initial plasma or so-called pre-plasma. The impact of several significant factors, such as magnetic field configuration, pre-fill gas pressure, ECH toroidal injection angle and ECH power on the evolution of pre-plasma are systematically studied, aiming to identify shared features, clarify their potential relationship and optimize the discharge parameters to generate a rather high pre-plasma density. By separating ECH power from the inductive start-up, the effect of pre-plasma on tokamak start-up can be observed. A dynamic magnetic configuration facilitates the transition of pre-plasma to tokamak plasma. To assess the influence of pre-plasma density on tokamak start-up, two kinds of magnetic field configurations are examined. While the effect of different pre-plasma densities on tokamak start-up is negligible, a significant difference is observed between pure ohmic start-up and start-up with pre-ionization. The studies presented here show evolutionary trends and threshold values needed to optimize ECH pre-ionization and also a feasible way to improve pre-plasma density and a configuration to stabilize pre-plasma for transition. Eventually, these results may contribute to multi-machine research and physics analysis that can assist ITER with its optimal preparation for first plasma operation.

ECH↗

Development of a Fast Spectrum Test Loop at the Advanced Test Reactor

The United States currently has no infrastructure to provide a fast neutron flux for research purposes. With the increased interest in Sodium Fast Reactor (SFR) technology, this type of infrastructure is necessary for further development of key components such as Versatile Test Reactor. The purpose of this research is to characterize the thermal hydraulics of a test loop capable of providing SFR like operating conditions. This includes showcasing feasibility in mechanical design, reaching a high fast to thermal neutron flux ratio ideally near 15, an inlet temperature of 600 K, a temperature increase of 50 K from inlet to outlet, and test chamber radial temperature profiles similar to SFR fuel rods. Existing infrastructure must be used to provide the neutron flux for the test loop, and it was decided that the Advanced Test Reactor (ATR) was the optimal choice due to its availability and previous work conducted similar to this project. To achieve the goals of this project, several mechanical designs were put forward with the most feasible design being chosen. The neutronics and thermal hydraulics analysis were conducted in tandem, with an iterative process used to maximize efficiency of both systems. Analysis of the neutronics system was done using Attila and MCNP while the thermal hydraulic analysis was done using RELAP5. This document outlines the thermal hydraulic analysis, though it includes the results from the mechanical design and neutronics analysis. The loop design includes a primary system with molten sodium, a set of “booster fuel plates” to increase the fast flux within the test section, a water system to cool the booster fuel plates, a helium system to separate the sodium and water systems, and a thermal neutron shield to decrease the thermal flux. Each of these systems is modeled in the neutronics and thermal hydraulic analyses. The RELAP5 model showcased the feasibility of the design proposed by this report. For all designs tested, the desired 50 K temperature increase from inlet to outlet was reached. The fuel rod radial temperature profile was lower than was expected, likely due to a lower rod thermal power than expected. A parametric study was also conducted to characterize the system with three fuel pins. This study showed the effect of thermal power, fluid pressure, fluid velocity, and inlet temperature on the system. ©Copyright by Anton Higgins December 11, 2019 All Rights Reserved

99 GENERAL AND MISCELLANEOUS↗

Computational methods of robust controller design for aerodynamic flutter suppression

The development of Riccati iteration, a tool for the design and analysis of linear control systems is examined. First, Riccati iteration is applied to the problem of pole placement and order reduction in two-time scale control systems. Order reduction, yielding a good approximation to the original system, is demonstrated using a 16th order linear model of a turbofan engine. Next, a numerical method for solving the Riccati equation is presented and demonstrated for a set of eighth order random examples. A literature review of robust controller design methods follows which includes a number of methods for reducing the trajectory and performance index sensitivity in linear regulators. Lastly, robust controller design for large parameter variations is discussed.

Anderson, L. R.↗

ESTIMATING THE DEEP SOLAR MERIDIONAL CIRCULATION USING MAGNETIC OBSERVATIONS AND A DYNAMO MODEL: A VARIATIONAL APPROACH

We show how magnetic observations of the Sun can be used in conjunction with an axisymmetric flux-transport solar dynamo model in order to estimate the large-scale meridional circulation throughout the convection zone. Our innovative approach rests on variational data assimilation, whereby the distance between predictions and observations (measured by an objective function) is iteratively minimized by means of an optimization algorithm seeking the meridional flow that best accounts for the data. The minimization is performed using a quasi-Newton technique, which requires knowledge of the sensitivity of the objective function to the meridional flow. That sensitivity is efficiently computed via the integration of the adjoint flux-transport dynamo model. Closed-loop (also known as twin) experiments using synthetic data demonstrate the validity and accuracy of this technique for a variety of meridional flow configurations, ranging from unicellular and equatorially symmetric to multicellular and equatorially asymmetric. In this well-controlled synthetic context, we perform a systematic study of the behavior of our variational approach under different observational configurations by varying their spatial density, temporal density, and noise level, as well as the width of the assimilation window. We find that the method is remarkably robust, leading in most cases to a recovery of the true meridional flow to within better than 1%. These encouraging results are a first step toward using this technique to (i) better constrain the physical processes occurring inside the Sun and (ii) better predict solar activity on decadal timescales.

79 ASTRONOMY AND ASTROPHYSICS↗