Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic boundaries”

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,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Large eddy simulation of a boundary layer with concave streamwise curvature

One of the most exciting recent developments in the field of large eddy simulation (LES) is the dynamic subgrid-scale model. The dynamic model concept is a general procedure for evaluating model constants by sampling a band of the smallest scales actually resolved in the simulation. To date, the procedure has been used primarily in conjunction with the Smagorinsky model. The dynamic procedure has the advantage that the value of the model constant need not be specified a priori, but rather is calculated as a function of space and time as the simulation progresses. This feature makes the dynamic model especially attractive for flows in complex geometries where it is difficult or impossible to calibrate model constants. The dynamic model was highly successful in benchmark tests involving homogeneous and channel flows. Having demonstrated the potential of the dynamic model in these simple flows, the overall direction of the LES effort at CTR shifted toward an evaluation of the model in more complex situations. The current test cases are basic engineering-type flows for which Reynolds averaged approaches were unable to model the turbulence to within engineering accuracy. Flows currently under investigation include a backward-facing step, wake behind a circular cylinder, airfoil at high angles of attack, separated flow in a diffuser, and boundary layer over a concave surface. Preliminary results from the backward-facing step and cylinder wake simulations are encouraging. Progress on the LES of a boundary layer on a concave surface is discussed. Although the geometry of a concave wall is not very complex, the boundary layer that develops on its surface is difficult to model due to the presence of streamwise Taylor-Gortler vortices. These vortices arise as a result of a centrifugal instability associated with the convex curvature.

Lund, Thomas S.↗
NASA NTRS2022

A Multimodel Investigation of Asian Summer Monsoon UTLS Transport Over the Western Pacific

The Asian summer monsoon (ASM) as a chemical transport system is investigated using a suite of models in preparation for an airborne field campaign over the Western Pacific. Results show that the dynamical process of anticyclone eddy shedding in the upper troposphere rapidly transports convectively uplifted Asian boundary layer air masses to the upper troposphere and lower stratosphere over the Western Pacific. The models show that the transported air masses contain significantly enhanced aerosol loading and a complex chemical mixture of trace gases that are relevant to ozone chemistry. The chemical forecast models consistently predict the occurrence of the shedding events, but the predicted concentrations of transported trace gases and aerosols often differ between models. The airborne measurements to be obtained in the field campaign are expected to help reduce the model uncertainties. Furthermore, the large-scale seasonal chemical structure of the monsoon system is obtained from modeled carbon monoxide, a tracer of the convective transport of pollutants, which provides a new perspective of the ASM circulation, complementing the dynamical characterization of the monsoon.

Laura L. Pan↗
NASA NTRS1975

Thermospheric composition

The global view of the atmosphere provided by satellite observations and radar backscatter measurements has shown a pattern of composition variability resulting largely from dynamic processes. Important progress has been made in defining composition variability at 120 km, which has greatly improved models which previously had fixed boundary conditions at that altitude. The discovery of the early morning maximum in helium abundance and the better definition of diurnal variation in relative composition have provided a basis for more sophisticated dynamical models. Similarly, models of meridional wind systems have been developed which can account for the better-defined winter enhancements in helium and atomic oxygen. Finally, composition variability has indicated the presence of a high-latitude heat source which during periods of magnetic activity profoundly modifies the state of the atmosphere.

Carignan, G. R.↗
NASA NTRS2018

Simulating the Effects of Carbon Nanotube Continuity and Interfacial Bonding on Composite Strength and Stiffness

Molecular dynamics simulations of carbon nanotube (CNT) composites, in which the CNTs are continuous across the periodic boundary, overestimate the experimentally measured mechanical properties of CNT composites along the fiber direction. Since the CNTs in these composites are much shorter than the composite dimensions, load must be transferred either directly between CNTs or through the matrix, a mechanism that is absent in simulations of effectively continuous CNTs. In this study, the elastic and fracture properties of high volume fraction discontinuous carbon nanotube/amorphous carbon composite systems were compared to those of otherwise equivalent continuous CNT composites using ReaxFF reactive molecular dynamics simulations. These simulations were used to show how the number of nanotube-matrix interfacial covalent bonds affect composite mechanical properties. Furthermore, the mechanical impact of interfacial bonding was decomposed to reveal its effect on the properties of the CNTs, the interfacial layer of matrix, and the bulk matrix. For the composites with continuous reinforcement, it was found that any degree of interfacial bonding has a negative impact on axial tensile strength and stiffness. This is due to disruption of the structure of the CNTs and interfacial matrix layer by the interfacial bonds. For the discontinuous composites, the modulus was maximized between 4%-7% interfacial bonding and the strength continues to increase up to the highest levels of interfacial bonding studied. Areas of low stress and voids were observed in the simulated discontinuous composites at the ends of the tubes, from which fracture was observed to initiate. Experimental carbon nanotube yarn composites were fabricated and tested. The results were used to illustrate knockdown factors relative to the mechanical performance of the tubes themselves.

Benjamin D Jensen↗
NASA NTRS1990

Nonlinear dynamics of global atmospheric and Earth-system processes

Researchers are continuing their studies of the nonlinear dynamics of global weather systems. Sensitivity analyses of large-scale dynamical models of the atmosphere (i.e., general circulation models i.e., GCM's) were performed to establish the role of satellite-signatures of soil moisture, sea surface temperature, snow cover, and sea ice as crucial boundary conditions determining global weather variability. To complete their study of the bimodality of the planetary wave states, they are using the dynamical systems approach to construct a low-order theoretical explanation of this phenomenon. This work should have important implications for extended range forecasting of low-frequency oscillations, elucidating the mechanisms for the transitions between the two wave modes. Researchers are using the methods of jump analysis and attractor dimension analysis to examine the long-term satellite records of significant variables (e.g., long wave radiation, and cloud amount), to explore the nature of mode transitions in the atmosphere, and to determine the minimum number of equations needed to describe the main weather variations with a low-order dynamical system. Where feasible they will continue to explore the applicability of the methods of complex dynamical systems analysis to the study of the global earth-system from an integrative viewpoint involving the roles of geochemical cycling and the interactive behavior of the atmosphere, hydrosphere, and biosphere.

Saltzman, Barry↗
NASA NTRS1988

Aerodynamics of vortex generators

An experimental and theoretical study was undertaken of the separation delay and dramatic boundary-layer thinning that can occur in vortex-generator installations. Wind tunnel measurements of the dynamic-pressure profile downstream of a vortex generator were found to compare under certain conditions with that downstream of a suction slit, while water-tunnel visualization studies of vortex-generator height and geometry suggested optimum configurations, and only a minor effect of base porosity. A series of progressively more complex inviscid flow models was developed to be applied to a 3-D integral boundary-layer code. This code predicted layer thinning downstream of the suction site of the vortex models, and other observed features. Thin-layer Navier-Stokes equations are now being used with the ultimate goal of clarifying the physical processes involved in vortex generator performance and developing calculational procedures capable of predicting it.

Breidenthal, Robert E., Jr.↗
NASA NTRS1987

Experimental studies on the tripping behavior of narrow T-stiffened flat plates subjected to hydrostatic pressure and underwater shock

An experimental investigation was conducted to determine the static and dynamic responses of a specific stiffened flat plate design. The air-backed rectangular flat plates of 6061-T6 aluminum with an externally machined longitudinal narrow-flanged T-stiffener and clamped boundary conditions were subjected to static loading by water hydropump pressure and shock loading from an eight pound TNT charge detonated underwater. The dynamic test plate was instrumented to measure transient strains and free field pressure. The static test plate was instrumented to measure transient strains, plate deflection, and pressure. Emphasis was placed upon forcing static and dynamic stiffener tripping, obtaining relevant strain and pressure data, and studying the associated plate-stiffener behavior.

Budweg, H. L.↗
DOE OSTI2020

The mechanisms and seasonal differences of the impact of aerosols on daytime surface urban heat island effect

The urban heat island intensity (UHII) is the temperature difference between urban areas and their rural surroundings. It is commonly attributed to changes in the underlying surface structure caused by urbanization. Air pollution caused by aerosol particles can affect the UHII through changing (1) the surface energy balance by the aerosol radiative effect (ARE) and (2) planetary-boundary-layer (PBL) stability and airflow intensity by modifying thermodynamic structure, which is referred to as the aerosol dynamic effect (ADE). By analyzing satellite data and ground-based observations collected from 2001 to 2010 at 35 cities in China and using the WRF-Chem model, we find that the impact of aerosols on UHII differs considerably: reducing the UHII in summer but increasing the UHII in winter. This seasonal contrast is proposed to be caused by the different strengths of the ARE and ADE between summer and winter. In summer, the ARE on UHII is dominant over the ADE, cooling down surface temperature more strongly in urban areas than in rural areas because of much higher aerosol loading, and offsets the urban heating, therefore weakening UHII. In winter, however, the ADE is more dominant, because aerosols stabilize the PBL more in the polluted condition, weakening the near-surface heat transport over urban areas in both vertical and horizontal directions. This means that the heat accumulated in urban areas is dispersed less effectively, and thus the UHII is enhanced. These findings shed new light on the impact of the interaction between urbanization-induced surface changes and air pollution on urban climate.urban areas and their rural surroundings. It is commonly attributed to changes in the underlying surface structure caused by urbanization. Air pollution caused by aerosol particles can affect the UHII through changing (1) the surface energy balance by the aerosol radiative effect (ARE) and (2) planetary-boundary-layer (PBL) stability and airflow intensity by modifying thermodynamic structure, which is referred to as the aerosol dynamic effect (ADE). By analyzing satellite data and ground-based observations collected from 2001 to 2010 at 35 cities in China and using the WRF-Chem model, we find that the impact of aerosols on UHII differs considerably: reducing the UHII in summer but increasing the UHII in winter. This seasonal contrast is proposed to be caused by the different strengths of the ARE and ADE between summer and winter. In summer, the ARE on UHII is dominant over the ADE, cooling down surface temperature more strongly in urban areas than in rural areas because of much higher aerosol loading, and offsets the urban heating, therefore weakening UHII. In winter, however, the ADE is more dominant, because aerosols stabilize the PBL more in the polluted condition, weakening the near-surface heat transport over urban areas in both vertical and horizontal directions. This means that the heat accumulated in urban areas is dispersed less effectively, and thus the UHII is enhanced. These findings shed new light on the impact of the interaction between urbanization-induced surface changes and air pollution on urban climate.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2020

X-56A Structural Dynamics Ground Testing Overview and Lessons Learned

The X-56A Multi-Utility Technology Testbed (MUTT) is a subscale, fixed-wing aircraft designed for high-risk aeroelastic flight demonstration and research. Structural dynamics ground testing for model validation was especially important for this vehicle because the structural model was directly used in the development of a flight control system with active flutter suppression capabilities. Structural dynamics ground tests of the X-56A MUTT with coupled rigid-body and structural modes provided a unique set of challenges. An overview of the ground vibration test (GVT) and moment of inertia (MOI) test setup and execution is presented. The series of GVTs included the wing by itself attached to a strongback and complete vehicle at two mass conditions: empty and full fuel. Two boundary conditions for the complete-vehicle test were studied: on landing gear and suspended free-free. Pitch MOI tests were performed using a compound pendulum method and repeated with two different pendulum lengths for independent verification. The original soft-support test configuration for the GVT used multiple bungees, resulting in unforeseen coupling interactions between the soft-support bungees and the vehicle structural modes. To resolve this problem, the soft-support test setup underwent multiple iterations. The various GVT configurations and boundary-condition modifications are highlighted and explained. Lessons learned are captured for future consideration when performing structural dynamics testing with similar vehicles.

Chin, Alexander W.↗
DOE OSTI2022

Improved Representation of Horizontal Variability and Turbulence in Mesoscale Simulations of an Extended Cold-Air Pool Event

Abstract Cold-air pools (CAPs), or stable atmospheric boundary layers that form within topographic basins, are associated with poor air quality, hazardous weather, and low wind energy output. Accurate prediction of CAP dynamics presents a challenge for mesoscale forecast models in part because CAPs occur in regions of complex terrain, where traditional turbulence parameterizations may not be appropriate. This study examines the effects of the planetary boundary layer (PBL) scheme and horizontal diffusion treatment on CAP prediction in the Weather Research and Forecasting (WRF) Model. Model runs with a one-dimensional (1D) PBL scheme and Smagorinsky-like horizontal diffusion are compared with runs that use a new three-dimensional (3D) PBL scheme to calculate turbulent fluxes. Simulations are completed in a nested configuration with 3-km/750-m horizontal grid spacing over a 10-day case study in the Columbia River basin, and results are compared with observations from the Second Wind Forecast Improvement Project. Using event-averaged error metrics, potential temperature and wind speed errors are shown to decrease both with increased horizontal grid resolution and with improved treatment of horizontal diffusion over steep terrain. The 3D PBL scheme further reduces errors relative to a standard 1D PBL approach. Error reduction is accentuated during CAP erosion, when turbulent mixing plays a more dominant role in the dynamics. Last, the 3D PBL scheme is shown to reduce near-surface overestimates of turbulence kinetic energy during the CAP event. The sensitivity of turbulence predictions to the master length-scale formulation in the 3D PBL parameterization is also explored. Significance Statement In this article, we demonstrate how a new framework for modeling atmospheric turbulence improves cold pool predictions, using a case study from January 2017 in the Columbia River basin (U.S. Pacific Northwest). Cold pools are regions of cold, stagnant air that form within valleys or basins, and improved forecasts could help to mitigate the risks they pose to air quality, transportation, and wind energy production. For the chosen case study, our tests show a reduction in temperature and wind speed errors by up to a factor of 2–3 relative to standard model options. These results strongly motivate continued development of the framework as well as its application to other complex weather events.

17 WIND ENERGY↗
NASA NTRS1992

Numerical analysis of hypersonic turbulent film cooling flows

As a building block, numerical capabilities for predicting heat flux and turbulent flowfields of hypersonic vehicles require extensive model validations. Computational procedures for calculating turbulent flows and heat fluxes for supersonic film cooling with parallel slot injections are described in this study. Two injectant mass flow rates with matched and unmatched pressure conditions using the database of Holden et al. (1990) are considered. To avoid uncertainties associated with the boundary conditions in testing turbulence models, detailed three-dimensional flowfields of the injection nozzle were calculated. Two computational fluid dynamics codes, GASP and FDNS, with the algebraic Baldwin-Lomax and k-epsilon models with compressibility corrections were used. It was found that the B-L model which resolves near-wall viscous sublayer is very sensitive to the inlet boundary conditions at the nozzle exit face. The k-epsilon models with improved wall functions are less sensitive to the inlet boundary conditions. The testings show that compressibility corrections are necessary for the k-epsilon model to realistically predict the heat fluxes of the hypersonic film cooling problems.

Chen, Y. S.↗
NASA NTRS1981

Unsteady environments and responses of the shuttle combined loads orbiter test

Both separate and combined wind tunnel and vibration shaker tests were conducted on two structural panels representative of the Shuttle orbiter in the NASA LaRC 8-foot transonic pressure tunnel (TPT) to determine the effects of combined loads on the thermal protection system (TPS). The primary objective of this test was to provide a combined full-scale load environment and realistic time history of the dynamic pressures, Mach numbers (through transonic), and dynamic structural responses of these panels. The panels were selected from orbiter locations where interactive load sources such as aerodynamic shock waves, turbulent boundary layers, strut-induced vorticity, and substrate deformation combined to provide high bonding loads between the TPS and the orbiter structure. The test panels for this program were highly instrumented with static and dynamic pressure gages, accelerometers, deflectometers, strain gages, Schileren and high-speed photography, and special instrumentation necessary to determine TPS/structure interface loads and tile motions. Two test specimens of each orbiter panel were utilized. Both were high-fidelity representations of the selected orbiter location. The first panel (calibration panel) for each test was utilized to calibrate the tunnel and shaker system as well as provide a data-gathering source for applied static and dynamic loads. The second panel (test panel) for each test was subjected to flight time exposures only so that TPS characteristics as a function of flight mission exposure could be determined. These tests were the last in a series that was suce6ssfully completed and necessary to certify the TPS prior to the first launching of the Shuttle.

P. H. Schuetz↗
DOE OSTI2021

Computational Characterization and Boundary Condition Models of a Jet Interaction Fluidic Oscillator

Fluidically oscillating jet actuators that generate sweeping jets when supplied with a pressurized fluid have been used to mitigate separation and reduce drag in a range of flow control applications. Their implementation in future flight vehicles would require fundamental understanding and accurate predictive techniques of the physics of their internal flow and jet formation. The present investigations focus on high-fidelity, time-accurate simulations to characterize the flow physics of the actuator in quiescent conditions. An important element of the present simulations is to demonstrate the ability of the computational fluid dynamics (CFD) solver to predict the jet characteristics and provide a basis for the development of improved boundary conditions (BC) without entirely resolving the geometrical features of the fluidic device. The CFD-predicted oscillation frequencies of the engendered jets were found to be in excellent agreement with experiments, even on two-dimensional meshes. The study revealed that three-dimensional simulations are required to capture some of the flow features of the sweeping jet such as the double peak in time-averaged velocity distributions downstream of the actuator’s orifice that were measured in experiments. Furthermore, several approaches for modeling the actuator were implemented and assessed in quiescent conditions. The evaluation of a boundary condition at the device throat, based on the phase-averaged flow variables, provides the basis for devising surface-based boundary conditions. The influence and necessity of including turbulent characteristics as part of the boundary conditions have also been identified.

42 ENGINEERING↗
DOE OSTI2021

Fast Na diffusion and anharmonic phonon dynamics in superionic Na 3 PS 4

The design of new solid electrolytes (SEs) hinges on identifying and tuning relevant descriptors. Phonons describe the atomic dynamics in crystalline materials and provide a basis to encode possible minimum energy pathways for ion migration, but anharmonic effects can be large in SEs. Identifying and controlling the pertinent phonon modes coupled most strongly with ionic conductivity, and assessing the role of anharmonicity, could therefore pave the way for discovering and designing new SEs via phonon engineering. In this work, we investigate phonons in Na 3 PS 4 and their coupling to fast Na diffusion, using a combination of neutron scattering, ab initio molecular dynamics (AIMD), and extended molecular dynamics based on machine-learned potentials. We identify that anharmonic soft-modes at the Brillouin zone boundary of the anharmonically stabilized cubic phase constitute key phonon modes that control the Na diffusion process in Na 3 PS 4 . We demonstrate how these strongly anharmonic phonon modes enable Na-ions to hop along the minimum energy pathways. Further, the quasi-elastic neutron scattering (QENS) measurements, supplemented with ab initio and machine-learned molecular dynamics simulations, probe the Na diffusivity and diffusion mechanism. These results offer detailed microscopic insights into the dynamic mechanism of fast Na diffusion and provide an avenue to search for further Na solid electrolytes.

42 ENGINEERING↗
DOE OSTI2022

Anharmonic thermodynamic properties and phase boundary across the postperovskite transition in MgSi O 3

To address the effects of lattice anharmonicity across the perovskite to postperovskite transition in MgSiO 3 , we conduct calculations using the phonon quasiparticle (PHQ) approach. The PHQ is based on abinitio molecular dynamics and, in principle, captures full anharmonicity. Free energies in the thermodynamic limit (N → ∞) are computed using temperature-dependent quasiparticle dispersions within the phonon gas model. Systematic results on anharmonic thermodynamic properties and phase boundary are reported. Both the local density approximation and the generalized gradient approximation calculations are performed to provide confident constraints on these properties. Anharmonic effects are demonstrated by comparing results with those obtained using the quasiharmonic approximation (QHA). The inadequacy of the QHA is indicated by its overestimation of thermal expansivity and thermodynamic Grüneisen parameter and its converged isochoric heat capacity in the high-temperature limit. The PHQ phase boundary has a Clapeyron slope (dP/dT) that increases with temperature. This result contrasts with the nearly zero curvature of the QHA phase boundary. Anharmonicity bends the phase boundary to lower temperatures at high pressures. Implications for the double-crossing of the phase boundary by the mantle geotherm are discussed.

36 MATERIALS SCIENCE↗
NASA NTRS2003

Discrete Abstractions of Hybrid Systems: Verification of Safety and Application to User-Interface Design

Human interaction with a complex control system involves the user, the automation’s discrete mode logic, and the underlying continuous dynamics of the physical system. The user-interface of such systems always displays a reduced set of information about the entire system. Designing interfaces such that all the pertinent information is available and assuring that this information is correct is important for any user-interface, but especially so for safety-critical systems such as automotive systems and autopilots. Here we describe a methodology for the analysis of hybrid control systems that incorporate user interaction, with the goal of assuring that the information provided to the user is correct. That is, the user-interface must contain all information necessary to safely complete a desired procedure or task. We begin with a hybrid system model which incorporates discrete mode logic as well as nonlinear continuous dynamics. Using a hybrid computational tool for reachability, we find the largest region of the state-space in which we can guarantee the state of the system can always remain – this is the safe region of operation. By implementing a controller for safety which arises from this computation, we mathematically guarantee that this safe region is invariant, meaning that the system will always remain within the safe region if the determined controller is used on the boundary of the safe region. Verification within a hybrid framework allows us to account for the continuous dynamics underlying the discrete representations displayed to the user. Using the computed invariant regions as discrete states, we can abstract a discrete event system from this hybrid system with safety restrictions. This abstraction can be used to determine what information must be provided on the display. Furthermore, in cases in which an interface already exists, the abstraction provides the necessary input into existing interface verification methods. We provide two examples: a car traveling through a yellow light at an intersection and an aircraft autopilot in an automatic landing/go-around maneuver. The examples demonstrate the applicability of this methodology to hybrid systems that have operational constraints we can pose in terms of safety. This methodology differs from existing work in hybrid system verification in that we directly account for the user’s interactions with the system.

Meeko Oishi↗
NASA NTRSDate not supplied

Generalized Predictive Control for Active Stability Augmentation and Vibration Reduction on an Aeroelastic Tiltrotor Model

Tiltrotor aircraft are defining the state-of-the-art in vertical lift technology as they have the potential to greatly expand rotary-wing operational boundaries. However, they are often limited in forward flight speed due to complex coupled rotor and wing dynamic instabilities. The U.S. Army and NASA have been developing a new wind tunnel model, the TiltRotor Aeroelastic Stability Testbed(TRAST), to test proprotors in the NASA Langley Research Center Transonic Dynamics Tunnel (TDT) to investigate aeroelastic stability in cruise. The test is intended to provide high-quality research data for analytical tool development and validation. In addition, the TRAST model will support, develop, and mature new technologies for the design of advanced proprotor aircraft. Stability augmentation and vibration reduction during testing is planned with the use of an active control methodology known as Generalized Predictive Control(GPC). GPC is an autoregressive control law that experimentally acquires a system identification to derive the input-output relation of controls and corresponding sensors. This type of control law is especially useful for complex dynamic interactions that are difficult to explicitly model such as proprotor pylon instability, often referred to as whirl flutter. GPC has been successfully employed on other tiltrotor vehicles to suppress whirl flutter instabilities and vibrations. To aid in the characterization of the wind-tunnel model and in tool development, an analytical representation of the wind-tunnel model was developed using the rotorcraft comprehensive analysis system (RCAS) that simulates structural dynamics and aerodynamics. RCAS was used to derive state-space estimates of the physical plant at various flight conditions to test control law effectiveness. This paper will present an overview of the test article development, a description of RCAS, an explanation of the GPC methodology, and results of GPC being applied to state-space plant estimates of the TRAST model. In these simulations, GPC was effective at stabilizing the aircraft beyond the whirl-flutter boundary while simultaneously reducing vibrations across the flight regime. Additionally, a modern advancement to GPC, termed advanced GPC (AGPC), is introduced that enables a self-adapting system identification. Preliminary results show that AGPC is successful at self-correction as the plant changes from what was used for system identification.

tiltrotor↗
NASA NTRS2013

Approach to Integrate Global-Sun Models of Magnetic Flux Emergence and Transport for Space Weather Studies

The Sun lies at the center of space weather and is the source of its variability. The primary input to coronal and solar wind models is the activity of the magnetic field in the solar photosphere. Recent advancements in solar observations and numerical simulations provide a basis for developing physics-based models for the dynamics of the magnetic field from the deep convection zone of the Sun to the corona with the goal of providing robust near real-time boundary conditions at the base of space weather forecast models. The goal is to develop new strategic capabilities that enable characterization and prediction of the magnetic field structure and flow dynamics of the Sun by assimilating data from helioseismology and magnetic field observations into physics-based realistic magnetohydrodynamics (MHD) simulations. The integration of first-principle modeling of solar magnetism and flow dynamics with real-time observational data via advanced data assimilation methods is a new, transformative step in space weather research and prediction. This approach will substantially enhance an existing model of magnetic flux distribution and transport developed by the Air Force Research Lab. The development plan is to use the Space Weather Modeling Framework (SWMF) to develop Coupled Models for Emerging flux Simulations (CMES) that couples three existing models: (1) an MHD formulation with the anelastic approximation to simulate the deep convection zone (FSAM code), (2) an MHD formulation with full compressible Navier-Stokes equations and a detailed description of radiative transfer and thermodynamics to simulate near-surface convection and the photosphere (Stagger code), and (3) an MHD formulation with full, compressible Navier-Stokes equations and an approximate description of radiative transfer and heating to simulate the corona (Module in BATS-R-US). CMES will enable simulations of the emergence of magnetic structures from the deep convection zone to the corona. Finally, a plan will be summarized on the development of a Flux Emergence Prediction Tool (FEPT) in which helioseismology-derived data and vector magnetic maps are assimilated into CMES that couples the dynamics of magnetic flux from the deep interior to the corona.

Mansour, Nagi N.↗