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 739 records · Page 41

Integration of Online Cross-Section Generation Capability with Depletion and Transient Solvers in Griffin

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE)-based reactor multiphysics analysis application jointly developed by Argonne and Idaho National Laboratories under the DOENE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. In FY25, an online crosssection generation capability based on the Self-Shielding Application Programming Interface (SSAPI) was demonstrated for TRISO-fueled reactor problems under steady-state conditions. This fiscal year, that capability was extended to support depletion and transient multiphysics calculations, enabling high-fidelity analyses that generate self-shielded cross sections on the fly from the actual evolving composition and temperature states rather than from pre-tabulated libraries. For depletion, a two-way coupling was established in which SSAPI computes compact-averaged self-shielded cross sections that the depletion solver then uses to advance the Bateman equations, with the updated compositions returned to SSAPI at each step; the depletion module was refactored to support both library-based and SSAPI-based cross sections, and additional logic was added to track daughter isotopes and to exclude minor isotopes for efficiency. For transient analysis, the SSAPI multigroup library was extended with the kinetics data required for time-dependent calculations, the Improved Quasi-Static (IQS) scheme was coupled with SSAPI, and several supporting capabilities were implemented, including a self-shielding treatment that lets control rods and drums move within a self-shielded model, which had previously been impossible and had ruled out rod- and drum-movement transients with on-the-fly cross sections altogether, a new mixing scheme for delayed-neutron precursor decay constants, a checkpoint-based restart workflow, and performance improvements such as pointwise cross-section interpolation and the bypassing of unnecessary Dancoff factor calculations. The implemented capabilities were verified against Serpent Monte Carlo solutions. For depletion, a prismatic pin-cell problem based on a Next Generation Nuclear Plant (NGNP) Very High Temperature Reactor benchmark showed excellent agreement, with eigenvalue differences within 200 pcm over the entire burnup range (up to 140 MWD/kgU) and fission-product and actinide inventories agreeing to within 0.8% and 2.5%, respectively; a heat-pipe microreactor assembly problem with a much higher fuel loading confirmed the same behavior and quantified the bias introduced when the multigroup equivalence effect is neglected. For transient analysis, a pin-cell problem with a step reactivity insertion and temperature feedback reproduced the analytically expected asymptotic power and showed close agreement between the direct and IQS solutions, and a two-dimensional microreactor core problem with control-drum rotation exercised the new moving-drum self-shielding treatment and demonstrated successful coupling of the online crosssection generation with both the direct and IQS transient methods. The capability was further exercised on a full-core pebble-bed problem, in which Griffin was coupled with the System Analysis Module (SAM) to simulate load-following operation of the gPBR with the Doppler feedback resolved at the TRISO fuel kernel temperature. These developments in Griffin provide a convenient, high-fidelity approach to cross-section generation for advanced thermal reactors with geometrically complex and highly heterogeneous configurations, including TRISO-fueled prismatic and pebble-bed systems, and support steady-state, depletion, and transient multiphysics calculations. They also enable self-shielded cross sections to be evaluated directly at the actual coupled state of the system, thereby establishing a foundation for high-fidelity, fully coupled multiphysics analysis of advanced reactors

Park, H.↗

LCTV Holographic Imaging

Astronauts are required to interface with complex systems that require sophisticated displays to communicate effectively. Lightweight, head-mounted real-time displays that present holographic images for comfortable viewing may be the ideal solution. We describe an implementation of a liquid crystal television (LCTV) as a spatial light modulator (SLM) for the display of holograms. The implementation required the solution of a complex set of problems. These include field calculations, determination of the LCTV-SLM complex transmittance characteristics and a precise knowledge of the signal mapping between the LCTV and frame grabbing board that controls it. Realizing the hologram is further complicated by the coupling that occurs between the phase and amplitude in the LCTV transmittance. A single drive signal (a gray level signal from a framegrabber) determines both amplitude and phase. Since they are not independently controllable (as is true in the ideal SLM) one must deal with the problem of optimizing (in some sense) the hologram based on this constraint. Solutions for the above problems have been found. An algorithm has been for field calculations that uses an efficient outer product formulation. Juday's MEDOF 7 (Minimum Euclidean Distance Optimal Filter) algorithm used for originally for filter calculations has been successfully adapted to handle metrics appropriate for holography. This has solved the problem of optimizing the hologram to the constraints imposed by coupling. Two laboratory methods have been developed for determining an accurate mapping of framegrabber pixels to LCTV pixels. A friendly software system has been developed that integrates the hologram calculation and realization process using a simple set of instructions. The computer code and all the laboratory measurement techniques determining SLM parameters have been proven with the production of a high quality test image.

Knopp, Jerome↗

Effect of Plasma Sheaths on Earth Re-entry MHD Processes

The first study of the full coupling between the aerothermodynamics, the magnetohydrodynamics, and the plasma sheaths within earth entry flows is here per-formed. The problem tackled herein is representative of a force-generating MHD patch located between the plenum and the frustum of a capsule entering the earth’s atmosphere at Mach 34. The reactions are obtained from a modified 11-species Park chemical solver. The physical model fully couples the drift-diffusion model for the sheaths to the multispecies Navier-Stokes equations for the plasma flow. The Hall and ion slip effects are taken into consideration within the plasma flow and within the sheaths. Results obtained indicate that the cathode sheath can exert significant resistance to current flow and thus reduce the MHD-generated force by 40% or more. The effect of the sheath on the MHD force is found to be more pronounced at lower applied magnetic field.

Earth Re-entry↗

Effect of Plasma Sheaths on Earth Re-entry MHD Processes

The first study of the full coupling between the aerothermodynamics, the magnetohydrodynamics, and the plasma sheaths within earth entry flows is here per-formed. The problem tackled herein is representative of a force-generating MHD patch located between the plenum and the frustum of a capsule entering the earth’s atmosphere at Mach 34. The reactions are obtained from a modified 11-species Park chemical solver. The physical model fully couples the drift-diffusion model for the sheaths to the multispecies Navier-Stokes equations for the plasma flow. The Hall and ion slip effects are taken into consideration within the plasma flow and within the sheaths. The effect of the electrode material on the MHD process is studied. Using thoriated tungsten instead of graphite leads to a four-fold increase in the Lorentz forces and also leads to significantly reduced heat fluxes on the cathode. This is attributed to the much higher electrical conductivity of the thoriated tungsten sheath reducing by orders of magnitude the plasma electrical resistance near the surfaces.

Bernard Parent↗

Predictive dynamic wetting, fluid–structure interaction simulations for braze run-out

Brazing and soldering are metallurgical joining techniques that use a wetting molten metal to create a joint between two faying surfaces. Here, the quality of the brazing process depends strongly on the wetting properties of the molten filler metal, namely the surface tension and contact angle, and the resulting joint can be susceptible to various defects, such as run-out and underfill, if the material properties or joining conditions are not suitable. In this work, we implement a finite element simulation to predict the formation of such defects in braze processes. This model incorporates both fluid–structure interaction through an arbitrary Eulerian–Lagrangian technique and free surface wetting through conformal decomposition finite element modeling. Upon validating our numerical simulations against experimental run-out studies on a silver-Kovar system, we then use the model to predict run-out and underfill in systems with variable surface tension, contact angles, and applied pressure. Finally, we consider variable joint/surface geometries and show how different geometrical configurations can help to mitigate run-out. This work aims to understand how brazing defects arise and validate a coupled wetting and fluid–structure interaction simulation that can be used for other industrial problems.

36 MATERIALS SCIENCE↗

Neutron stars with exceptionally light QCD axions

We present a comprehensive study of axion condensed neutron stars that arise in models of an exceptionally light axion that couples to quantum chromodynamics (QCD). These axions solve the strong-charge-parity (𝐶⁢𝑃) problem, but have a mass-squared lighter than that due to QCD by a factor of 𝜖 < 1. Inside dense matter, the axion potential is altered, and much of the matter in neutron stars resides in the axion condensed phase where the strong-𝐶⁢𝑃 parameter 𝜃 = 𝜋 and 𝐶⁢𝑃 remains a good symmetry. In these regions, masses and interactions of nucleons are modified, in turn changing the equation of state (EOS), structure, and phenomenology of the neutron stars. We take the first steps toward the study of the EOS of neutron star matter at 𝜃 = 𝜋 within chiral effective field theory and use relativistic mean field theory to deduce the resulting changes to nuclear matter and the neutron star low-density EOS. We derive constraints on the exceptionally light axion parameter space based on observations of the thermal relaxation of accreting neutron stars, isolated neutron star cooling, and pulsar glitches, excluding the region up to 5 ×10 −7 ≲ 𝜖 ≲ 0.2 for 𝑚 𝑎 ≳ 2 × 10 −9 eV. We comment on potential changes to the neutron star mass-radius relationship, and discuss the possibility of novel, nuclear-density compact objects with 𝜃 = 𝜋 that are stabilized not by gravity but by the axion potential.

axions↗

Multi-Period Optimization of Multi-Timescale Energy Systems: Application to Solid-Oxide Electrolysis Cells

A presentation detailing the recent progress in the optimization of Solid-Oxide Cells under Chemical Degradation. The existing methods typically involve a quasi-steady state assumption to make this large-scale problem tractable. Here, we introduce an extension to this method that allows for adaptive coupling of degradation with the fast-timescale process based on a variety of error thresholds. This reduces the error accumulated during long-term optimization of the SOC under degradation and allows for decision making at multiple timescales. The results presented include long-term operating profiles of the SOC under degradation with both steady-state and fluctuating day to day operation.

Giridhar, Nishant↗

The ORNL Moderator Test Station Science Case

Oak Ridge National Laboratory (ORNL) hosts two world-leading slow neutron sources, the Spallation Neutron Source (SNS) and the High Flux Isotope Reactor (HFIR), and is currently developing the technical design for a Second Target Station (STS) for the SNS. Upon completion of the STS project, ORNL will be uniquely positioned to optimize each of its three neutron sources, the SNS First Target Station (FTS), the STS, and HFIR, in a complementary way. Among the essential aspects of a re-imagined FTS and the current STS design are high-brightness parahydrogen moderators—moderators which are optimized for high per-unit-area neutron brightness rather than integrated-across-large-area neutron intensity. The high-brightness moderators proposed for the STS will, for the brightness metric, significantly outperform the coupled moderators currently on the FTS for appropriately optimized neutron beamlines, provided the moderating hydrogen is converted to near-equilibrium levels of parahydrogen (approximately 99.8% at 20 K). The original FTS moderators, by contrast, were conservatively designed to be relatively insensitive to the exact ortho:para ratio, with a consequent loss in performance. As a result, a redesign of the FTS moderators assuming fully converted parahydrogen could result in significant performance improvements on the FTS coupled moderators, and more consistent performance over time for all hydrogen moderators. This “parahydrogen problem” is a long-standing challenge for the effective implementation of hydrogen cold moderators at high-power neutron sources. In addition, the development of new moderator concepts, whether based on previously unused materials, structured heterogeneous arrays, or even simply on changes in overall shape and size is significantly restricted at a large-scale production facility intended to use the resulting neutron beams. Accordingly, moderators for production neutron sources are often designed in a very conservative, low-risk fashion, even though this compromises the absolute neutronic performance. Advanced moderator concepts worthy of study include features that could not be tested without redesigning and redeploying the entire existing reflector, shielding, and neutron beamline installation, making such development efforts far more expensive than building a stand-alone test facility.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Charged particle motion in the vicinity of a neutral plane

The motion of a charged particle in a 2-dimensional neutral sheet with linear magnetic field variation is analyzed by Hamiltonian methods. From the basic Hamiltonian formulation, results of Sonnerup are recovered and given intuitive interpretation. The transformed Hamiltonian is used to derive the correct frequency of oscillation and serves as the basis of analytical treatment of perturbed versions of the motion, e.g. with electric field added. For 2-dimensional fields with slightly differing configurations - in particular, with added small field component orthogonal to the sheet, with an X-type null or with a chain of null points - two alternative methods are developed, reducing the problem either to motion in a 2-dimensional potential or to a pair of coupled oscillators.

Stern, D.↗

Low-cost calibration of acoustic locators

Method uses modified commerically-available piezoelectric-torch lighter. Handheld lighter has controlled spark gap that can be easily adjusted to produce repeatable short-duration high-amplitude voltage spikes. Pulser and lighter are coupled via short axial cable, eliminating long cable run variations in cable attenuation, and problem with cable entangling with anything in its path.

Berry, R. F.↗

Development of basic theories and techniques for determining stresses in rotating turbine or compressor blades

A method for measuring in-plane displacement of a rotating structure by using two laser speckle photographs is described. From the displacement measurements one can calculate strains and stresses due to a centrifugal load. This technique involves making separate speckle photographs of a test model. One photograph is made with the model loaded (model is rotating); the second photograph is made with no load on the model (model is stationary). A sandwich is constructed from the two speckle photographs and data are recovered in a manner similar to that used with conventional speckle photography. The basic theory, experimental procedures of this method, and data analysis of a simple rotating specimen are described. In addition the measurement of in-plane surface displacement components of a deformed solid, and the application of the coupled laser speckle interferometry and boundary-integral solution technique to two dimensional elasticity problems are addressed.

Chien, C. H.↗

Identification of helicopter rotor dynamic models

A recursive, extended Kalman-filter approach is applied to the identifiction of rotor damping levels of representative helicopter dynamic systems. The general formulation of the approach is presented in the context of a typically posed stochastic estimation problem, and the method is analytically applied to determining the damping levels of a coupled rotor-body system. The identified damping covergence characteristics are studied for sensitivity to both constant-coefficient and periodic-coefficient measurement models, process-noise covariance levels, and specified initial estimates of the rotor-system damping. A second application of the method to identifying the plant model for a highly damped, isolated flapping blade with a constant-coefficient state model (hover) and a periodic-coefficient state model (forward flight) is also investigated. The parameter-identification capability is evaluated for the effect of periodicity on the plant model coefficients and the influence of different measurement noise levels.

Molusis, J. A.↗

Analysis of rectangular microstrip antennas

The problem of microstrip antennas covered by a dielectric substrate is formulated in terms of coupled integro-differential equations with the current distribution on the conducting patch as an unknown quantity. The Galerkin method is used to solve for the unknown patch current. Using the present formulation, the radiation pattern, the resonant frequency, and the bandwidth of a rectangular microstrip antenna are computed. Design data for a rectangular microstrip antenna are also presented.

Bailey, M. C.↗

A high reliability battery management system

Over a period of some 5 years Canadian Astronautics Limited (CAL) has developed a system to autonomously manage, and thus prolong the life of, secondary storage batteries. During the development, the system was aimed at the space vehicle application using nickel cadmium batteries, but is expected to be able to enhance the life and performance of any rechargeable electrochemical couple. The system handles the cells of a battery individually and thus avoids the problems of over, and under, drive that inevitably occur in a battery of cells managed by an averaging system. This individual handling also allow cells to be totally bypassed in the event of failure, thus avoiding the losses associated with low capacity, partial short circuit, and the catastrophe of open circuit. The system has an optional capability of managing redundant batteries simultaneously, adding the advantage of on line reconditioning of one battery, while the other maintains the energy storage capability of the overall system. As developed, the system contains a dedicated, redundant, microprocessor, but the capability exists to have this computing capability time shared, or remote, and operating through a data link. As adjuncts to the basic management system CAL has developed high efficiency, polyphase, power regulators for charge and discharge power conditioning.

Moody, M. H.↗

Recent trends in rotary-wing aeroelasticity

The purpose of this paper is to survey the principal developments which have occurred in the field of rotary-wing aeroelasticity during the past five-year period. This period has been one of considerable activity and approximately one hundred papers have been published on this topic. To facilitate this review the field has been divided into a number of areas in which concentrated research activity has taken place. The main areas in which recent research is reviewed are: (1) structural modeling; (2) aerodynamic modeling; (3) aeroelastic problem formulation using automated or computerized methods; (4) aeroelastic analyses in forward flight; (5) coupled rotor/fuselage analyses; (6) active controls and their application to aeroelastic response and stability; (7) application of structural optimization to vibration response; and (8) aeroelastic analysis and testing of special configurations. These areas are reviewed with different levels of detail and some useful observations regarding potentially rewarding areas of future research are made.

Friedmann, Peretz P.↗

Recent trends in rotary-wing aeroelasticity

The purpose of this paper is to survey the principal developments which have occurred in the field of rotary-wing aeroelasticity during the past five year period. This period has been one of considerable activity and approximately one hundred papers have been published on this topic. To facilitate this review the field has been divided into a number of areas in which concentrated research activity has taken place. The main areas in which recent research is reviewed are: (1) structural modeling; (2) aerodynamic modeling; (3) aeroelastic problem formulation using automated or computerized methods; (4) aeroelastic analyses in forward flight; (5) coupled rotor/fuselage analyses; (6) active controls and their application to aeroelastic response and stability; (7) application of structural optimization to vibration reduction; and (8) aeroelastic analysis and testing of special configurations. These areas are reviewed with different levels of detail and some useful observations regarding potentially rewarding areas of future research are made.

Friedmann, Peretz P.↗

A surface crack in shells under mixed-mode loading conditions

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

Joseph, P. F.↗

Heat pipe dynamic behavior

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

Issacci, F.↗