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 595 records · Page 33

Intricacies of frustrated magnetism in the Kondo metal YbAgGe

The combination of localized magnetic moments, their frustration and interaction with itinerant electrons is a key challenge of condensed matter physics. Frustrated magnetic interactions promote degenerate ground states with enhanced fluctuations, a topic that is predominantly studied in magnetic insulators. The coupling between itinerant and localized electrons in metals add complexity to the problem, and is presently formulated only for extreme cases in which the itinerant electrons mediate exchange between localized spins (RKKY interaction) or suppress the formation of magnetic moments (Kondo screening). Here, we report an in-depth experimental study of the distorted Kagome metal YbAgGe, unravelling the open questions of how frustration, localized magnetism and itinerant electrons are intertwined in frustrated Kondo metals. We find that coupled itinerant and localized electrons give rise to dynamic magnetic correlations below T* ≈ 20 K. At lower temperature, frustrated magnetic interactions establish anisotropic magnetic short-range correlations that culminate into antiferromagnetic long-range order below T N = 0.68 K with a significantly reduced modulated magnetic moment. We show that local moment Hamiltonians can yield limited understanding of the microscopic behaviour in frustrated metals, and prompt the extension of more sophisticated model Hamiltonians incorporating itinerant effects.

Mazzone, Daniel G. [PSI Center for Neutron and Muo↗

Reliability and availability modeling of coupled communication networks - A simplified modeling approach

The network-complexity of LANs and of LANs that are interconnected by bridges and routers poses a challenging reliability-modeling problem. The present effort toward these problems' solution attempts to simplify them by reducing their number of states through truncation and state merging, as suggested by Shooman and Laemmel (1990). Through the use of state merging, it becomes possible to reduce the Bateman-Cortes 161 state model to a two state model with a closed-form solution. In the case of coupled networks, a technique which allows for problem-decomposition must be used.

Shooman, Martin L.↗

Stiffness coupling application to modal synthesis program, users guide

A FORTRAN IV computer program used to perform modal synthesis of structures by stiffness coupling, using the dynamic transformation method is described. The program was named SCAMP (Stiffness Coupling Approach Modal-Synthesis Program). The program begins with the entry of a substructure's physical mode shapes and eigenvalues or a substructure's mass and stiffness matrix. If the mass and stiffness matrices are entered, the eigen problem for the individual substructure is solved. Provisions are included for a maximum of 20 substructures which are coupled by stiffness matrix springs. Each substructure has a number degrees of freedom (DOF), except that for DOF greater than 100; vector sets having maximum row and column size of 100 were generated prior to entering SCAMP. The substructures are then coupled together via coupling springs, and the dynamic transformation is used to reduce the size of the eigen problem.

Kuhar, E. J.↗

A selective modal scheme for the analysis of EM coupling into or radiation from large open-ended waveguides

A selective modal scheme is proposed to efficiently analyze the problem of high-frequency (HF) electromagnetic (EM) coupling/penetration into or radiation from open-ended waveguides. This scheme is based on the phenomenon that at sufficiently high frequencies, the modes which contribute most significantly to the fields coupled into the waveguide are those whose modal ray directions are most nearly parallel to the incident-wave direction. This concept is illustrated by calculating the EM radiation and backscattering from open-ended parallel-plate, rectangular, circular, and sectoral waveguide geometries. The calculations use the usual geometrical optics, aperture field, and Ufimtsev edge current techniques. Also included are some measured results which further verify the accuracy of the above computations.

Altintas, Ayhan↗

Use of the NLP10x10 Sequential Quadratic Programming Algorithm To Solve Rotorcraft Hub Loads Minimisation Problems

Previous research and experimentation on the use of a non-linear programming constrained optimisation technique to define an optimal control vector for rotorcraft applications indicated that use of this methodology was feasible and desirable in many cases. In particular, use of non-linear programming methods that solve a sequence of related quadratic-programming sub-problems were used successfully to solve these problems. Accordingly, a licence for one of the latest versions of Professor Klaus Schittkowskis very successful Sequential Quadratic Programming NLPQLP software was obtained and used to experiment with and analyse typical optimisation problems of the type encountered in various rotorcraft wind tunnel and flight tests. This research resulted in the development of the general NLPQLP Computation System that could be used to solve problems of the type encountered in various rotorcraft applications where there is a linear dependence of the measurement vector on the control vector, and where equality andor inequality constraints might be imposed. This development was accomplished on a mainframe computer not part of actual wind tunnel andor flight-test experiment, but in a format which was transferable to wind tunnel lap-top computers. Emphasis was directed toward obtaining efficiency, robustness and speed in computation.The System was developed in support of the five-bladed SMART Rotor Active Flap Rotor Hub Loads analytical minimisation research. The design and development of the Computation System was tailored to address the particular requirements of the problem to minimise a performance metric function of measured hub load harmonic angular couple components by optimising the control vector harmonic flap angular couple components subject to constraints on the amplitudes of these control vector harmonic flap angular couple components. In addition, to facilitate real time wind tunnel experimentation, the ability to rapidly selectchange the particular hub load harmonic angular couple components andor the particular control vector harmonic angular couple components to be considered in the optimisation procedure was provided in the System. This capability allows the singling out of particular hub load frequencies andor particular flap angle frequencies to be analysed during testing operations. The System was used very successfully for the SMART Active Flap Rotor Hub minimisation problems considered in the study, the results of which were presented at the American Helicopter Society Fifth Decennial Aeromechanics Specialist Conference in January 2014. Excellent agreement between cases initiated with best guess starting estimates for the control vector elements and cases initiated with zero control vector starting element estimates resulted, indicating the robustness of the NLP10x10 algorithm.

Rotorcraft Hub Loads↗

Finite difference time domain calculations of antenna mutual coupling

The Finite Difference Time Domain (FDTD) technique has been applied to a wide variety of electromagnetic analysis problems, including shielding and scattering. However, the method has not been extensively applied to antennas. In this short paper calculations of self and mutual admittances between wire antennas are made using FDTD and compared with results obtained using the Method of Moments. The agreement is quite good, indicating the possibilities for FDTD application to antenna impedance and coupling.

Luebbers, Raymond J.↗

Venus tectonic styles and crustal differentiation

Two of the most important constraints are known from Pioneer Venus data: the lack of a system of spreading rises, indicating distributed deformation rather than plate tectonics; and the high gravity/topography ratio, indicating the absence of an asthenosphere. In addition, the high depth/diameter ratios of craters on Venus indicate that Venus probably has no more crust than Earth. The problems of the character of tectonics and crustal formation and recycling are closely coupled. Venus appears to lack a recycling mechanism as effective as subduction, but may also have a low rate of crustal differentiation because of a mantle convection pattern that is more distributed, less concentrated, than Earth's. Distributed convection, coupled with the nonlinear dependence of volcanism on heat flow, would lead to much less magmatism, despite only moderately less heat flow, compared to Earth. The plausible reason for this difference in convective style is the absence of water in the upper mantle of Venus. We have applied finite element modeling to problems of the interaction of mantle convection and crust on Venus. The main emphasis has been on the tectonic evolution of Ishtar Terra, as the consequence of convergent mantle flow. The early stage evolution is primarily mechanical, with crust being piled up on the down-stream side. Then the downflow migrates away from the center. In the later stages, after more than 100 m.y., thermal effects develop due to the insulating influence of the thickened crust. An important feature of this modeling is the entrainment of some crustal material in downflows. An important general theme in both convergent and divergent flows is that of mixing vs. stratification. Models of multicomponent solid-state flow obtain that lower-density crustal material can be entrained and recycled, provided that the ration of low-density to high-density material is small enough (as in subducted slabs on Earth). The same considerations should apply in upflows; a small percent of partial melt may be carried along with its matrix and never escape to the surface. Models that assume melt automatically rising to the crust and no entrainment or other mechanism of recycling lower-density material obtain oscillatory behavior, because it takes a long time for heat to build up enough to overcome a Mg-rich low-density residuum. However, these models develop much thicker crust than consistent with estimates from crater depth/diameter ratios.

Kaula, W. M.↗

Convergence Estimates for Multidisciplinary Analysis and Optimization

A quantitative analysis of coupling between systems of equations is introduced. This analysis is then applied to problems in multidisciplinary analysis, sensitivity, and optimization. For the sensitivity and optimization problems both multidisciplinary and single discipline feasibility schemes are considered. In all these cases a "convergence factor" is estimated in terms of the Jacobians and Hessians of the system, thus it can also be approximated by existing disciplinary analysis and optimization codes. The convergence factor is identified with the measure for the "coupling" between the disciplines in the system. Applications to algorithm development are discussed. Demonstration of the convergence estimates and numerical results are given for a system composed of two non-linear algebraic equations, and for a system composed of two PDEs modeling aeroelasticity.

Arian, Eyal↗

Coupled Marangoni-Benard/Rayleigh-Benard Instability with Temperature Dependent Viscosity

The onset of convection induced by coupled surface tension gradient and buoyancy forces is investigated with temperature dependent viscosity. Both surface tension and viscosity are assumed to vary linearly with temperature. The limiting case, Ma=0, is the buoyancy driven convection problem typically referred to as the Rayleigh-Benard stability problem. The other limiting case, Ra=0, is the surface tension (gradient) driven flow problem referred to as the Marangoni-Benard problem. The equations and boundary conditions obtained from the linear analysis are solved numerically as a generalized eigenvalue problem. Neutral stability curves for different viscosity slopes have been generated for the Marangoni-Benard and Rayleigh-Benard problems. It is shown that the curves can be collapsed to a single curve by appropriately scaling the results for each of the limiting cases. The critical Marangoni number is determined as a function of the slope of the viscosity temperature variation, epsilon, for different values of the Rayleigh number. When the viscosity decreases linearly with temperature, the coupled buoyancy-surface tension problem, including the limiting cases of Ra=0 and Ma=0, is found to be more stable than the constant viscosity case.

Skarda, J. Raymond Lee↗

Simulation of Acoustic Scattering from a Trailing Edge

Three model problems were examined to assess the difficulties involved in using a hybrid scheme coupling flow computation with the the Ffowcs Williams and Hawkings equation to predict noise generated by vortices passing over a sharp edge. The results indicate that the Ffowcs Williams and Hawkings equation correctly propagates the acoustic signals when provided with accurate flow information on the integration surface. The most difficult of the model problems investigated inviscid flow over a two-dimensional thin NACA airfoil with a blunt-body vortex generator positioned at 98 percent chord. Vortices rolled up downstream of the blunt body. The shed vortices possessed similarities to large coherent eddies in boundary layers. They interacted and occasionally paired as they convected past the sharp trailing edge of the airfoil. The calculations showed acoustic waves emanating from the airfoil trailing edge. Acoustic directivity and Mach number scaling are shown.

Singer, Bart A.↗

Concurrent two-way coupling of global and local models across internal boundaries with non-matching discretizations

Coupling local and global models enables efficient simulation of multiscale systems, where global models capture large-scale behavior and local models, with enhanced physics, resolve finer details over a smaller region. Here, this paper presents a mathematically consistent method for coupling physics-based models of varying fidelity across adjacent, non-overlapping subdomains, even when discretizations do not match at the immersed interdomain interfaces. Incompressible Navier-Stokes equations (NSE) constitute the global model while residual-based turbulence model serves as the local high-fidelity model. In addition, a scalar advection-diffusion equation that models the convection of an active scalar field is appended to the turbulence model in the local domain. This scalar field does not have its complement in the global model, giving rise to unequal number of equations at the immersed boundary between local and global models. Interdomain coupling terms are derived via the Variational Multiscale Discontinuous Galerkin (VMDG) method with new developments in scale representation and efficient fine-scale estimation. While transient laminar flows modeled with NSE in the global domain can be resolved with relatively coarse mesh, turbulent flow calculations in the local model require much finer spatial discretizations as well as smaller time-step for appropriately resolving the turbulent flow physics. The proposed framework also accommodates non-matching meshes at the immersed boundaries. Test problems in 2D and 3D numerically showcase the concurrent two-way coupling of unknown fields across the immersed boundaries. The 3D test presents a case with an unequal number of equations, where the scalar field represents the convection of contaminant concentration. This provides more detailed physics in the local region and highlights its application in climate modeling and atmospheric sciences.

Variational Multiscale Discontinuous Galerkin (VMD↗

A Galerkin method for the estimation of parameters in hybrid systems governing the vibration of flexible beams with tip bodies

An approximation scheme is developed for the identification of hybrid systems describing the transverse vibrations of flexible beams with attached tip bodies. In particular, problems involving the estimation of functional parameters are considered. The identification problem is formulated as a least squares fit to data subject to the coupled system of partial and ordinary differential equations describing the transverse displacement of the beam and the motion of the tip bodies respectively. A cubic spline-based Galerkin method applied to the state equations in weak form and the discretization of the admissible parameter space yield a sequence of approximating finite dimensional identification problems. It is shown that each of the approximating problems admits a solution and that from the resulting sequence of optimal solutions a convergent subsequence can be extracted, the limit of which is a solution to the original identification problem. The approximating identification problems can be solved using standard techniques and readily available software.

Banks, H. T.↗

Some problems in astroelasticity.

The author's studies in astroelasticity are reviewed, covering the dynamics of coupled systems of dynamic solids, the vibrations of a paraboloidal antenna, the thermal flutter of spacecraft boom and the effect of large structural deflection on control stability. A method proposed by Popov in nonlinear automatic control theory (1961) is used as the basis of the studies.

Yu, Y.-Y.↗

An efficient model for coupling structural vibrations with acoustic radiation

The scattering of an incident wave by a flexible panel is studied. The panel vibration is governed by the nonlinear plate equations while the loading on the panel, which is the pressure difference across the panel, depends on the reflected and transmitted waves. Two models are used to calculate this structural-acoustic interaction problem. One solves the three dimensional nonlinear Euler equations for the flow-field coupled with the plate equations (the fully coupled model). The second uses the linear wave equation for the acoustic field and expresses the load as a double integral involving the panel oscillation (the decoupled model). The panel oscillation governed by a system of integro-differential equations is solved numerically and the acoustic field is then defined by an explicit formula. Numerical results are obtained using the two models for linear and nonlinear panel vibrations. The predictions given by these two models are in good agreement but the computational time needed for the 'fully coupled model' is 60 times longer than that for 'the decoupled model'.

Frendi, Abdelkader↗

Linear stochastic optimal control and estimation problem

Problem involves design of controls for linear time-invariant system disturbed by white noise. Solution is Kalman filter coupled through set of optimal regulator gains to produce desired control signal. Key to solution is solving matrix Riccati differential equation. LSOCE effectively solves problem for wide range of practical applications. Program is written in FORTRAN IV for batch execution and has been implemented on IBM 360.

Geyser, L. C.↗

thornado+FLASH-X: A Hybrid Discontinuous Galerkin–Implicit-explicit and Finite-volume Framework for Neutrino-radiation Hydrodynamics in Core-collapse Supernovae

We present neutrino-transport algorithms implemented in the toolkit for high-order neutrino-radiation hydrodynamics (thornado) and their coupling to self-gravitating hydrodynamics within the adaptive mesh refinement–based multiphysics simulation framework FLASH-X. thornado, developed primarily for simulations of core-collapse supernovae (CCSNe), employs a spectral, six-species two-moment formulation with algebraic closure and special-relativistic observer corrections accurate to $\mathcal{O}(v/c)$, and uses discontinuous Galerkin (DG) methods for phase-space discretization combined with implicit-explicit time stepping. A key development is a nonlinear neutrino–matter coupling algorithm based on nested fixed-point iteration with Anderson acceleration, enabling fully implicit treatment of collisional processes, including energy-coupling interactions such as neutrino–electron scattering and pair production. Coupling to finite-volume (FV) hydrodynamics is achieved through a hybrid DG-FV representation of the fluid variables and operator-split evolution within FLASH-X. The implementation is verified using basic transport tests with idealized opacities and relaxation and deleptonization problems with tabulated microphysics. Spherically symmetric CCSN simulations demonstrate accuracy and robustness of the coupled scheme, including close agreement with the CCSN simulation code Chimera. An axisymmetric CCSN simulation further demonstrates the viability of DG-based neutrino transport for multidimensional supernova modeling within FLASH-X. thornado’s neutrino-transport solver is GPU-enabled using OpenMP offloading or OpenACC, and all CCSN applications included in this work use the GPU implementation. Together, these results establish a foundation for future enhancements in physics fidelity, numerical algorithms, and computational performance, for increasingly realistic large-scale CCSN simulations.

Endeve, Eirik [Oak Ridge National Laboratory (ORNL↗

The design, development, and flight testing of a modern-control-designed autoland system

This paper discusses the design, development, and flight test results of a full-state feedback modern-controls-designed autoland system - the Digital Integrated Automatic Landing System (DIALS). The lateral and longitudinal control laws were designed by formulating a linear quadratic regulator with disturbances problem. Although the designs were independent of each other, in implementation cross-coupling of the lateral and longitudinal variables does occur. A brief discussion of the control modes - localizer capture, localizer track, decrab, glideslope capture, glideslope track, and flare - and the modifications made to the basic design during the simulation phase to achieve desired performance is given. Some of the efforts required to implement the system in the flight computers and some problems encountered in the ground hardware simulation checkout are discussed. Finally, flight test data is presented for this system which performed 10 'hands off' automatic landings. The flight test data includes the performance of the aircraft in mild, gust, and wind shear conditions.

Hueschen, R. M.↗

Multidisciplinary Design Techniques Applied to Conceptual Aerospace Vehicle Design

Multidisciplinary design optimization (MDO) is an emerging discipline within aerospace engineering. Its goal is to bring structure and efficiency to the complex design process associated with advanced aerospace launch vehicles. Aerospace vehicles generally require input from a variety of traditional aerospace disciplines - aerodynamics, structures, performance, etc. As such, traditional optimization methods cannot always be applied. Several multidisciplinary techniques and methods were proposed as potentially applicable to this class of design problem. Among the candidate options are calculus-based (or gradient-based) optimization schemes and parametric schemes based on design of experiments theory. A brief overview of several applicable multidisciplinary design optimization methods is included. Methods from the calculus-based class and the parametric class are reviewed, but the research application reported focuses on methods from the parametric class. A vehicle of current interest was chosen as a test application for this research. The rocket-based combined-cycle (RBCC) single-stage-to-orbit (SSTO) launch vehicle combines elements of rocket and airbreathing propulsion in an attempt to produce an attractive option for launching medium sized payloads into low earth orbit. The RBCC SSTO presents a particularly difficult problem for traditional one-variable-at-a-time optimization methods because of the lack of an adequate experience base and the highly coupled nature of the design variables. MDO, however, with it's structured approach to design, is well suited to this problem. The result of the application of Taguchi methods, central composite designs, and response surface methods to the design optimization of the RBCC SSTO are presented. Attention is given to the aspect of Taguchi methods that attempts to locate a 'robust' design - that is, a design that is least sensitive to uncontrollable influences on the design. Near-optimum minimum dry weight solutions are determined for the vehicle. A summary and evaluation of the various parametric MDO methods employed in the research are included. Recommendations for additional research are provided.

Olds, John Robert↗