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 379 records · Page 21

Prediction and Analysis of Utah FORGE Injection Activities using a Coupled Thermo-hydro-mechanical and Earthquake (THM+E) Modeling Workflow

A coupled thermo-hydro-mechanical (THM) numerical workflow that is capable of modeling seismic slip is critical for the successful development of enhanced geothermal systems (EGS). By integrating key physical processes, this workflow enables accurate simulation of temperature and pressure diffusions, stress changes, and induced seismicity. As a result, it serves as a vital tool for predicting induced seismicity and optimizing reservoir stimulation strategies. The Utah FORGE (Frontier Observatory for Research in Geothermal Energy) project, located near Milford, Utah, is a U.S. Department of Energy initiative aimed at advancing EGS technology. In April 2024, eight new stimulation stages (Stages 3R-10) were conducted in well 16A (injection well) subsequent to the first series of stimulation (Stages 1-3) performed in April, 2022. To monitor the induced seismicity, geophones were deployed in wells 58-32, 56-32, and 78B-32, while fiber optic cables were also installed in wells 16B, 78-32, and 78B-32 to collect microseismic data and detect frac hits Preliminary analyses of microseismic catalogs and fiber optic data suggest that the stimulated fractures in Stages 3R–6 closely align with that generated during Stage 3, indicating that the new stimulations were likely reactivating the previously stimulated fracture. To better understand the underlying process, a comprehensive modeling approach that can accurately capture thermal, hydrological, mechanical, and seismic responses is essential. In this work, we propose and utilize a coupled thermo-hydro-mechanical and earthquake (THM+E) simulation workflow to numerically investigate the stimulation activities on well 16A. The specific objective is to confirm whether the new stimulation stages (Stages 3R–6) reactivated fractures previously stimulated during Stage 3. For this purpose, we perform THM+E simulations individually for Stages 3, 3R, 4, and 5, incorporating the discrete fracture networks (DFNs) created by the plane-fitting technique based on the microseismic catalogs. The simulation workflow consists of two separate models: a THM model and an earthquake model, coupled in a one-way manner. Detailed descriptions of the workflow are provided in Section 3. Simulation results are presented in terms of injection pressure, permeability evolution, and predicted seismic catalogs, which are then compared with field data for further analyses. This report is structured as follows. In Section 2, we present detailed analyses of the field data and propose the hypothesis that the new stimulation stages (Stages 3R–6) were probably reactivating the previously stimulated fractures in Stage 3. In Section 3, we introduce the coupled THM+E workflow and the problem setup to validate our hypothesis, followed by the simulation results for each stage in Section 4. Meanwhile, discussions are included to analyze the model predictions and their comparison with field data. Lastly, we conclude the report and outline future plans in Section 5.

15 GEOTHERMAL ENERGY↗

Repeated cascade theory of homogeneous turbulence.

The problem of turbulent spectrum engenders two coupled hierarchies: one originates from the development of stress, leading to a transfer function, and the other from the development of an eddy viscosity. In order to incorporate physical roles among scales, the turbulent velocity fluctuation is decomposed into a series of ranks in the increasing order of randomness, contributing successively to energy or stress, eddy viscosity, relaxation frequency, and higher-rank frequencies in the memory chain. As a result, the first hierarchy mentioned above becomes closed at the quadrupole correlation. The second hierarchy governs the eddy viscosities of different ranks, related to relaxation frequencies of such ranks, in the form of a memory chain. It is cut off by an implicit viscous mechanism. For zero wind gradient, the spectrum in the inertial subrange recovers the Kolmogoroff k to the minus 5/3 law with a numerical constant 1.58, in good agreement with experiments. For a strong wind gradient, the spectrum in the production subrange has a k to the minus 1 law.

Tchen, C. M.↗

A direct method for calculating flutter speeds

The basic objective of this work is to provide a direct and simple means for calculating classical flutter. The analysis makes no use of free vibration modes, although the method is readily adaptable to a modal type approach. Only subsonic aerodynamics are considered; however, the method applies to any speed regime. The problems of following modes and modal coupling are completely avoided. The simplicity and accuracy of the direct solution is stressed. The analysis is applicable to both straight and swept wings although only straight wings, uniform as well as nonuniform, are presented.

Beres, D. P.↗

A nonlinear theory of cosmic ray pitch angle diffusion in homogeneous magnetostatic turbulence

A plasma strong turbulence, weak coupling, theory is applied to the problem of cosmic ray pitch angle scattering in magnetostatic turbulence. The theory used is a rigorous generalization of Weinstock's resonance-broadening theory and contains no ad hoc approximations. A detailed calculation is presented for a model of slab turbulence with an exponential correlation function. The results agree well with numerical simulations. The rigidity dependence of the pitch angle scattering coefficient differs from that found by previous researchers. The differences result from an inadequate treatment of particle trajectories near 90 deg pitch angle in earlier work.

Goldstein, M. L.↗

Fiber distributed feedback laser

Utilizing round optical fibers as communication channels in optical communication networks presents the problem of obtaining a high efficiency coupling between the optical fiber and the laser. A laser is made an integral part of the optical fiber channel by either diffusing active material into the optical fiber or surrounding the optical fiber with the active material. Oscillation within the active medium to produce lasing action is established by grating the optical fiber so that distributed feedback occurs.

Elachi, C.↗

A nonlinear theory of cosmic-ray pitch-angle diffusion in homogeneous magnetostatic turbulence

A plasma strong turbulence, weak coupling, theory is applied to the problem of cosmic-ray pitch-angle scattering in magnetostatic turbulence. The theory used is a rigorous generalization of Weinstock's 'resonance broadening' theory and contains no ad hoc approximations. A detailed calculation is presented for a model of 'slab' turbulence with an exponential correlation function. The results agree well with numerical simulations. The rigidity dependence of the pitch-angle scattering coefficient differs from that found by previous researchers. The differences result from an inadequate treatment of particle trajectories near 90 deg pitch angle in earlier work.

Goldstein, M. L.↗

Initial values for the integration scheme to compute the eigenvalues for propagation in ducts

A scheme for the calculation of eigenvalues in the problem of acoustic propagation in a two-dimensional duct is described. The computation method involves changing the coupled transcendental nonlinear algebraic equations into an initial value problem involving a nonlinear ordinary differential equation. The simplest approach is to use as initial values the hardwall eigenvalues and to integrate away from these values as the admittance varies from zero to its actual value with a linear variation. The approach leads to a powerful root finding routine capable of computing the transverse and axial wave numbers for two-dimensional ducts for any frequency, lining, admittance and Mach number without requiring initial guesses or starting points.

Eversman, W.↗

Engineering description of the OMS/RCS/DAP modes used in the HP-9825A High Fidelity Relative Motion Program (HFRMP)

Simplified mathematical models are reported for the space shuttle's Orbital Maneuvering System (OMS), Reaction Control System (RCS), and on-orbit Digital Autopilot (DAP) that have been incorporated in the High-Fidelity Relative Motion Program (HFRMP) for the HP-9825A desk-top calculator. Comparisons were made between data generated by the HFRMP and by the Space Shuttle Functional Simulator (SSFS), which models the cited shuttle systems in much greater detail. These data include propellant requirements for representative translational maneuvers, rotational maneuvers, and attitude maintenance options. Also included are data relating to on-orbit trajectory deviations induced by RCS translational cross coupling. Potential close-range stationkeeping problems that are suggested by HFRMP simulations of 80 millisecond (as opposed to 40 millisecond) DAP cycle effects are described. The principal function of the HFRMP is to serve as a flight design tool in the area of proximity operations.

Wilson, S. W.↗

Sudden stretching of a four layered composite plate

An approximate theory of laminated plates is developed by assuming that the extensioral and thickness mode of vibration are coupled. The mixed boundary value crack problem of a four layered composite plate is solved. Dynamic stress intensity factors for a crack subjected to suddenly applied stress are found to vary as a function of time and depend on the material properties of the laminate. Stress intensification in the region near the crack front can be reduced by having the shear modulus of the inner layers to be larger than that of the outer layers.

Sih, G. C.↗

Formulation of coupled rotor/fuselage equations of motion

The governing equations of motion of a helicopter rotor coupled to a rigid body fuselage are derived. A consistent formulation is used to derive nonlinear periodic coefficient equations of motion which can be used to study coupled rotor/fuselage dynamics in forward flight. The methodology of rotor/fuselage coupling is clearly described and the importance of an ordering scheme in deriving consistent nonlinear equations of motion is reviewed. The final equations which are presented in partial differential form can be used to model coupled rotor/fuselage aeroelastic response or stability problems.

Warmbrodt, W.↗

The quantum limit for gravitational-wave detectors and methods of circumventing it

The Heisenberg uncertainty principle prevents the monitoring of the complex amplitude of a mechanical oscillator more accurately than a certain limit value. This 'quantum limit' is a serious obstacle to the achievement of a 10 to the -21st gravitational-wave detection sensitivity. This paper examines the principles of the back-action evasion technique and finds that this technique may be able to overcome the problem of the quantum limit. Back-action evasion does not solve, however, other problems of detection, such as weak coupling, large amplifier noise, and large Nyquist noise.

Thorne, K. S.↗

Analysis of the Earth's variable rotation

The objective was to analyze the Earth's rotation, and to compare the observations with the atmospheric and seismic excitation functions. The variation in the Earth's rotation can be separated into the length-of-day (LOD) variation and the polar motion. Possible driving mechanisms include: atmospheric/oceanic circulations, seismic activities, solar-lunar tides, mantle convection, core-mantle coupling, and solar activities. The major problem of concern was the identification of these dynamical processes as primary driving mechanisms for the various features in the variation of the Earth's rotation.

Chao, B. F.↗

Buckling and vibration of any prismatic assembly of shear and compression loaded anisotropic plates with an arbitrary supporting structure

The computer program designated 'VIPASA', which accurately treats buckling and vibration in prismatic plate assemblies with a response that varies sinusoidally in the longitudinal direction, has been found to be limited by the production of an in-plane shear loading of component plates that produces skewed mode shapes. These do not conform to desired support conditions. This problem is presently overcome through a coupling of the VIPASA stiffness matrices for different wavelength responses by means of the Lagrangian Multipliers method. The theory extends to supports at arbitrary locations, and even to the support provided by any elastic structure. The generality and capabilities of VIPASA have been retained in the computer program designated 'VICON', which permits constraints and a supporting structure consisting of any number of transverse beam columns.

Anderson, M. S.↗

Physical processes in the strong magnetic fields of accreting neutron stars

Analytical formulae are fitted to observational data on physical processes occurring in strong magnetic fields surrounding accreting neutron stars. The propagation of normal modes in the presence of a quantizing magnetic field is discussed in terms of a wave equation in Fourier space, quantum electrodynamic effects, polarization and mode ellipticity. The results are applied to calculating the Thomson scattering, bremsstrahlung and Compton scattering cross-sections, which are a function of the frequency, angle and polarization of the magnetic field. Numerical procedures are explored for solving the radiative transfer equations. When applied to modeling X ray pulsars, a problem arises in the necessity to couple the magnetic angle and frequency dependence of the cross-sections with the hydrodynamic equations. The use of time-dependent averaging and approximation techniques is indicated.

Meszaros, P.↗

Development of an integrated BEM approach for hot fluid structure interaction

The development of a boundary element formulation for the study of hot fluid-structure interaction in earth-to-orbit engine hot section components is described. The initial primary thrust of the program to date was directed quite naturally toward the examination of fluid flow, since boundary element methods for fluids are at a much less developed state. This required the development of integral formulations for both the solid and fluid, and some preliminary infrastructural enhancements to a boundary element code to permit coupling of the fluid-structure problem. Boundary element formulations are implemented in two dimensions for both the solid and the fluid. The solid is modeled as an uncoupled thermoelastic medium under plane strain conditions, while several formulations are investigated for the fluid. For example, both vorticity and primitive variable approaches are implemented for viscous, incompressible flow, and a compressible version is developed. All of the above boundary element implementations are incorporated in a general purpose two-dimensional code. Thus, problems involving intricate geometry, multiple generic modeling regions, and arbitrary boundary conditions are all supported.

Dargush, Gary F.↗

Numerical procedures for three-dimensional computational surface thermochemistry

Models and equations for surface thermochemistry and near-surface thermophysics of aerodynamically-heated thermal protection materials are reviewed, with particular emphasis on computational boundary conditions for surface mass and energy transfer. The surface energy and mass balances, coupled with an appropriate ablation or surface catalysis model, provide complete thermochemical boundary conditions for a true multidisciplinary solution of the fully coupled fluid-dynamics/solid mechanics problem. Practical approximate solutions can be obtained by using a detailed model with full thermophysics for either the solid or fluid phase amd a semianalytic method for the other half of the problem. A significant increase in the state-of-the-art in aerothermal computational fluid dynamics is possible by uniting CFD methodology with surface thermochemistry boundary conditions and the heat-balance-integral method.

Milos, Frank S.↗

Flex-gear electrical power transmission

This study was conducted to develop an alternative way of transferring electricity across a continuously rotating joint, with little wear and the potential for low electrical noise. The problems with wires, slip rings, electromagnetic couplings, and recently invented roll-rings are discussed. Flex-gears, an improvement of roll-rings, are described. An entire class of flexgear devices is developed. Finally, the preferred flex-gear device is optimized for maximum electrical contact and analyzed for average mechanical power loss and maximum stress. For a device diameter of six inches, the preferred device is predicted to have a total electrical contact area of 0.066 square inches. In the preferred device, a small amount of internal sliding produces a 0.003 inch-pound torque that resists the motion of the device.

Vranish, John↗

Radiatively driven interactions between stratocumulus and synoptic waves

Quasi-geostrophic disturbances on a midlatitude beta-plane channel forced by radiative heating perturbations due to synoptic-scale variations of low-level stratiform cloud are considered. The longitudinal phase of the cloud is linked to that of the low-level streamfunction. Cloudiness is wavelike with cooling above cloudy areas that decreases exponentially with height. No perturbation cooling occurs above cloud-free areas. Two background flows are considered: one is constant and the other has a jet centered near the tropopause. Forced linear steady waves are found. Though infinitesimal amplitude disturbances are considered, the problem is nonlinear because of the coupling between cloud and winds. The resulting structures are sensitive to the phase shift between cloud and streamfunction with strongest coupling for cloud to the west of surface troughs. The waves have vertical scales on the order of the troposphere depth. The stationary structures capture the summertime pattern of stratocumulus off California and its linkage to the mid-Pacific ridge. When zonal-mean cloud cooling is allowed, the mean westerlies are strengthened in the northern half of the channel near the lower boundary. Synoptic-scale amplitudes respond to the mean cloud cooling by increasing (decreasing) just above the cloud in the northern (southern) part of the channel.

Clark, John H. E.↗