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 649 records · Page 36

A surface crack in shells under mixed-mode loading conditions

The present consideration of a shallow shell's surface crack under general loading conditions notes that while the mode I state can be separated, modes II and III remain coupled. A line spring model is developed to formulate the part-through crack problem under mixed-mode conditions, and then to consider a shallow shell of arbitrary curvature having a part-through crack located on the outer or the inner surface of the shell; Reissner's transverse shear theory is used to formulate the problem under the assumption that the shell is subjected to all five moment and stress resultants.

Joseph, P. F.↗

Heat pipe dynamic behavior

The vapor flow in a heat pipe was mathematically modeled and the equations governing the transient behavior of the core were solved numerically. The modeled vapor flow is transient, axisymmetric (or two-dimensional) compressible viscous flow in a closed chamber. The two methods of solution are described. The more promising method failed (a mixed Galerkin finite difference method) whereas a more common finite difference method was successful. Preliminary results are presented showing that multi-dimensional flows need to be treated. A model of the liquid phase of a high temperature heat pipe was developed. The model is intended to be coupled to a vapor phase model for the complete solution of the heat pipe problem. The mathematical equations are formulated consistent with physical processes while allowing a computationally efficient solution. The model simulates time dependent characteristics of concern to the liquid phase including input phase change, output heat fluxes, liquid temperatures, container temperatures, liquid velocities, and liquid pressure. Preliminary results were obtained for two heat pipe startup cases. The heat pipe studied used lithium as the working fluid and an annular wick configuration. Recommendations for implementation based on the results obtained are presented. Experimental studies were initiated using a rectangular heat pipe. Both twin beam laser holography and laser Doppler anemometry were investigated. Preliminary experiments were completed and results are reported.

Issacci, F.↗

Unified Optimization Of A Structure And Its Control System

Method for unified optimization of structure and its vibration-suppressing control system involves extension of established techniques of optimization via approach based on homotopy, in which cost function used as criterion for optimization of control system coupled to corresponding cost function of structure via new parameter. Applicable to such terrestrial problems as designing buildings with control actuators to prevent damage by earthquakes.

Milman, Mark H.↗

Control/structure interaction during Space Station Freedom-Orbiter berthing

The berthing maneuver is essential for the construction and assembly of Space Station Freedom (SSF) and has a direct effect on the SSF assembly build up and SSF/Orbiter operations. The effects of flexible body dynamics coupled with the available control system may impose new requirements on the maneuver. The problem is further complicated by the effect of the SSF control system on the Shuttle Remote Manipulator System (SRMS). These effects will play a major role in the development of operational requirements which need to be identified and validated in order to assure total safety and maneuver execution during SSF construction. This paper presents the results of ongoing studies to investigate the Control/Structure Interaction (CSI) during the berthing operations. The problem is formulated in terms of multi-flex body equations of motion for SSF and the SRMS and on-orbit flight control systems for the SRMS and the SSF, which includes the Control Moment Gyro (CMG) and Reaction Control System (RCS) Attitude Control Systems (ACS). The SSF control system designs are based on the Preliminary Design Review (PDR) version of the Honeywell design. The simulation tool used for the analysis is briefly described and the CSI results are presented for given berthing scenarios.

Hua, Tuyen↗

An analysis of supersonic flows with low-Reynolds number compressible two-equation turbulence models using LU finite volume implicit numerical techniques

A generalized flow solver using an implicit Lower-upper (LU) diagonal decomposition based numerical technique has been coupled with three low-Reynolds number kappa-epsilon models for analysis of problems with engineering applications. The feasibility of using the LU technique to obtain efficient solutions to supersonic problems using the kappa-epsilon model has been demonstrated. The flow solver is then used to explore limitations and convergence characteristics of several popular two equation turbulence models. Several changes to the LU solver have been made to improve the efficiency of turbulent flow predictions. In general, the low-Reynolds number kappa-epsilon models are easier to implement than the models with wall-functions, but require much finer near-wall grid to accurately resolve the physics. The three kappa-epsilon models use different approaches to characterize the near wall regions of the flow. Therefore, the limitations imposed by the near wall characteristics have been carefully resolved. The convergence characteristics of a particular model using a given numerical technique are also an important, but most often overlooked, aspect of turbulence model predictions. It is found that some convergence characteristics could be sacrificed for more accurate near-wall prediction. However, even this gain in accuracy is not sufficient to model the effects of an external pressure gradient imposed by a shock-wave/ boundary-layer interaction. Additional work on turbulence models, especially for compressibility, is required since the solutions obtained with base line turbulence are in only reasonable agreement with the experimental data for the viscous interaction problems.

Lee, J.↗

Single-drop reactive extraction/extractive reaction with forced convective diffusion and interphase mass transfer

An algorithm has been developed for the forced convective diffusion-reaction problem for convection inside and outside a droplet by a recirculating flow field hydrodynamically coupled at the droplet interface with an external flow field that at infinity becomes a uniform streaming flow. The concentration field inside the droplet is likewise coupled with that outside by boundary conditions at the interface. A chemical reaction can take place either inside or outside the droplet or reactions can take place in both phases. The algorithm has been implemented and results are shown here for the case of no reaction and for the case of an external first order reaction, both for unsteady behavior. For pure interphase mass transfer, concentration isocontours, local and average Sherwood numbers, and average droplet concentrations have been obtained as a function of the physical properties and external flow field. For mass transfer enhanced by an external reaction, in addition to the above forms of results, we present the enhancement factor, with the results now also depending upon the (dimensionless) rate of reaction.

Kleinman, Leonid S.↗

Distributed Adaptive Control: Beyond Single-Instant, Discrete Variables

In extensive form noncooperative game theory, at each instant t, each agent i sets its state x, independently of the other agents, by sampling an associated distribution, q(sub i)(x(sub i)). The coupling between the agents arises in the joint evolution of those distributions. Distributed control problems can be cast the same way. In those problems the system designer sets aspects of the joint evolution of the distributions to try to optimize the goal for the overall system. Now information theory tells us what the separate q(sub i) of the agents are most likely to be if the system were to have a particular expected value of the objective function G(x(sub 1),x(sub 2), ...). So one can view the job of the system designer as speeding an iterative process. Each step of that process starts with a specified value of E(G), and the convergence of the q(sub i) to the most likely set of distributions consistent with that value. After this the target value for E(sub q)(G) is lowered, and then the process repeats. Previous work has elaborated many schemes for implementing this process when the underlying variables x(sub i) all have a finite number of possible values and G does not extend to multiple instants in time. That work also is based on a fixed mapping from agents to control devices, so that the the statistical independence of the agents' moves means independence of the device states. This paper also extends that work to relax all of these restrictions. This extends the applicability of that work to include continuous spaces and Reinforcement Learning. This paper also elaborates how some of that earlier work can be viewed as a first-principles justification of evolution-based search algorithms.

Wolpert, David H.↗

Broadband Phase Retrieval for Image-Based Wavefront Sensing

A focus-diverse phase-retrieval algorithm has been shown to perform adequately for the purpose of image-based wavefront sensing when (1) broadband light (typically spanning the visible spectrum) is used in forming the images by use of an optical system under test and (2) the assumption of monochromaticity is applied to the broadband image data. Heretofore, it had been assumed that in order to obtain adequate performance, it is necessary to use narrowband or monochromatic light. Some background information, including definitions of terms and a brief description of pertinent aspects of image-based phase retrieval, is prerequisite to a meaningful summary of the present development. Phase retrieval is a general term used in optics to denote estimation of optical imperfections or aberrations of an optical system under test. The term image-based wavefront sensing refers to a general class of algorithms that recover optical phase information, and phase-retrieval algorithms constitute a subset of this class. In phase retrieval, one utilizes the measured response of the optical system under test to produce a phase estimate. The optical response of the system is defined as the image of a point-source object, which could be a star or a laboratory point source. The phase-retrieval problem is characterized as image-based in the sense that a charge-coupled-device camera, preferably of scientific imaging quality, is used to collect image data where the optical system would normally form an image. In a variant of phase retrieval, denoted phase-diverse phase retrieval [which can include focus-diverse phase retrieval (in which various defocus planes are used)], an additional known aberration (or an equivalent diversity function) is superimposed as an aid in estimating unknown aberrations by use of an image-based wavefront-sensing algorithm. Image-based phase-retrieval differs from such other wavefront-sensing methods, such as interferometry, shearing interferometry, curvature wavefront sensing, and Shack-Hartmann sensing, all of which entail disadvantages in comparison with image-based methods. The main disadvantages of these non-image based methods are complexity of test equipment and the need for a wavefront reference.

Dean, Bruce H.↗

Methane on Mars: New Insights into the Sensitivity of CH4 with the NOMAD/ExoMars Spectrometer Through Its First In-Flight Calibration

The Nadir and Occultation for MArs Discovery instrument (NOMAD), onboard the ExoMars Trace Gas Orbiter (TGO) spacecraft was conceived to observe Mars in solar occultation, nadir, and limb geometries, and will be able to produce an outstanding amount of diverse data, mostly focused on properties of the atmosphere. The infrared channels of the instrument operate by combining an echelle grating spectrometer with an Acousto- Optical Tunable Filter (AOTF). Using in-flight data, we characterized the instrument performance and parameterized its calibration. In particular: an accurate frequency calibration was achieved, together with its variability due to thermal effects on the grating. The AOTF properties and transfer function were also quantified, and we developed and tested a realistic method to compute the spectral continuum transmitted through the coupled grating and AOTF system. The calibration results enabled unprecedented insights into the important problem of the sensitivity of NOMAD to methane abundances in the atmosphere. We also deeply characterized its performance under realistic conditions of varying aerosol abundances, diverse albedos and changing illumination conditions as foreseen over the nominal mission. The results show that, in low aerosol conditions, NOMAD single spectrum, 1σ sensitivity to CH4 is around 0.33 ppbv at 20 km of altitude when performing solar occultation's, and better than 1 ppbv below 30 km. In dusty conditions, we show that the sensitivity drops to 0 below 10 km. In Nadir geometry, results demonstrate that NOMAD will be able to produce seasonal maps of CH4 with a sensitivity around 5 ppbv over most of planet's surface with spatial integration over 5×5° bins. Results show also that such numbers can be improved by a factor of ~10 to ~30 by data binning. Overall, our results quantify NOMAD's capability to address the variable aspects of Martian climate.

Liuzzi, Giuliano↗

Updraft Dynamics and Microphysics: On the Added Value of The Cumulus Thermal Reference Frame in Simulations of Aerosol–Deep Convection Interactions

One fundamental question about atmospheric moist convection processes that remains debated is whether, or under which conditions, a relevant variability in background aerosol concentrations may have a significant dynamical impact on convective clouds and their associated precipitation. Furthermore, current climate models must parameterize both the microphysical and the cumulus convection processes, but this is usually implemented separately, whereas in nature there is a strong coupling between them. As a first step to improve our understanding of these two problems, we investigate how aerosol concentrations modify key properties of updrafts in eight large-eddy-permitting regional simulations of a case study of scattered convection over Houston, Texas, in which convection is explicitly simulated and microphysical processes are parameterized. Dynamical and liquid-phase microphysical responses are investigated using the following two different reference frames: static cloudy updraft grid cells versus tracked cumulus thermals. In both frameworks, we observe the expected microphysical responses to higher aerosol concentrations, such as higher cloud number concentrations and lower rain number concentrations. In terms of the dynamical responses, both frameworks indicate weak impacts of varying aerosol concentrations relative to the noise between simulations over the observationally derived range of aerosol variability for this case study. On the other hand, results suggest that thermals are more selective than cloudy updraft grid cells in terms of sampling the most active convective air masses. For instance, vertical velocity from thermals is significantly higher at upper levels than when sampled from cloudy updraft grid points, and several microphysical variables have higher average values in the cumulus thermal framework than in the cloudy updraft framework. In addition, the thermal analysis is seen to add rich quantitative information about the rates and covariability of microphysical processes spatially and throughout tracked thermal lifecycles, which can serve as a stronger foundation for improving subgrid-scale parameterizations.

aerosols↗

Updraft Dynamics and Microphysics: On the Added Value of the Cumulus Thermal Reference Frame in Simulations of Aerosol–Deep Convection Interactions

One fundamental question about atmospheric moist convection processes that remains debated is whether or under what conditions a relevant variability in background aerosol concentrations may have a significant dynamical impact on convective clouds and their associated precipitation. Furthermore, current climate models must parameterize both the microphysical and the cumulus convection processes, but this is usually implemented separately, whereas in nature there is a strong coupling between them. As a first step to improve our understanding of these two problems, we investigate how aerosol concentrations modify key properties of updrafts in eight large-eddy permitting regional simulations of a case study of scattered convection over Houston, Texas, in which convection is explicitly simulated and microphysical processes are parameterized. Dynamical and liquid-phase microphysical responses are investigated using two different reference frames: static cloudy-updraft grid cells versus tracked cumulus thermals. In both frameworks we observe the expected microphysical responses to higher aerosol concentrations, such as higher cloud number concentrations and lower rain number concentrations. In terms of the dynamical responses, both frameworks indicate weak impacts of varying aerosol concentrations relative to the noise between simulations over the observationally derived range of aerosol variability for this case study. On the other hand, results suggest that thermals are more selective than cloudy-updraft grid cells in terms of sampling the most active convective air masses. For instance, vertical velocity from thermals is significantly higher at upper levels than when sampled from cloudy-updraft grid points, and several microphysical variables have higher average values in the cumulus thermal framework than in the cloudy-updraft framework. In addition, the thermal analysis is seen to add rich quantitative information about the rates and covariability of microphysical processes spatially and throughout tracked thermal lifecycles, which can serve as a stronger foundation for improving subgrid-scale parameterizations.

Daniel Hernandez-Deckers↗

Hardware-Efficient Quantum Optimization Layered Algorithms and Experiments

Quantum optimization algorithms, such as QAOA, that implement parametrized stochastic optimization solvers attempt to identify low-energy solutions of Ising systems by exploiting available quantum effects in noisy-intermediate scale machines. Engineering a well-performing parametrized quantum optimization circuit is indeed an exercise in balancing the trade-off between expressivity and implementation complexity. We show that, for MaxCut QAOA circuits defined on native hardware topology (Rigetti’s Aspen Quantum Processors), error-mitigation techniques recover simulated features of the noiseless theory. Moreover, we explore a design space for QAOA-like ansatze that perform well in theory as well as in hardware for fully-connected problems. We also discuss how efficient coherence and entanglement detection methods that could be coupled with quantum optimization experiments require only linear overhead in benchmarking time.

quantum computing↗

A Variational Assimilation Method for Satellite and Conventional Data: a Revised Basic Model 2B

A variational objective analysis technique that modifies observations of temperature, height, and wind on the cyclone scale to satisfy the five 'primitive' model forecast equations is presented. This analysis method overcomes all of the problems that hindered previous versions, such as over-determination, time consistency, solution method, and constraint decoupling. A preliminary evaluation of the method shows that it converges rapidly, the divergent part of the wind is strongly coupled in the solution, fields of height and temperature are well-preserved, and derivative quantities such as vorticity and divergence are improved. Problem areas are systematic increases in the horizontal velocity components, and large magnitudes of the local tendencies of the horizontal velocity components. The preliminary evaluation makes note of these problems but detailed evaluations required to determine the origin of these problems await future research.

Achtemeier, Gary L.↗

Actuator placement for active sound and vibration control of cylinders

Active structural acoustic control is a method in which the control inputs (used to reduce interior noise) are applied directly to a vibrating structural acoustic system. The control concept modeled in this work is the application of in-plane force inputs to piezoceramic patches bonded to the wall of a vibrating cylinder. The cylinder is excited by an exterior noise source -- an acoustic monopole -- located near the outside of the cylinder wall. The goal is to determine the force inputs and sites for the piezoelectric actuators so that (1) the interior noise is effectively damped; (2) the level of vibration of the cylinder shell is not increased; and (3) the power requirements needed to drive the actuators are not excessive. We studied external monopole excitations at two frequencies. A cylinder resonance of 100 Hz, where the interior acoustic field is driven in multiple, off-resonance cylinder cavity modes, and a cylinder resonance of 200 Hz are characterized by both near and off-resonance cylinder vibration modes which couple effectively with a single, dominant, low-order acoustic cavity mode at resonance. Previous work has focused almost exclusively on meeting objective (1) and solving a complex least-squares problem to arrive at an optimal force vector for a given set of actuator sites. In addition, it has been noted that when the cavity mode couples with cylinder vibration modes (our 200 Hz case) control spillover may occur in higher order cylinder shell vibrational modes. How to determine the best set of actuator sites to meet objectives (1)-(3) is the main contribution of our research effort. The selection of the best set of actuator sites from a set of potential sites is done via two metaheuristics -- simulated annealing and tabu search. Each of these metaheuristics partitions the set of potential actuator sites into two disjoint sets: those that are selected to control the noise (on) and those that are not (off). Next, each metaheuristic attempts to improve this initial solution by calculating the change in the objective value when one selected actuator site is turned off and one actuator site that previously was not selected is turned on. All such pairwise exchanges are performed and the exchange that improves the objective the most is made. Eventually the search is unable to improve the objective value and a local optimum (with respect to pairwise exchanges) is reached. Both simulated annealing and tabu search provide mechanisms to escape local optima and allow the search to continue until (hopefully) a global optimum is found. Our experiments with the 100 Hz and 200 Hz cases confirm that both metaheuristics are able to uncover better solutions than those selected based upon engineering judgement alone. In addition, the high quality solutions generated by these metaheuristics, when minimizing interior noise, do not further excite the cylinder shell. Thus, we are able to meet objective (2) without imposing an additional constraint or forming a multiobjective performance measure. An additional observation is that in many cases the amplitude and phase values for several chosen actuator sites were nearly identical. This natural grouping means that fewer control channels are needed and the resulting control system is simpler. Currently no power requirements have been set, so objective (3) cannot be addressed. A set of experiments is planned with a laboratory test article (a cylinder). For these experiments the transfer matrices will be generated experimentally. It is hoped that the predicted performance of the best actuator sites found by our metaheuristics will correlate well with the measured performance.

Kincaid, Rex K.↗

Computational Issues Associated with Temporally Deforming Geometries Such as Thrust Vectoring Nozzles

During the past decade, computational simulation of fluid flow around complex configurations has progressed significantly and many notable successes have been reported, however, unsteady time-dependent solutions are not easily obtainable. The present effort involves unsteady time dependent simulation of temporally deforming geometries. Grid generation for a complex configuration can be a time consuming process and temporally varying geometries necessitate the regeneration of such grids for every time step. Traditional grid generation techniques have been tried and demonstrated to be inadequate to such simulations. Non-Uniform Rational B-splines (NURBS) based techniques provide a compact and accurate representation of the geometry. This definition can be coupled with a distribution mesh for a user defined spacing. The present method greatly reduces cpu requirements for time dependent remeshing, facilitating the simulation of more complex unsteady problems. A thrust vectoring nozzle has been chosen to demonstrate the capability as it is of current interest in the aerospace industry for better maneuverability of fighter aircraft in close combat and in post stall regimes. This current effort is the first step towards multidisciplinary design optimization which involves coupling the aerodynamic heat transfer and structural analysis techniques. Applications include simulation of temporally deforming bodies and aeroelastic problems.

Boyalakuntla, Kishore↗

On a stochastic control method for weakly coupled linear systems

The stochastic control of two weakly coupled linear systems with different controllers is considered. Each controller only makes measurements about his own system; no information about the other system is assumed to be available. Based on the noisy measurements, the controllers are to generate independently suitable control policies which minimize a quadratic cost functional. To account for the effects of weak coupling directly, an approximate model, which involves replacing the influence of one system on the other by a white noise process is proposed. Simple suboptimal control problem for calculating the covariances of these noises is solved using the matrix minimum principle. The overall system performance based on this scheme is analyzed as a function of the degree of intersystem coupling.

Kwong, R. H.↗

How Certain Physical Considerations Impact Aerostructural Wing Optimization

Wing design optimization has been studied extensively and is of continued interest as optimization tools are developed and become more accessible. In each of these studies, certain assumptions and simplifications are made to make the design problem tractable. However, it is difficult to find systematic studies in which several considerations are added or removed one at a time to study how much impact they have. In this work, we examine how certain physical considerations (viscous drag, wave drag, thrust loads, and inertial relief from structural, fuel, and engine masses), impact the aerostructural optimization results for three distinct aircraft wings. The goal is to help develop a rough idea of how important these physical considerations are. We do this using gradient-based optimization and a multidisciplinary design optimization framework, OpenMDAO. We use the open-source tool OpenAeroStruct that couples a vortex lattice method to a finite element method. We establish a baseline aerostructural design optimization problem then perform a series of optimizations, each with one physical consideration removed from the baseline case. We find that depending on the size of the aircraft and flight conditions, the importance of some of these physical considerations varies considerably whereas the importance of others do not. Specifically, the optimal designs change radically without proper viscous and wave drag considerations and smaller aircraft with more distributed propulsion are more affected by the inclusion of engine loads.

Aeropropulsion↗

Comparison of approximate and numerical analyses of nonlinear combustion instability

At the present time, there are three general analytical techniques available to study problems of unsteady motions in rocket motors: linear stability analysis; approximate nonlinear analysis, founded on examining the behavior of coupled normal modes; and numerical calculations based on the conservation equations for one-dimensional flows. The last two yield the linear results as a limit. It is the main purpose of this paper to check the accuracy of the approximate analysis against the numerical analysis for some special cases. The results provide some justification for using the approximate analysis to study three-dimensional problems.

Culick, F. E. C.↗