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 487 records · Page 27

Field Effects of Buoyancy on Lean Premixed Turbulent Flames

Buoyancy affects the entire flowfield of steady turbulent flames and this aspect of flame buoyancy coupling is largely unexplored by experiments or by theory. Open flames and flames within large confinements are free to expand and interact with the surrounding environment. In addition to fluid and combustion conditions, their aerodynamic flowfields are determined by the flame brush orientation and geometry, wake of the stabilizer, enclosure size, and of course, the gravitational field. Because the flowfield consists mainly of cold reactants (mostly in the nearfield) and hot products (mostly in the farfield), buoyancy effects are manifested in the farfield region. In upward pointing flames, an obvious effect is a favorable axial pressure gradient that accelerates the products thereby increasing the axial aerodynamic stretch rate. Intrinsic to turbulent flows, changes in mean aerodynamic stretch also couple to the fluctuating pressure field. Consequently, buoyancy can influence the turbulence intensities upstream and downstream of the flame. Flame wrinkling process, and heat release rate are also directly affected. This backward coupling mechanism is the so-called elliptic problem. To resolve the field effects of buoyancy would require the solution of three-dimensional non-linear Navier Stokes equations with full specification of the upstream, wall and downstream boundary conditions.

Cheng, R. K.↗

Three Program Architecture for Design Optimization

In this presentation, I would like to review historical perspective on the program architecture used to build design optimization capabilities based on mathematical programming and other numerical search techniques. It is rather straightforward to classify the program architecture in three categories as shown above. However, the relative importance of each of the three approaches has not been static, instead dynamically changing as the capabilities of available computational resource increases. For example, we considered that the direct coupling architecture would never be used for practical problems, but availability of such computer systems as multi-processor. In this presentation, I would like to review the roles of three architecture from historical as well as current and future perspective. There may also be some possibility for emergence of hybrid architecture. I hope to provide some seeds for active discussion where we are heading to in the very dynamic environment for high speed computing and communication.

Miura, Hirokazu↗

Control of Shock Loading from a Jet in the Presence of a Flexible Structure

The control of shock noise or screech from a jet near a flexible structure is discussed. The pressure from the supersonic jet consists of a shock with spiral and flapping nonaxisymmetric modes superimposed on broadband response. This shock induces a nonlinear-nonstationary loading problem associated with acoustic wave generation and propagation coupled with structural vibration. Control of the shock is achieved by placing a ring at the nozzle lip oscillating at the shock fundamental frequency. The ring prevents the shock characteristics originating in the column of the shear layer from sustaining connection with the out-of-phase surface vibration. Shock-free flow is maintained over a large pressure ratio. The peak power pressure level is reduced by 40 dB.

Maestrello, Lucio↗

Development of Probabilistic Life Prediction Methodologies and Testing Strategies for MEMS and CMC's

This effort is to investigate probabilistic life prediction methodologies for ceramic matrix composites and MicroElectroMechanical Systems (MEMS) and to analyze designs that determine stochastic properties of MEMS. For CMC's this includes a brief literature survey regarding lifing methodologies. Also of interest for MEMS is the design of a proper test for the Weibull size effect in thin film (bulge test) specimens. The Weibull size effect is a consequence of a stochastic strength response predicted from the Weibull distribution. Confirming that MEMS strength is controlled by the Weibull distribution will enable the development of a probabilistic design methodology for MEMS - similar to the GRC developed CARES/Life program for bulk ceramics. A main objective of this effort is to further develop and verify the ability of the Ceramics Analysis and Reliability Evaluation of Structures/Life (CARES/Life) code to predict the time-dependent reliability of MEMS structures subjected to multiple transient loads. A second set of objectives is to determine the applicability/suitability of the CARES/Life methodology for CMC analysis, what changes would be needed to the methodology and software, and if feasible, run a demonstration problem. Also important is an evaluation of CARES/Life coupled to the ANSYS Probabilistic Design System (PDS) and the potential of coupling transient reliability analysis to the ANSYS PDS.

Jadaan, Osama↗

Analyzing Aeroelastic Stability of a Tilt-Rotor Aircraft

Proprotor Aeroelastic Stability Analysis, now at version 4.5 (PASTA 4.5), is a FORTRAN computer program for analyzing the aeroelastic stability of a tiltrotor aircraft in the airplane mode of flight. The program employs a 10-degree- of-freedom (DOF), discrete-coordinate, linear mathematical model of a rotor with three or more blades and its drive system coupled to a 10-DOF modal model of an airframe. The user can select which DOFs are included in the analysis. Quasi-steady strip-theory aerodynamics is employed for the aerodynamic loads on the blades, a quasi-steady representation is employed for the aerodynamic loads acting on the vibrational modes of the airframe, and a stability-derivative approach is used for the aerodynamics associated with the rigid-body DOFs of the airframe. Blade parameters that vary with the blade collective pitch can be obtained by interpolation from a user-defined table. Stability is determined by examining the eigenvalues that are obtained by solving the coupled equations of motions as a matrix eigenvalue problem. Notwithstanding the relative simplicity of its mathematical foundation, PASTA 4.5 and its predecessors have played key roles in a number of engineering investigations over the years.

Kvaternil, Raymond G.↗

A Simple Data-Centric Methodology for Producible Geothermal Well Determinations: Preprint

The Bureau of Land Management (BLM) has traditionally lacked a standardized methodology for determining if a newly drilled geothermal well is "producible," a designation essential for deciding whether a lease should be "held by production." This is a straightforward problem to solve in oil and gas: Demonstrate that a well is economically viable, meaning it produces sufficient oil or gas to exceed direct operating costs and lease-related expenses, such as rentals or minimum royalties. In geothermal, the problem is more complex: Geothermal wells are tightly coupled with the downstream infrastructure - specifically, the power plant, which is often not designed until well after a lease is deemed as "held by production." Although this designation is critical for advancing geothermal power plant development on BLM-managed lands, current geothermal well assessments often rely on ad hoc approaches that can be complex, operator-biased, and heavy in assumptions related to economic viability. To address this, we have developed two complementary methodologies: a minimum power requirement-based approach and a productivity index (PI)-based approach. These methods leverage key flow test data - pressure, temperature, flow rate, and specific enthalpy - to provide reliable and standardized producible well determinations. The minimum power requirement-based approach evaluates wells against specific power output thresholds informed by reservoir experts and the associated temperature requirements. The PI-based approach assesses well productivity using widely accepted reservoir engineering metrics, proposing a threshold of 2.5 kg/s/bar. Both methods are data-driven and grounded in empirical production data from operational geothermal wells, avoiding uncertain economic assumptions while maintaining decision-making accuracy. Wells falling below key performance thresholds (i.e., PI, specific power) are deemed non-producible. These methodologies aim to streamline BLM's decision-making process, reduce nontechnical barriers to geothermal energy adoption, and enable regulatory expansion into states lacking geothermal expertise. Preliminary results indicate clear trends and thresholds in production data that provide actionable insights for evaluating well producibility. Validation using well completion report (WCR) data is ongoing, with promising results demonstrating the potential for these standardized methodologies to impact geothermal development significantly.

15 GEOTHERMAL ENERGY↗

Distributed optimization for multi-commodity urban traffic control

A distributed method for concurrent traffic signal and routing control of traffic networks is proposed. The method is based on the multi-commodity store-and-forward model, in which the destinations are the commodities. The system benefits from the communication between vehicles and infrastructure, providing optimal signal timings to intersections and routes to vehicles on a link-by-link basis. Using the augmented Lagrangian to model the constraints into the objective, the baseline centralized problem is decomposed into a set of objective-coupled subproblems, one for each intersection, enabling the solution to be computed by a distributed- gradient projection algorithm. Further, the intersection agents only need to communicate and coordinate with neighboring intersections to ensure convergence to the optimal solution while tolerating suboptimal iterations that offer more flexibility, unlike other distributed approaches. Through microsimulation, we demonstrate the effectiveness of the proposed algorithm in traffic networks with time-varying demand. Computational analysis shows that the distributed problem is suitable for real-time applications. A robustness analysis show that the distributed formulation enables a graceful degradation of the system in case of failure.

Augmented Lagrangian↗

A Nonlinear Schur Complement Solver for CFD-Based Multidisciplinary Models

CFD-based multidisciplinary models are the fundamental building blocks of multidis- ciplinary design optimization frameworks. Linear and nonlinear solutions of these coupled models are difficult, especially when the Jacobian matrices represent a saddle point problem, where a block-diagonal corresponding to a discipline is non-invertible. These scenarios necessitate the use of a coupled solver algorithm such as the Newton’s method instead of the popular block Gauss– Seidel-based methods because of this non-invertible block. To address this challenge, we introduce a nonlinear Schur complement solver suitable for CFD-based multidisciplinary models. The solver leverages the specialized linear and nonlinear solvers of the CFD code, and therefore, does not require the solution of a large coupled linear system as the coupled Newton’s method. Further- more, because the solver primarily uses the specialized linear and nonlinear solvers of the CFD code, it does not suffer from the same robustness limitations as the coupled Newton’s method. In this work, we will implement this solver in NASA’s OpenMDAO framework and demonstrate its effectiveness using a CFD-based aeropropulsive model. The solver will contribute to the develop- ment of aeropropulsive design optimization and CFD-based design optimization methods with the ultimate goal of accelerating the design and integration of advanced propulsion systems.

Nonlinear Solvers↗

A phase-field diffraction model for thermo-hydro-mechanical propagating fractures

This paper introduces a novel diffraction based thermo-hydraulic–mechanical (THM) model for fracture propagation using a phase-field fracture (PFF) approach. The key innovation of the THM-PFF model lies in its integrated treatment of four solution variables—displacements, phase-field, pressure, and temperature—each governed by a combination of conservation of momentum (mechanics problem), a variational inequality (constrained minimization problem), mass conservation (pressure problem), and energy conservation (temperature problem). This leads to a new formulation of a coupled variational inequality system. A major advancement is the development of an extended fixed-stress algorithm, where displacements, phase-field, pressures, and temperatures are solved in a staggered sequence. An important aspect of this work is the global coupling of pressures and temperatures across the domain using diffraction systems, with diffraction coefficients defined by material parameters weighted by the diffusive phase-field variable. To ensure robust local mass conservation, we employ enriched Galerkin finite elements (EG) for both pressure and temperature diffraction equations. By enriching the continuous Galerkin basis functions with discontinuous piecewise constants, EG accurately represents solution and parameter discontinuities while preserving local mass and energy conservation—crucial aspects for THM problems and realistic behavior. Moreover, the use of a predictor–corrector local mesh adaptivity scheme is employed, allowing the model to handle small phase-field length-scale parameters while maintaining high numerical accuracy and reasonable computational cost. Furthermore, these new model and algorithmic developments represent significant advances in the field and have been substantiated through rigorous numerical tests.

Diffraction systems↗

On Developing a Taxonomy for Multidisciplinary Design Optimization: A Decision-Based Perspective

In this paper, we approach MDO from a Decision-Based Design (DBD) perspective and explore classification schemes for designing complex systems and processes. Specifically, we focus on decisions, which are only a small portion of the Decision Support Problem (DSP) Technique, our implementation of DBD. We map coupled nonhierarchical and hierarchical representations from the DSP Technique into the Balling-Sobieski (B-S) framework (Balling and Sobieszczanski-Sobieski, 1994), and integrate domain-independent linguistic terms to complete our taxonomy. Application of DSPs to the design of complex, multidisciplinary systems include passenger aircraft, ships, damage tolerant structural and mechanical systems, and thermal energy systems. In this paper we show that Balling-Sobieski framework is consistent with that of the Decision Support Problem Technique through the use of linguistic entities to describe the same type of formulations. We show that the underlying linguistics of the solution approaches are the same and can be coalesced into a homogeneous framework with which to base the research, application, and technology MDO upon. We introduce, in the Balling-Sobieski framework, examples of multidisciplinary design, namely, aircraft, damage tolerant structural and mechanical systems, and thermal energy systems.

Lewis, Kemper↗

Integral dependent spin couplings in CI calculations

Although the number of ways to combine Slater determinants to form spin eigenfunctions increases rapidly with the number of open shells, most of these spin couplings will make only a small contribution to a given state, provided the spin coupling is chosen judiciously. The technique of limiting calculations to the interacting subspace pioneered by Bunge (1970) was employed by Munch and Davidson (1975) to the vanadium atom. The use of an interacting space looses its advantage in more complex cases. However, the problem can always be reduced to only one interacting spin coupling by making the coefficients integral dependent. The present investigation is concerned with the performance of integral dependent interacting couplings, taking into account the results of three test calculations.

Iberle, K.↗

Benchmarking Gas Path Diagnostic Methods: A Public Approach

Recent technology reviews have identified the need for objective assessments of engine health management (EHM) technology. The need is two-fold: technology developers require relevant data and problems to design and validate new algorithms and techniques while engine system integrators and operators need practical tools to direct development and then evaluate the effectiveness of proposed solutions. This paper presents a publicly available gas path diagnostic benchmark problem that has been developed by the Propulsion and Power Systems Panel of The Technical Cooperation Program (TTCP) to help address these needs. The problem is coded in MATLAB (The MathWorks, Inc.) and coupled with a non-linear turbofan engine simulation to produce "snap-shot" measurements, with relevant noise levels, as if collected from a fleet of engines over their lifetime of use. Each engine within the fleet will experience unique operating and deterioration profiles, and may encounter randomly occurring relevant gas path faults including sensor, actuator and component faults. The challenge to the EHM community is to develop gas path diagnostic algorithms to reliably perform fault detection and isolation. An example solution to the benchmark problem is provided along with associated evaluation metrics. A plan is presented to disseminate this benchmark problem to the engine health management technical community and invite technology solutions.

Simon, Donald L.↗

CELFE/NASTRAN Code for the Analysis of Structures Subjected to High Velocity Impact

CELFE (Coupled Eulerian Lagrangian Finite Element)/NASTRAN Code three-dimensional finite element code has the capability for analyzing of structures subjected to high velocity impact. The local response is predicted by CELFE and, for large problems, the far-field impact response is predicted by NASTRAN. The coupling of the CELFE code with NASTRAN (CELFE/NASTRAN code) and the application of the code to selected three-dimensional high velocity impact problems are described.

Chamis, C. C.↗

CELFE/NASTRAN Code for the analysis of structures subjected to high velocity impact

The CELFE (Coupled Eulerian Lagrangian Finite Element)/NASTRAN Code three-dimensional finite element code has the capability for analyzing of structures subjected to high velocity impact. The local response is predicted by CELFE and, for large problems, the far-field impact response is predicted by NASTRAN. The coupling of the CELFE code with NASTRAN (CELFE/NASTRAN code) and the application of the code to selected three-dimensional high velocity impact problems are described.

Chamis, C. C.↗

A CFD Approach to Modeling Spacecraft Fuel Slosh

Energy dissipation and resonant coupling from sloshing fuel in spacecraft fuel tanks is a problem that occurs in the design of many spacecraft. In the case of a spin stabilized spacecraft, this energy dissipation can cause a growth in the spacecrafts' nutation (wobble) that may lead to disastrous consequences for the mission. Even in non-spinning spacecraft, coupling between the spacecraft or upper stage flight control system and an unanticipated slosh resonance can result in catastrophe. By using a Computational Fluid Dynamics (CFD) solver such as Fluent, a model for this fuel slosh can be created. The accuracy of the model must be tested by comparing its results to an experimental test case. Such a model will allow for the variation of many different parameters such as fluid viscosity and gravitational field, yielding a deeper understanding of spacecraft slosh dynamics. In order to gain a better understanding of the dynamics behind sloshing fluids, the Launch Services Program (LSP) at the NASA Kennedy Space Center (KSC) is interested in finding ways to better model this behavior. Thanks to past research, a state-of-the-art fuel slosh research facility was designed and fabricated at Embry Riddle Aeronautical University (ERAU). This test facility has produced interesting results and a fairly reliable parameter estimation process to predict the necessary values that accurately characterize a mechanical pendulum analog model. The current study at ERAU uses a different approach to model the free surface sloshing of liquid in a spherical tank using Computational Fluid Dynamics (CFD) methods. Using a software package called Fluent, a model was created to simulate the sloshing motion of the propellant. This finite volume program uses a technique called the Volume of Fluid (VOF) method to model the interaction between two fluids [4]. For the case of free surface slosh, the two fluids are the propellant and air. As the fuel sloshes around in the tank, it naturally displaces the air. Using the conservation of mass, momentum, and energy equations, as well as the VOF equations, one can predict the behavior of the sloshing fluid and calculate the forces, pressure gradients, and velocity field for the entire liquid as a function of time.

Marsell, Brandon↗

A linear decomposition method for large optimization problems. Blueprint for development

A method is proposed for decomposing large optimization problems encountered in the design of engineering systems such as an aircraft into a number of smaller subproblems. The decomposition is achieved by organizing the problem and the subordinated subproblems in a tree hierarchy and optimizing each subsystem separately. Coupling of the subproblems is accounted for by subsequent optimization of the entire system based on sensitivities of the suboptimization problem solutions at each level of the tree to variables of the next higher level. A formalization of the procedure suitable for computer implementation is developed and the state of readiness of the implementation building blocks is reviewed showing that the ingredients for the development are on the shelf. The decomposition method is also shown to be compatible with the natural human organization of the design process of engineering systems. The method is also examined with respect to the trends in computer hardware and software progress to point out that its efficiency can be amplified by network computing using parallel processors.

Sobieszczanski-Sobieski, J.↗

Experimental Verification of the Use of Metal Filled Via Hole Fences for Crosstalk Control of Microstrip Lines in LTCC Packages

Coupling between microstrip lines in dense RF packages is a common problem that degrades circuit performance. Prior three-dimensional-finite element method (3-D-FEM) electromagnetic simulations have shown that metal filled via hole fences between two adjacent microstrip lines actually Increases coupling between the lines: however, if the top of the via posts are connected by a metal strip, coupling is reduced. In this paper, experimental verification of the 3-D-FEM simulations is demonstrated for commercially fabricated low temperature cofired ceramic (LTCC) packages. In addition, measured attenuation of microstrip lines surrounded by the shielding structures is presented and shows that shielding structures do not change the attenuation characteristics of the line.

Ponchak, George E.↗

On the Contribution of Turbulent Boundary Layers to the Noise Inside a Fuselage

The following report deals in preliminary fashion with the transmission through a fuselage of random noise generated on the fuselage skin by a turbulent boundary layer. The concept of attenuation is abandoned and instead the problem is formulated as a sequence of two linear couplings: the turbulent boundary layer fluctuations excite the fuselage skin in lateral vibrations and the skin vibrations induce sound inside the fuselage. The techniques used are those required to determine the response of linear systems to random forcing functions of several variables. A certain degree of idealization has been resorted to. Thus the boundary layer is assumed locally homogeneous, the fuselage skin is assumed flat, unlined and free from axial loads and the 'cabin' air is bounded only by the vibrating plate so that only outgoing waves are considered. Some of the details of the statistical description have been simplified in order to reveal the basic features of the problem. The results, strictly applicable only to the limiting case of thin boundary layers, show that the sound pressure intensity is proportional to the square of the free stream density, the square of cabin air density and inversely proportional to the first power of the damping constant and to the second power of the plate density. The dependence on free stream velocity and boundary layer thickness cannot be given in general without a detailed knowledge of the characteristics of the pressure fluctuations in the boundary layer (in particular the frequency spectrum). For a flat spectrum the noise intensity depends on the fifth power of the velocity and the first power of the boundary layer thickness. This suggests that boundary layer removal is probably not an economical means for decreasing cabin noise. In general, the analysis presented here only reduces the determination of cabin noise intensity to the measurement of the effect of any one of our variables (free stream velocity, boundary layer thickness, plate thickness or the characteristic velocity of propagation in the plate). The plate generates noise by vibrating in resonance over a wide range of frequencies and increasing the damping constant is consequently an effective method of decreasing noise generation. One of the main features of the results is that the relevant quantities upon which noise intensity depends are non-dimensional numbers in which boundary layer and plate properties enter as ratios. This is taken as an indication that in testing models of structures for boundary layer noise it is not sufficient to duplicate in the model the structural characteristics of the fuselage. One must match properly the characteristics of the exicitng pressure fluctuations to that of the structure.

Corcos, G. M.↗