Search NASA⌕ Search

SEARCH · Search NASA

Results for “coupled problems”

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 451 records · Page 25

The Regolith Biters: A Divide-And-Conquer Architecture for Sample-Return Missions

The collective interaction of simple systems can be leveraged to attain complex goals. Based on this principle, we envision space system architectures where the core functional components are decoupled, autonomous, and cooperative. We aim to pursue this vision in the context of small-body sample-return missions. After all, no experimental study sheds more light into our understanding of the origin and evolution of the Solar System than the analysis of samples from asteroids and comets. We also believe that their study is important from a strategic perspective: meteorite impacts pose a direct and credible threat to life on Earth, and the development of contingency small-body deflection missions presupposes some knowledge of the target body. The current architectural paradigm for sample-return missions is centered around a design where spacecraft and sampling device are merged into a single, complex system. We argue that this monolithic approach couples the navigation and sample-collection problems, making both more difficult. In contrast, we propose a decoupled system based on the coordinated interaction between a spacecraft and a collective of small, simple devices - the Regolith Biters (RBs). A spacecraft carrying a number of RBs would travel to the vicinity of a small body. From a favorable vantage point, and while remaining within a safe distance in a non-colliding trajectory, it would release the RBs towards the target body. Upon encountering the body, they would bite the regolith (thus retaining a sample), and eject back to orbit. The spacecraft, being endowed with appropriate navigation and tracking capabilities, would rendezvous with and collect those RBs within its reach, and bring them back to Earth. Separating the navigation and sampling concerns removes the need for proximity operations with the small body-the stage in current architectures that carries the most challenges and risks. Eliminating the need for proximity operations brings back to the discussion the exploration of exciting prospects, like highly active comets, fast-rotating bodies, and binary systems. Distributing the sampling problem among a collective of agents provides the opportunity to sample multiple regions in a single mission. It also provides robustness to various environmental conditions, and may enable the distributed, in situ characterization of the body. In the search for reliability, current architectures rely on complexity: an elaborate system should succeed at once. We rely on numbers: a given agent may fail at any stage, but success is attained by the collective.

Asteroids↗

Integrated Optimization of Mars Hybrid Solar-Electric/Chemical Propulsion Trajectories

NASA’s Human Exploration and Operation Mission Directorate is developing a reusable hybrid transportation architecture in which both chemical and solar-electric propulsion systems are used to deliver crew and cargo to the Martian sphere of influence. By combining chemical and solar-electric propulsions into a single spacecraft and applying each where it is the most effective, the hybrid architecture enables a series of Mars trajectories that are more fuel efficient than an all chemical propulsion architecture without significant increase to trip time. Solving the complex problem of low-thrust trajectory optimization coupled with the vehicle sizing requires development of an integrated trajectory analysis frame- work. Previous studies have utilized a more segmented optimization framework due to the limitation of the tools available. A new integrated optimization framework was recently developed to address the deficiencies of the previous methods that enables higher fidelity analysis to be performed and increases the efficiency of large design space explorations.

Chai, Patrick R.↗

Coupled Aeropropulsive Design Optimization of a Podded Electric Propulsor

New aircraft concepts are increasingly relying on non-traditional propulsion systems to achieve lower energy consumption. These non-traditional technologies, such as boundary layer ingestion or distributed electric propulsion, require a tight integration of the propulsion sys- tem to the airframe, and therefore, a traditional design approach where the aerodynamics and propulsion disciplines are considered separately is likely to result in sub-optimal designs. As a result, we have to rely on coupled aeropropulsive design optimization, in which fully coupled aeropropulsive models are optimized to maximize the advantages of these tightly integrated propulsion systems. Despite its advantages, aeropropulsive design optimization is a relatively new field, and significant advancements are required for wide adoption of this approach, especially in the context of robustness. In this work, we introduce two fully coupled aeropropulsive design optimization approaches that are compatible with CFD models with body-force terms and boundary conditions to model the effects of the propulsion system in the flow field. To test these new approaches, we developed a podded electric fan model based on the aft-propulsor of NASA’s STARC-ABL concept. This simple design problem enables us to rapidly develop and test different aeropropulsive coupling approaches, while including the challenging physical interactions that arise from coupling CFD and propulsion models. The coupled aeropropulsive model is implemented using the MPhys library, which is built with NASA’s OpenMDAO frame- work. Using this benchmark design problem, we demonstrate the robustness of the approaches and our aeropropulsive design optimization framework by performing a sweep of optimizations for a total of 50 CFD-based aeropropulsive design optimizations. Furthermore, the new aero- propulsive coupling approaches enable multi-point design optimizations. We demonstrate this capability in a multi-point design optimization problem where we optimize cruise performance subject to a fan-face distortion constraint at rolling take-off conditions. The aeropropulsive modeling approaches we present in this work represent a significant milestone in the field of aeropropulsive design optimization. The framework we developed is extremely robust and flexible, and these developments will be crucial for future aeropropulsive design optimizations of complete turbofan engines.

Optimization propulsion aerodynamics↗

Convergence Criteria for Multiphysics Simulations

The behavior of engineered systems is often influenced by multiple physical phenomena, such as mechanical deformation, heat transfer, and chemical species transport and reactions. There are often strong interactions between these phenomena, and there is increasing interest in applying coupled-physics models to improve understanding of physical behavior under complex environmental conditions. Multiple simulation frameworks that facilitate coupled-physics simulations are in widespread use, and these employ a variety of techniques to account for interactions between those physics. Many frameworks solve the physics models independently and transfer results between them. Alternatively, a single monolithic system of equations for every physics model can be formed and solved. Each of these approaches has its benefits and drawbacks, and the optimal approach varies depending on the nature of the problem. The open-source MOOSE framework was developed targeting solution of large-scale multiphysics problems. Although it provides options for all these coupling approaches, its standard approach for multiphysics solutions is to form and solve a single monolithic system of equations containing the unknowns for all physics models. MOOSE provides a streamlined approach for users to define the solution variables, the terms in the partial differential equations pertaining to each variable, and interactions between solution variables. One aspect of the monolithic solution approach that can be problematic, however, is defining appropriate convergence criteria for the nonlinear system. A standard approach is to determine convergence is to simply take a norm of the residual vector corresponding to the full vector of unknowns. However, if the residual vector contains variables for multiple physics models, the magnitudes of those variables can differ significantly, and the variables can converge at significantly different rates from each other. It is important to ensure that the variables for each of the physics are converged, and also ensure that the convergence criteria are not excessively stringent in cases when there is little change in the solution. This talk presents representative multiphysics problems to highlight these issues, and shows strategies for convergence criteria in MOOSE that are robust for multiphysics models under a variety of conditions.

97 - MATHEMATICS AND COMPUTING↗

Dynamics in Layer Models of Solid Flame Propagation

Self-propagating high-temperature synthesis (SHS) is a process in which combustion waves, e.g., "solid flames", which are considered here, are employed to synthesize desired materials. Like many other systems, SHS is a pattern forming system. The problem of describing experimentally observed patterns and of predicting new, as yet unobserved, patterns continues to attract the attention of scientists and mathematicians due to the fundamental significance of the phenomena in combustion in particular, and in nonlinear science in general. Here, we analyze the dynamics of solid flame propagation in a 2D region by considering the region to be composed of parallel, identical layers aligned along the direction of propagation and having thermal contact. Each layer is then described by wave propagation in 1D, with the transverse Laplacian replaced by a term describing heat exchange between neighboring layers. This configuration is the simplest model of a 2D system because it accounts, in a simple way, for the principal feature of the problem, i.e., heat exchange between neighbors in the transverse direction. For simplicity, we describe the situation for two layers. Because the layers are identical, uniformly propagating waves in each layer must be identical, independent of the heat exchange rate alpha. When the Zeldovich number Z exceeds a critical value Z(sub c), which depends on alpha, uniformly propagating waves become unstable. The stability diagram for the two coupled layers reproduces that for the full 2D problem after appropriate identification of parameters in the two problems. Depending on parameter values, we determine three different steady-state dynamical behaviors (though additional behaviors are also expected to occur). The three behaviors are: (i) waves in each layer which pulsate in phase as they propagate, so that together they form a single pulsating propagating wave; (ii) the waves in each layer are no longer identical, and antiphase pulsations occur, with the two waves alternately advancing and receding as they propagate. This wave is the analog of the spinning wave on the surface of a circular cylinder; (iii) finally, there is a region of bistability between the in phase and antiphase waves. with each having its own domain of attraction, so that which of the two behaviors occur depends on the condition of initiation of the wave. The results of our computations indicate a qualitative similarity in the behavior of combustion waves in the layer model and in the full 2D model. Specifically, due to the similarity between the small alpha wave behavior in the layer model and the large diameter behavior in the model of waves on the surface of a cylinder, we are able to predict the behavior of the mean velocity for the waves on the cylinder, where computations of the full problem can be rather difficult. We also compute solutions for three or more layers. The results of our computations prompt us to predict that, while planar uniformly propagating waves are unstable, the wave will be quasiplanar, i.e., the resulting spinning waves have very low amplitude hot spots, and travel with the velocity close to that of the uniformly propagating wave. Such waves may be difficult to distinguish from uniformly propagating waves in experiments. We also find that for both the layer model and full 2D problem, steady-state time-dependent waves, e.g., pulsating and spinning wave, have a conserved quantity H which characterizes the excess energy in the wave, just as in the case of uniformly propagating waves. The quantity H, which is generated by dissipation, does not vary in time and is proportional to the diffusivity and caloricity of the system, and inversely proportional to the mean wave velocity.

Aldushin, A. P.↗

Ares-I-X Stability and Control Flight Test: Analysis and Plans

The flight test of the Ares I-X vehicle provides a unique opportunity to reduce risk of the design of the Ares I vehicle and test out design, math modeling, and analysis methods. One of the key features of the Ares I design is the significant static aerodynamic instability coupled with the relatively flexible vehicle - potentially resulting in a challenging controls problem to provide adequate flight path performance while also providing adequate structural mode damping and preventing adverse control coupling to the flexible structural modes. Another challenge is to obtain enough data from the single flight to be able to conduct analysis showing the effectiveness of the controls solutions and have data to inform design decisions for Ares I. This paper will outline the modeling approaches and control system design to conduct this flight test, and also the system identification techniques developed to extract key information such as control system performance (gain/phase margins, for example), structural dynamics responses, and aerodynamic model estimations.

Brandon, Jay M.↗

A Finite Element Procedure for Calculating Fluid-Structure Interaction Using MSC/NASTRAN

This report is intended to serve two purposes. The first is to present a survey of the theoretical background of the dynamic interaction between a non-viscid, compressible fluid and an elastic structure is presented. Section one presents a short survey of the application of the finite element method (FEM) to the area of fluid-structure-interaction (FSI). Section two describes the mathematical foundation of the structure and fluid with special emphasis on the fluid. The main steps in establishing the finite element (FE) equations for the fluid structure coupling are discussed in section three. The second purpose is to demonstrate the application of MSC/NASTRAN to the solution of FSI problems. Some specific topics, such as fluid structure analogy, acoustic absorption, and acoustic contribution analysis are described in section four. Section five deals with the organization of the acoustic procedure flowchart. Section six includes the most important information that a user needs for applying the acoustic procedure to practical FSI problems. Beginning with some rules concerning the FE modeling of the coupled system, the NASTRAN USER DECKs for the different steps are described. The goal of section seven is to demonstrate the use of the acoustic procedure with some examples. This demonstration includes an analytic verification of selected FE results. The analytical description considers only some aspects of FSI and is not intended to be mathematically complete. Finally, section 8 presents an application of the acoustic procedure to vehicle interior acoustic analysis with selected results.

Chargin, Mladen↗

Models of Small-Scale Patchiness

Patchiness is perhaps the most salient characteristic of plankton populations in the ocean. The scale of this heterogeneity spans many orders of magnitude in its spatial extent, ranging from planetary down to microscale. It has been argued that patchiness plays a fundamental role in the functioning of marine ecosystems, insofar as the mean conditions may not reflect the environment to which organisms are adapted. Understanding the nature of this patchiness is thus one of the major challenges of oceanographic ecology. The patchiness problem is fundamentally one of physical-biological-chemical interactions. This interconnection arises from three basic sources: (1) ocean currents continually redistribute dissolved and suspended constituents by advection; (2) space-time fluctuations in the flows themselves impact biological and chemical processes, and (3) organisms are capable of directed motion through the water. This tripartite linkage poses a difficult challenge to understanding oceanic ecosystems: differentiation between the three sources of variability requires accurate assessment of property distributions in space and time, in addition to detailed knowledge of organismal repertoires and the processes by which ambient conditions control the rates of biological and chemical reactions. Various methods of observing the ocean tend to lie parallel to the axes of the space/time domain in which these physical-biological-chemical interactions take place. Given that a purely observational approach to the patchiness problem is not tractable with finite resources, the coupling of models with observations offers an alternative which provides a context for synthesis of sparse data with articulations of fundamental principles assumed to govern functionality of the system. In a sense, models can be used to fill the gaps in the space/time domain, yielding a framework for exploring the controls on spatially and temporally intermittent processes. The following discussion highlights only a few of the multitude of models which have yielded insight into the dynamics of plankton patchiness. In addition, this particular collection of examples is intended to furnish some exposure to the diversity of modeling approaches which can be brought to bear on the problem. These approaches range from abstract theoretical models intended to elucidate specific processes, to complex numerical formulations which can be used to actually simulate observed distributions in detail.

McGillicuddy, D. J.↗

Receptivity in parallel flows: An adjoint approach

Linear receptivity studies in parallel flows are aimed at understanding how external forcing couples to the natural unstable motions which a flow can support. The vibrating ribbon problem models the original Schubauer and Skramstad boundary layer experiment and represents the classic boundary layer receptivity problem. The process by which disturbances are initiated in convectively-unstable jets and shear layers has also received attention. Gaster was the first to handle the boundary layer analysis with the recognition that spatial modes, rather than temporal modes, were relevant when studying convectively-unstable flows that are driven by a time-harmonic source. The amplitude of the least stable spatial mode, far downstream of the source, is related to the source strength by a coupling coefficient. The determination of this coefficient is at the heart of this type of linear receptivity study. The first objective of the present study was to determine whether the various wave number derivative factors, appearing in the coupling coefficients for linear receptivity problems, could be reexpressed in a simpler form involving adjoint eigensolutions. Secondly, it was hoped that the general nature of this simplification could be shown; indeed, a rather elegant characterization of the receptivity properties of spatial instabilities does emerge. The analysis is quite distinct from the usual Fourier-inversion procedures, although a detailed knowledge of the spectrum of the Orr-Sommerfeld equation is still required. Since the cylinder wake analysis proved very useful in addressing control considerations, the final objective was to provide a foundation upon which boundary layer control theory may be developed.

Hill, D. Christopher↗

A boundary integral method for an inverse problem in thermal imaging

An inverse problem in thermal imaging involving the recovery of a void in a material from its surface temperature response to external heating is examined. Uniqueness and continuous dependence results for the inverse problem are demonstrated, and a numerical method for its solution is developed. This method is based on an optimization approach, coupled with a boundary integral equation formulation of the forward heat conduction problem. Some convergence results for the method are proved, and several examples are presented using computationally generated data.

Bryan, Kurt↗

Control considerations for CCV fighters at high angles of attack

Wind tunnel and piloted simulation studies were conducted to investigate the potential high angle of attack control problems that are introduced by the use of the CCV concept of relaxed static pitch stability (RSS) on fighter aircraft. A conventional wing/aft tail design incorporating modest levels of static instability and a close-coupled canard/wing design exhibiting very high levels of instability was investigated. Two types of high angle of attack control problems can result from the use of RSS: pitch departures caused by coupling and deep stall trim. Avoidance of these problems requires that the airplane have sufficient nose-down pitch control at high angles of attack. The effectiveness of several pitch control configurations were investigated including conventional aft-mounted stabilators, wing-mounted elevators, canard-mounted flaps, and all-moveable canards. Varying the incidence of the canards was the most effective scheme; however, very large deflections may be required on highly unstable configurations to prevent pitch departure without sacrificing roll performance and to avoid deep stall trim. For situations where the high angle of attack pitch control requirement is not met, control laws were developed to inhibit the departure and to allow deep stall recovery. However, these schemes involve limiting airplane roll capability and therefore can potentially compromise maneuverability.

Nguyen, L. T.↗

Moment method with isoparametric elements for three-dimensional anisotropic scatterers

A novel method for computing the frequency-domain electromagnetic fields scattered from, and penetrating into, arbitrarily shaped, three-dimensional, lossy, inhomogeneous anisotropic scatters is presented. The method is based on a general volume integrodifferential formulation of the scattering problem and consists of the numerical solution of the coupled integral equations by the moment method and point matching. The numerical model of the scatterer is obtained by parametric volume elements, and the basis functions used to represent the field within each element are the same used in the finite-element method. Element integration problems due to the singular kernel of the integral equations are treated in some detail. Numerical results for both the isotropic and the anisotropic spherical scatterer are presented, including comparisons with results obtained by different numerical methods for the isotropic cases considered. The capability of the numerical code presented to deal with cases where the material parameters of the scatterer are given by singular matrices is discussed for two particular examples.

Graglia, Roberto D.↗

Theoretical Study of CH3Cl-N2 Line Shapes in the ν1 Band. Line Mixing Effects in QR Doublets and QQk Sub-Branches

Line coupling and line mixing effects of CH3Cl lines in the ν1 band perturbed by N2 have been investigated. We have taken into account the non-diagonality of the exp(−S2) operator within specified line spaces as well as the k-degeneracy of the transitions (due to the double degeneracy of the j,k levels with k ≠ 0). These transitions should be considered as doublets coupled by the line mixing process. A new problem appears in the calculation when the atom-atom potential model is introduced. In order to overcome this difficulty, a pragmatic approach is proposed. Comparisons of theoretically calculated matrix elements of W with measurements of QR(j,k) doublets as well as some QQk sub-branches, which are strongly affected by line mixing, have been made. The results show that the formalism improved in this way leads to rather accurate predictions.

Methyl chloride↗

Coupled Lindblad Pseudomode Theory for Simulating Open Quantum Systems

Coupled Lindblad pseudomode theory is a promising approach for simulating non-Markovian quantum dynamics on both classical and quantum platforms, with dynamics that can be realized as a quantum channel. We provide theoretical evidence that the number of coupled pseudomodes only needs to scale as polylog⁡(𝑇/𝜖) in the simulation time 𝑇 and precision 𝜖. Inspired by the realization problem in control theory, we also develop a robust numerical algorithm for constructing the coupled modes that avoid the nonconvex optimization required by existing approaches. We demonstrate the effectiveness of our method by computing population dynamics and absorption spectra for the spin-boson model. Furthermore, this Letter provides a significant theoretical and computational improvement to the coupled Lindblad framework, which impacts a broad range of applications from classical simulations of quantum impurity problems to quantum simulations on near-term quantum platforms.

Anderson impurity model↗

Radiation Hardened Engineered Substrates for Time and Space Resolution

Low Gain Avalanche Diodes (LGADs) have become sensors of choice for fast timing of minimum ionizing particles. The current generation of these devices suffer from only moderate radiation hardness, low fill factor, and limited design options. Cactus Materials Inc. proposes to develop sensors with gain layer implants buried beneath the surface which will, with AC coupling developed by Brookhaven National Laboratory (BNL) and University of California Santa Cruz (UC Santa Cruz) solve these problems. Fermilab has a major role in the CMS timing upgrade using the current generation of LGADs and development of AC-coupled, buried layer LGADs will solve some of the central problems inherent in this technology. Annex B included a two-year work plan, deliverables, schedule and funding for the project’s SBIR Phase II award to develop and test prototype LGADs utilizing wafer bonding technology and epitaxial growth, that will allow for fabrication of radiation hard AC coupled devices with 100% fill factor and adjustable operating point.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Generalized multi-point inverse airfoil design

In a rather general sense, inverse airfoil design can be taken to mean the problem of specifying a desired set of airfoil characteristics, such as the airfoil maximum thickness ratio, pitching moment, part of the velocity distribution or boundary-layer development, etc., then from this information determine the corresponding airfoil shape. This paper presents a method which approaches the design problem from this perspective. In particular, the airfoil is divided into segments along which, together with the design conditions, either the velocity distribution or boundary-layer development may be prescribed. In addition to these local desired distributions, single parameters like the airfoil thickness can be specified. The problem of finding the airfoil shape is determined by coupling an incompressible, inviscid, inverse airfoil design method with a direct integral boundary-layer analysis method and solving the resulting nonlinear equations via a multidimensional Newton iteration technique. The approach is fast and easily allows for interactive design. It is also flexible and could be adapted to solving compressible, inverse airfoil design problems.

Selig, Michael S.↗

The development and preliminary application of an invariant coupled diffusion and chemistry model

In many real-world pollution chemical reaction problems, the rate of reaction problems, the rate of reaction may be greatly affected by unmixedness. An approximate closure scheme for a chemical kinetic submodel which conforms to the principles of invariant modeling and which accounts for the effects of inhomogeneous mixing over a wide range of conditions has been developed. This submodel has been coupled successfully with invariant turbulence and diffusion models, permitting calculation of two-dimensional diffusion of two reacting (isothermally) chemical species. The initial calculations indicate the ozone reactions in the wake of stratospheric aircraft will be substantially affected by the rate of diffusion of ozone into the wake, and in the early wake, by unmixedness.

Hilst, G. R.↗