Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “coupling constant”

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,243 records · Page 69

Aeroelastic modal characteristics of mistuned blade assemblies: Mode localization and loss of eigenstructure

An investigation of the effects of small mistuning on the aeroelastic modes of bladed disk assemblies with aerodynamic coupling between blades is presented. The cornerstone of the approach is the use and development of perturbation methods that exhibit the crucial role of the interblade coupling and yield general findings regarding mistuning effects. It is shown that blade assemblies with weak aerodynamic interblade coupling are highly sensitive to small blade mistuning, and that their dynamics is quantitatively altered in the following ways: the regular pattern that characterizes the root locus of the tuned aeroelastic eigenvalues in the complex plane is totally lost; the aeroelastic mode shapes becomes severely localized to only a few blades of the assembly and lose their constant interblade phase angle feature; and curve veering phenomena take place when the eigenvalues are plotted versus a mistuning parameter.

Pierre, Christophe↗

Aeroelastic modal characteristics of mistuned blade assemblies - Mode localization and loss of eigenstructure

An investigation of the effects of small mistuning on the aeroelastic modes of bladed-disk assemblies with aerodynamic coupling between blades is presented. The cornerstone of the approach is the use and development of perturbation methods that exhibit the crucial role of the interblade coupling and yield general findings regarding mistuning effects. It is shown that blade assemblies with weak aerodynamic interblade coupling are highly sensitive to small blade mistuning, and that their dynamics is qualitatively altered in the following ways: the regular pattern that characterizes the root locus of the tuned aeroelastic eigenvalues in the complex plane is totally lost; the aeroelastic mode shapes become severely localized to only a few blades of the assembly and lose their constant interblade phase angle feature; curve veering phenomena take place when the eigenvalues are plotted versus a mistuning parameter.

Pierre, Christophe↗

Millimeter-wave passive ultra-compact imaging technology for synthetic vision & mobile platforms

Substantial technical progress was made on all of the three high-risk subsystems of this program. The subsystems include dielectric antenna, G-band receiver, and electro-optic image processor. Progress is approximately on-schedule for both the receiver and the electro-optic processor development, while greater than anticipated challenges have been discovered in the dielectric antenna development. Much of the information in this report was covered in greater detail in the One-Year Review Meeting held at TTC on 22 February 1996. The performance goals of the dielectric antenna project are: Scan Angle -- 20 deg. desired; Loss -- 6 dB end to end (3 dB average); Frequency -- 206-218 GHz (6% bandwidth); Beam width -- 0.25 deg.; and Length -- 12 inches. The scan angle requirement was chosen to satisfy the needs of aircraft pilots. This requirement, coupled with the presently limited bandwidth processors (1 GHz state-of-the-art and 12 GHz in development in this program) forces the antenna to be dielectric (high scan angle air-filled waveguide-based antennas would be too lossy and their performance would vary too much as a function of frequency). A high dielectric constant (e.g., 10) was initially chosen for the dielectric material. This choice lead to the following fabrication challenges: total thickness variation (TTV) tolerance is 1 micrometer; coupler spacing tolerance is 1 micrometer; width tolerance is larger, but unknown, and the surfaces must have mirror finish. Also of importance is the difficulty in obtaining raw materials that satisfy the overall length requirement of 12 inches while simultaneously satisfying the above specifications.

Olsen, Randall↗

An Ultrasensitive Hot-Electron Bolometer for Low-Background SMM Applications

We are developing a hot-electron superconducting transition-edge sensor (TES) that is capable of counting THz photons and operates at T = 0.3K. The main driver for this work is moderate resolution spectroscopy (R approx. 1000) on the future space telescopes with cryogenically cooled (approx. 4 K) mirrors. The detectors for these telescopes must be background-limited with a noise equivalent power (NEP) approx. 10(exp -19)-10(exp -20) W/Hz(sup 1/2) over the range v = 0.3-10 THz. Above about 1 THz, the background photon arrival rate is expected to be approx. 10-100/s), and photon counting detectors may be preferable to an integrating type. We fabricated superconducting Ti nanosensors with a volume of approx. 3x10(exp -3) cubic microns on planar substrate and have measured the thermal conductance G to the thermal bath. A very low G = 4x10(exp -14) W/K, measured at 0.3 K, is due to the weak electron-phonon coupling in the material and the thermal isolation provided by superconducting Nb contacts. This low G corresponds to NEP(0.3K) = 3x10(exp -19) W/Hz(sup 1/2). This Hot-Electron Direct Detector (HEDD) is expected to have a sufficient energy resolution for detecting individual photons with v > 0.3 THz at 0.3 K. With the sensor time constant of a few microseconds, the dynamic range is approx. 50 dB.

bolometers↗

HYDRODYNAMIC ANALYSIS OF MACROALGAE LOCAL MODEL USING COMPUTATIONAL FLUID DYNAMICS

In this article, a local scale, fully nonlinear coupled fluid-structural interaction (FSI) sugar kelp model has been developed using a computational fluid dynamics (CFD) method. In this model, to be consistent with available experimental data, the sugar kelp is approximated as elongated rectangles with smoothed isosceles triangles at the ends and a single kelp model with one end fixed in a channel with constant current model is developed. Several different current speeds are simulated, and the resulting drag forces and calculated drag coefficients are validated by comparison with experimental data from the literature. In a previous study, a global scale model was developed using a computational structural dynamics (CSD) method to simulate macroalgae farming system and guide the system configuration design. In the global scale model, the hydrodynamic forces are calculated using Morison’s equation and the kinematics and dynamics of the sugar kelp are simplified and the group of kelps attached to the long line is modeled as a slender structure with the same length and an effective diameter such that the volumes are consistent with the real physical system. This simplified model matches the weight and buoyancy but adjusting the hydrodynamic properties when the general hydrodynamic coefficients are employed. Therefore, optimal hydrodynamic coefficients used in global scale model were determined to obtain the hydrodynamic force more accurately. The validated local scale model is then be applied to determine the hydrodynamic coefficients of the simplified sugar kelp model for global dynamic analysis.

hydrodynamic modeling, macroalgae, CFD model↗

Tensile Deformation of DED Graded IN718-René41 Superalloy at 650°C by in-situ Neutron Diffraction and Digital Image Correlation

A compositionally graded material comprised IN718 (718, γ'/γ" strengthened superalloy) and René41 (R41, a medium γ' strengthened superalloy) was fabricated using laser blown-powder directed energy deposition (DED/LB-M). Hot gas path components that require high and low temperature capabilities at different locations may benefit from such graded superalloy, potentially eliminating the use of welded or mechanically joined components and reducing cost. In the present work, microstructure and 650°C tensile behavior of 718-R41 graded superalloy in as-built and heat-treated states are investigated by Digital Image Correlation (DIC), in-situ neutron diffraction, coupled with conventional characterization techniques. The applied heat treatment generates static recrystallization and equiaxed grains in 718 rich compositions, while R41 rich compositions remain partially or un-recrystallized. Yield strengths of 718 and R41 sections in heat-treated state are comparable to wrought forms but the graded compositions show weakness. From neutron experiment, diffraction elastic constants and single crystal elastic constants are calculated, which show small difference between the as-built and heat-treated states and between the γ and γ' phases. Grain-level load transfer from the (220) grains to (200) grains is found to be composition dependent, and qualitatively agrees with DIC measured macroscopic yield strength. Within the (200) grains, the γ/γ' phases deform elastically until the γ phase yields and afterwards the γ' phase takes load from γ phase.

digital image correlation↗

Probing Axion-Like-Particles at the CERN Gamma Factory

We report the aim of the proposed CERN Gamma Factory is to produce ≈10 17 photons per second with energies up to 400 MeV. The photon beam intensity is expected to be a factor of $\mathcal{O}$(10 7 ) larger than that of the presently available photon beams in the MeV energy range. In this work, its potential to probe physics beyond the Standard Model is explored. In particular, searches for axion like particles (ALPs) with dominant couplings to photons are discussed and various production scenarios - fixed target, photon-photon collision, and conversion by a magnetic field - and detection schemes considered - via decay to photons or back-conversion. It is found that the Gamma Factory in a fixed target mode can probe ALPs with mass m a ≲ $\mathcal{O}$(100 MeV) and decay constants larger than 10 7 GeV, improving by an order of magnitude the discovery potential of previous beam dump experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Hydration of α-UO 3 following storage under controlled conditions of temperature and relative humidity

Changes in chemical speciation of uranium oxides following storage under varied conditions of temperature and relative humidity are valuable for characterizing material provenance. In this study, subsamples of high purity α-UO 3 were stored under four sets of controlled conditions of temperature and relative humidity over several years, and then measured periodically for chemical speciation. Powder X-ray diffraction (XRD) analysis and extended X-ray absorption fine structure spectroscopy confirm hydration of α-UO 3 to a schoepite-like end product following storage under each of the varied storage conditions, but the species formed during exposure to the lower relative humidity and lower temperature condition follows different trends from those formed under the other three storage conditions (high relative humidity with high or low temperatures, and low relative humidity with a high temperature). Thermogravimetry coupled with XRD analysis was carried out to distinguish desorption pathways of water from the hydrated end products. Density functional theory calculations discern changes in the structure of α-UO 3 following incorporation of 1, 2 or 3 H 2 O molecules or 1, 2 or 3 OH groups into the orthorhombic lattice, revealing differences in lattice constants, U–O bond lengths, and U–U distances. The collective results from this analysis are in contrast to analogous studies that report that U 3 O 8 is oxidized and hydrated in air during storage under high relative humidity conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Universal relations for neutron star 𝑓-mode and 𝑔-mode oscillations

Among the various oscillation modes of neutron stars, 𝑓- and 𝑔- modes are the most likely to be ultimately observed in binary neutron star mergers due to their relatively large coupling and shared frequencies with tidal excitations. The 𝑓-mode frequency and damping time are known to correlate in normal neutron stars with their compactness, and previous fits to hadronic stars are extended and shown to be valid for an extremely broad sampling of equations of state using a piecewise polytropic parametrization scheme for hadrons and a constant sound-speed parametrization for quark matter. Separate fits applicable to quark (self-bound) stars are improved. Much more significant correlations exist with tidal deformability, and therefore with moment of inertia and quadrupole moment. It is conclusively demonstrated that these correlations are the same for all types of stars, whether hadronic, hybrid, or pure quark, and its accuracy is quantified. A novel 1-node branch of the 𝑓-mode that occurs in low-mass hybrid stars in a narrow mass range just beyond the critical mass necessary for a phase transition to appear is identified. This 1-node branch shows the largest, but still small, deviations from the universal correlation for any configuration. It is characterized by a nonmonotonic relation between neutron star mass and 𝑓-mode frequency, in contrast to the behavior otherwise observed in normal, quark and hybrid stars. The 𝑔-mode only exists in matter with a nonbarotropic equation of state involving temperature, chemical potential or composition (such as being out of beta equilibrium), or a phase transition in barotropic matter. Here, the 𝑔-mode therefore could serve as a probe for studying phase transitions in hybrid stars. In contrast with the 𝑓-mode, 𝑔-mode frequencies do not correlate well with tidal deformability, but depend strongly on properties of the transition (the density and the magnitude of the discontinuity) at the transition. Imposing causality and maximum mass constraints, a fit involving neutron star and phase transition properties is found and the 𝑔-mode frequency is determined to have an upper bound of about 1.25 kHz. However, if the sound speed 𝑐𝑠 in the inner core at densities above the phase transition density is restricted to 𝑐$^{2}_{𝑠}$ ≤1/3, 𝑔-mode frequencies can only reach about 0.8 kHz, which are significantly lower than 𝑓-mode frequencies (1.3–2.8 kHz). 𝑔-mode gravitational wave damping times are found to be extremely long, >10 4 s (102 s) in the inner core with 𝑐$^{2}_{𝑠}$ ≤1/3 (1), in comparison with 𝑓-mode damping times (0.1–1 s).

Composition of astronomical objects↗

A helicopter flight investigation of roll-control sensitivity, damping and cross coupling in a low altitude lateral maneuvering task

A helicopter in-flight simulation was conducted to determine the effects of variations in roll damping, roll sensitivity, and pitch and roll rate cross-coupling on helicopter flying qualities in a low altitude maneuver. The experiment utilized the UH-1H helicopter in-flight simulator, which is equipped with the V/STOLAND avionics system. The response envelope of this vehicle allowed simulation of configurations with low to moderate damping and sensitivity. A visual, low level slalom course was set up, consisting of constant speed and constant altitude S-turns around the 1000 ft makers of an 8000 ft runway. Results are shown in terms of Cooper-Harper pilot ratings, pilot commentary, and statistical and frequency analyses of the lateral characteristics. These results show good consistency with previous ground simulator results and are compared with existing flying qualities criteria.

Corliss, L. D.↗

Experimental Observations on a Low Strain Counter-Flow Diffusion Flame: Flow and Bouyancy Effects

Diffusion flames are of great interest in fire safety and many industrial processes. The counter-flow configuration provides a constant strain flow, and therefore is ideal to study the structure of diffusion flames. Most studies have concentrated on the high velocity, high strain limit, since buoyantly induced instabilities will disintegrate the planar flame as the velocity decreases. Only recently, experimental studies in microgravity conditions have begun to explore the low strain regimes. Numerical work has shown the coupling between gas phase reaction rates, soot reaction rates, and radiation. For these programs, size, geometry and experimental conditions have been chosen to keep the flame unaffected by the physical boundaries. When the physical boundaries can not be considered infinitely far from the reaction zone discrepancies arise. A computational study that includes boundary effects and accounts for the deviations occurring when the major potential flow assumptions are relaxed was presented by Borlik et al. This development properly incorporates all heat loss terms and shows the possibility of extinction in the low strain regime. A major constraint of studying the low strain regime is buoyancy. Buoyant instabilities have been shown to have a significant effect on the nature of reactants and heat transport, and can introduce instabilities on the flow that result in phenomena such as flickering or fingering. The counter-flow configuration has been shown to provide a flame with no symmetry disrupting instabilities for inlet velocities greater than 50 mm/s. As the velocity approaches this limit, the characteristic length of the experiment has to be reduced to a few millimetres so as to keep the Rayleigh number (Ra L = (β g 0 L 3 ∆T)/(𝑎v)) below 2000. In this work, a rectangular counter-flow burner was used to study a two-dimensional counter-flow diffusion flame. Flow visualisation and Particle Image Velocimetry served to describe the nature of the stagnation plane for strain rates smaller than 100 (1/s). These experiments were conducted with a non-reacting flow. Video images of a propane air diffusion flame were used to describe the behaviour of a diffusion flame in this regime. Flame geometry and pulsation frequency are described.

J. A. Sutula↗

Investigation of Conjugate Heat Transfer in Turbine Blades and Vanes

We report on work carried out to develop a 3-D coupled Finite Volume/BEM-based temperature forward/flux back (TFFB) coupling algorithm to solve the conjugate heat transfer (CHT) which arises naturally in analysis of systems exposed to a convective environment. Here, heat conduction within a structure is coupled to heat transfer to the external fluid which is convecting heat into or out of the solid structure. There are two basic approaches to solving coupled fluid structural systems. The first is a direct coupling where the solution of the different fields is solved simultaneously in one large set of equations. The second approach is a loose coupling strategy where each set of field equations is solved to provide boundary conditions for the other. The equations are solved in turn until an iterated convergence criterion is met at the fluid-solid interface. The loose coupling strategy is particularly attractive when coupling auxiliary field equations to computational fluid dynamics codes. We adopt the latter method in which the BEM is used to solve heat conduction inside a structure which is exposed to a convective field which in turn is resolved by solving the NASA Glenn compressible Navier-Stokes finite volume code Glenn-HT. The BEM code features constant and bi-linear discontinuous elements and an ILU-preconditioned GMRES iterative solver for the resulting non-symmetric algebraic set arising in the conduction solution. Interface of flux and temperature is enforced at the solid/fluid interface, and a radial-basis function scheme is used to interpolated information between the CFD and BEM surface grids. Additionally, relaxation is implemented in passing the fluxes from the conduction solution to the fluid solution. Results from a simple test example are reported.

Kassab, A. J.↗

Analysis of Advanced Thermoelectric Materials and Their Functional Limits

The world's demand for energy is increasing dramatically, but the best energy conversion systems operate at approximately 30% efficiency. One way to decrease energy loss is in the recovery of waste heat using thermoelectric (TE) generators. A TE generator is device that generates electricity by exploiting heat flow across a thermal gradient. The efficiency of a TE material for power generation and cooling is determined by the dimensionless Figure of Merit (ZT): ZT = S(exp. 2)sigmaT/κ: where S is the Seebeck coefficient, sigma is the electrical conductivity, T is the absolute temperature, and κ is the thermal conductivity. The parameters are not physically independent, but intrinsically coupled since they are a function of the transport properties of electrons. Traditional research on TE materials has focused on synthesizing bulk semiconductor-type materials that have low thermal conductivity and high electrical conductivity affording ZT values of 1. The optimization of the σ/κ ratio is difficult to achieve using current material formats, as these material constants are complementary. Recent areas of research are focusing on using nanostructural artifacts that introduce specific dislocations and boundary conditions that scatter the phonons. This disrupts the physical link between thermal (phonon) and electrical (electron) transport. The result is that κ is decreased without decreasing σ. These material formats give ZT values of up to 2 which represent approximately 18% energy gain from waste heat recovery. The next challenge in developing the next generation of TE materials with superior performance is to tailor the interconnected thermoelectric physical parameters of the material system. In order to approach this problem, the fundamental physics of each parameter S, sigma, and κ need to be physically understood in their context of electron/phonon interaction for the construction of new high ZT thermoelectric devices. Is it possible to overcome the physical limit imposed by of the effect of phonon lattice oscillation and energetic electrons towards thermal conductivity? Is the Seebeck coefficient, based on the difference in voltage over temperature gradient ( deltaV/deltaT), an intrinsic parameter of each material? All these parameters were manipulated using nano-bridge and twin-lattice structural concepts at the NASA Langley Research Center. This talk will review the current trend of TE research to optimize the ZT and discuss about new approaches on increasing ZT within functional limits of each parameter.

Kim, Hyun Jung↗

Symmetry-driven persistent spin texture for the two-dimensional nonsymmorphic CdTe and ZnTe crystal structures

In this paper, two nonsymmorphic two-dimensional structures of CdTe and ZnTe are modeled, and using state-of-the-art density functional theory with the group theory of solids, their symmetry-enforced electronic properties are studied. The in-plane ferroelectricity coupled with strong spin-orbit coupling induces a unidirectional out-of-plane Rashba spin-orbit field that can host a momentum-independent uniform spin configuration known as persistent spin texture (PST) at the Brillouin zone center. PST in these structures is found to be robust against external perturbations such as strain, structural distortion, and independent of layer thickness. These unprecedented intrinsic spin transport properties hold utmost importance in spintronics, as the experimental stringent condition of equal Rashba and Dresselhaus constants [Phys. Rev. Lett. 90, 146801 (2003)] is eliminated. The calculated persistent spin helix wavelength of <~5 nm paves the way for developing next-generation nanosized nonballistic spin field-effect transistors compared with micrometer-sized GaAs/AlGaAs quantum wells. Further, these materials exhibit finite spin Hall conductivity at the band edges and hence can be used in ferromagnet-free spin Hall transistors. Although CdTe and ZnTe systems have been widely studied for photocatalysis and solar cell applications over the past few decades, their potential application in spintronic devices has not been explored. Mono/few layers of CdTe and ZnTe synthesized from (110) facets of bulk zinc-blende crystals [Nat. Commun. 3, 1057 (2012)] satisfy all symmetry operations of the nonsymmorphic space group and hence can be considered ideal materials to verify our theoretical results experimentally.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Pulsating Hydrodynamic Instability and Thermal Coupling in an Extended Landau/Levich Model of Liquid-Propellant Combustion

Hydrodynamic (Landau) instability in combustion is typically associated with the onset of wrinkling of a flame surface, corresponding to the formation of steady cellular structures as the stability threshold is crossed. In the context of liquid-propellant combustion, such instability has recently been shown to occur for critical values of the pressure sensitivity of the burning rate and the disturbance wavenumber, significantly generalizing previous classical results for this problem that assumed a constant normal burning rate. Additionally, however, a pulsating form of hydrodynamic instability has been shown to occur as well, corresponding to the onset of temporal oscillations in the location of the liquid/gas interface. In the present work, we consider the realistic influence of a nonzero temperature sensitivity in the local burning rate on both types of stability thresholds. It is found that for sufficiently small values of this parameter, there exists a stable range of pressure sensitivities for steady, planar burning such that the classical cellular form of hydrodynamic instability and the more recent pulsating form of hydrodynamic instability can each occur as the corresponding stability threshold is crossed. For larger thermal sensitivities, however, the pulsating stability boundary evolves into a C-shaped curve in the disturbance-wavenumber/ pressure-sensitivity plane, indicating loss of stability to pulsating perturbations for all sufficiently large disturbance wavelengths. It is thus concluded, based on characteristic parameter values, that an equally likely form of hydrodynamic instability in liquid-propellant combustion is of a nonsteady, long-wave nature, distinct from the steady, cellular form originally predicted by Landau.

Margolis, S. B.↗

Assessment of diffuse-interface methods for compressible multiphase fluid flows and elastic-plastic deformation in solids

This work describes three diffuse-interface methods for the simulation of immiscible, compressible multiphase fluid flows and elastic-plastic deformation in solids. The first method is the localized-artificial-diffusivity approach of Cook, Subramaniam et al., and Adler and Lele, in which artificial diffusion terms are added to the individual phase mass fraction transport equations and are coupled with the other conservation equations. The second method is the gradient-form approach that is based on the quasi-conservative method of Shukla et al., in which the diffusion and sharpening terms (together called regularization terms) are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. The third approach is the divergence-form approach that is based on the fully conservative method of Jain et al., in which the regularization terms are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. In the present study, all three diffuse-interface methods are used in conjunction with a four-equation, multicomponent mixture model, in which pressure and temperature equilibria are assumed among the various phases. The primary objective of this work is to compare these three methods in terms of their ability to: maintain constant interface thickness throughout the simulation; conserve mass, momentum, and energy; and maintain accurate interface shape for long-time integration. The second objective of this work is to consistently extend these methods to model interfaces between solid materials with strength. To assess and compare the methods, they are used to simulate a wide variety of problems, including (1) advection of an air bubble in water, (2) shock interaction with a helium bubble in air, (3) shock interaction and the collapse of an air bubble in water, and (4) Richtmyer–Meshkov instability of a copper–aluminum interface. The current work focuses on comparing these methods in the limit of relatively coarse grid resolution, which illustrates the true performance of these methods. In conclusion, this is because it is rarely practical to use hundreds of grid points to resolve a single bubble or drop in large-scale simulations of engineering interest.

97 MATHEMATICS AND COMPUTING↗

Oceanic lithosphere and asthenosphere: The thermal and mechanical structure

A coupled thermal and mechanical solid state model of the oceanic lithosphere and asthenosphere is presented. The model includes vertical conduction of heat with a temperature dependent thermal conductivity, horizontal and vertical advection of heat, viscous dissipation or shear heating, and linear or nonlinear deformation mechanisms with temperature and pressure dependent constitutive relations between shear stress and strain rate. A constant horizontal velocity u sub 0 and temperature t sub 0 at the surface and zero horizontal velocity and constant temperature t sub infinity at great depth are required. In addition to numerical values of the thermal and mechanical properties of the medium, only the values of u sub 0, t sub 0 and t sub infinity are specified. The model determines the depth and age dependent temperature horizontal and vertical velocity, and viscosity structures of the lithosphere and asthenosphere. In particular, ocean floor topography, oceanic heat flow, and lithosphere thickness are deduced as functions of the age of the ocean floor.

Schubert, G.↗

A Robust Control Design Framework for Substructure Models

A framework for designing control systems directly from substructure models and uncertainties is proposed. The technique is based on combining a set of substructure robust control problems by an interface stiffness matrix which appears as a constant gain feedback. Variations of uncertainties in the interface stiffness are treated as a parametric uncertainty. It is shown that multivariable robust control can be applied to generate centralized or decentralized controllers that guarantee performance with respect to uncertainties in the interface stiffness, reduced component modes and external disturbances. The technique is particularly suited for large, complex, and weakly coupled flexible structures.

Lim, Kyong B.↗