Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear stability analysis”

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,513 records · Page 84

Optimizing district energy systems by integrating Borehole Thermal Energy Storage Using a Mixed-Integer Linear Programming g-function framework with a Multi-Timescale Rolling Horizon method

Shallow geothermal has gained increasing attention in recent years; however, a reliable framework for its accurate incorporation into large-scale energy system optimization remains lacking. This study proposes a Mixed-Integer Linear Programming (MILP) framework combined with the g-function approach to integrate Borehole Thermal Energy Storage (BTES) technology into energy system optimization. Validation against a Modelica-based reservoir network simulation demonstrates that the proposed framework effectively captures the ground thermal response under varying energy loads and accurately estimates the borefield energy supply. To enhance scalability, a Rolling Horizon with Multi-Timescale (RH-MTS) method is further introduced, reducing computational time by 73 % for the 1-year optimization model with only minor loss of optimality. The framework is demonstrated through the case study of the UC Berkeley campus. Results indicate that BTES is a cost-effective and low-carbon solution: two borefields comprising 382 boreholes can meet 8.0 % and 6.6 % of the total campus heating and cooling demand, respectively, at an average energy rate of 0.70–0.77 USD/kWh and carbon intensity of 0.54 kg-CO2/kWh. Short-term analysis reveals a 35%–65% decline in BTES energy flow after 3–6 months of continuous heating/cooling operation, while long-term simulation shows that annual energy production of BTES can vary by up to 12.0 % after four years before stabilizing. Overall, this study develops a novel optimization framework that couples physics-based g-function method with MILP optimization framework, thereby advancing methodological development for shallow-geothermal integration and providing actionable guidance for BTES deployment in district-energy systems.

Yang, Jiahui↗

Clemson University Wind Turbine Drivetrain Test Facility

In November of 2009, Clemson University was awarded a competitive grant from the U.S. Department of Energy to design, build and operate a facility for full-scale, highly accelerated mechanical testing of next-generation wind turbine drivetrain technologies. The primary goal of the project was to design, construct, commission, and operate a state-of-the-art sustainable facility that permits full-scale highly accelerated testing of advanced drivetrain systems for large wind turbines. The secondary goal was to meet the objectives of the American Recovery and Reinvestment Act of 2009, especially in job creation, and provide a positive impact on economically distressed areas in the United States, and preservation and economic recovery in an expeditious manner. The project was executed according to a managed cooperative agreement with the Department of Energy and was an extraordinary success. The resultant new facility is located in North Charleston, SC, providing easy transportation access by rail, road or ship and operates on an open access model such that it is available to the U.S. Wind Industry for research, analysis, and evaluation activities. The 72 m by 97 m facility features two mechanical dynamometer test bays for evaluating the torque and blade dynamic forces experienced by the rotors of wind turbine drivetrains. The dynamometers are rated at 7.5 MW and 15 MW of low speed shaft power and are configured as independent test areas capable of simultaneous operation. All six degrees of freedom, three linear and three rotational, for blade and rotor dynamics are replicated through the combination of a drive motor, speed reduction gearbox and a controllable hydraulic load application unit (LAU). This new LAU setup readily supports accelerated lifetime mechanical testing and load analysis for the entire drivetrain system of the nacelle and easily simulates a wide variety of realistic operating scenarios in a controlled laboratory environment. The development of these two dynamometer test rigs is the first significant achievement for the project. These test rigs embody a new manner of test due to the system configuration and completely new design with a free floating loading hub in the LAU. This project provided the catalyst for the advancement to this new test rig configuration that has been adopted by every significant wind turbine test rig constructed since the inception of this project. There are currently two different vendors supplying these new systems. Catalyzing this new design is the second major success of the project. With the increased market penetration of wind energy over the past decade, many regions and countries have developed specific electrical grid specifications and performance codes for large wind farms to ensure operational reliability and stability. These grid codes provide requirements for interconnection with the grid during low or high voltage phenomena, typically encountered during and after system fault events. Given the installed infrastructure of the Wind Turbine Drivetrain Testing Facility (WTDTF), a natural expansion of facility capability was to include the necessary equipment for performing fault ride-through evaluations of wind turbines to the Low Voltage Ride Through (LVRT) codes. Once the decision was made to expand the scope of the original grant into fault ride-through testing, it was clear that there are several markets, not just wind, which could benefit from this type of test and that simple fault ride-through testing could be extended into a broader scope of electrical testing capabilities. It was at this point that Clemson University was awarded a second grant to build a 15 MW Hardware-In-the-Loop (HIL) Grid Simulator in order to establish world class electrical testing capabilities to compliment the mechanical testing at the WTDTF. This third significant achievement resulted in the 15 MW HIL Grid Simulator as the corner stone of the Duke Energy eGRID Center and is collocated with the WTDTF in the SCE&G Energy Innovation Center at the new Facility in North Charleston, SC, USA. Before the eGRID was completed, it was recognized that the ability to test solar farm equipment was but a small step away thru the addition of enhanced equipment to provide for DC testing. In yet another expansion/success, a 2.5 MW rectifier system was designed and implemented by Clemson staff to enhance the Center’s capabilities. The program required over 250,000 man-hours of on-site construction labor reworking the brownfield facility on the former Navy Base, clearly satisfying one of the major goals of the Reinvestment Act. This was done while winning numerous awards for design and construction of the facility, including the Top US Project for 2014 from the Trade Journal Engineering News Record. The project was a major collaborative developmental activity managed by Clemson University staff that involved the DOE and many partners and organizations.

17 WIND ENERGY↗

Organometallic Polymeric Conductors

For aerospace applications, the use of polymers can result in tremendous weight savings over metals. Suitable polymeric materials for some applications like EMI shielding, spacecraft grounding, and charge dissipation must combine high electrical conductivity with long-term environmental stability, good processability, and good mechanical properties. Recently, other investigators have reported hybrid films made from an electrically conductive polymer combined with insulating polymers. In all of these instances, the films were prepared by infiltrating an insulating polymer with a precursor for a conductive polymer (either polypyrrole or polythiophene), and oxidatively polymerizing the precursor in situ. The resulting composite films have good electrical conductivity, while overcoming the brittleness inherent in most conductive polymers. Many aerospace applications require a combination of properties. Thus, hybrid films made from polyimides or other engineering resins are of primary interest, but only if conductivities on the same order as those obtained with a polystyrene base could be obtained. Hence, a series of experiments was performed to optimize the conductivity of polyimide-based composite films. The polyimide base chosen for this study was Kapton. 3-MethylThiophene (3MT) was used for the conductive phase. Three processing variables were identified for producing these composite films, namely time, temperature, and oxidant concentration for the in situ oxidation. Statistically designed experiments were used to examine the effects of these variables and synergistic/interactive effects among variables on the electrical conductivity and mechanical strength of the films. Multiple linear regression analysis of the tensile data revealed that temperature and time have the greatest effect on maximum stress. The response surface of maximum stress vs. temperature and time (for oxidant concentration at 1.2 M) is shown. Conductivity of the composite films was measured for over 150 days in air at ambient temperature. The conductivity of the films dropped only half an order of magnitude in that time. Films aged under vacuum at ambient temperature diminished slightly in conductivity in the first day, but did not change thereafter. An experimental design approach will be applied to maximize the efficiency of the laboratory effort. The material properties (initial and long term) will also be monitored and assessed. The experimental results will add to the existing database for electrically conductive polymer materials. Attachments: 1) Synthesis Crystal Structure, and Polymerization of 1,2:5,6:9,10-Tribenzo-3,7,11,13-tetradehydro(14) annulene. 2) Reinvestigation of the Photocyclization of 1,4-Phenylene Bis(phenylmaleic anhydride): Preparation and Structure of (5)Helicene 5,6:9,10-Dianhydride. 3) Preparation and Structure Charecterization of a Platinum Catecholate Complex Containing Two 3-Ethynyltheophone Groups. and 4) Rigid-Rod Polymers Based on Noncoplanar 4,4'-Biphenyldiamines: A Review of Polymer Properties vs Configuration of Diamines.

Youngs, Wiley J.↗

AOCS Performance and Stability Validation for a 160-m Solar Sail with Control-Structure Interactions

Future solar sail missions, such as NASA's Solar Polar Imager Vision, will require sails with dimensions on the order of 50-500 m. We are examining a square sail design with moving mass (trim control mass, TCM) and quadrant rotation primary actuators plus pulsed plasma thrusters (PPTs) at the mast tips for backup attitude control. Quadrant rotation is achieved via roll stabilizer bars (RSB) at the mast tips. At these sizes, given the gossamer nature of the sail supporting structures, flexible modes may be low enough to interact with the control system, especially as these actuators are located on the flexible structure itself and not on the rigid core. This paper develops a practical analysis of the flexible interactions using state-space systems and modal data from finite element models of the system. Torsion and bending of the masts during maneuvers could significantly affect the function of the actuators while activation of the membrane modes could adversely affect the thrust vector direction and magnitude. Analysis of the RSB and TCM dynamics for developing high-fidelity simulations is included. For control analysis of the flexible system, standard finite-element models of the flexible sail body are loaded and the modal data is used to create a modal coordinate state-space system. Key parameters include which modes to include, which nodes are of interest for force inputs and displacement outputs, connecting nodes through which external forces and torques are applied from the flex body to the core, any nominal momentum in the system, and any steady rates. The system is linearized about the nominal attitude and rate. The state-space plant can then be analyzed with a state-space controller, and Bode, Nyquist, step and impulse responses generated. The approach is general for any rigid core with a flexible appendage. This paper develops a compensator for a simple two-mass flex system and extrapolates the results to the solar sail. A finite element model of the 20 m solar sail by ATK Space Systems, recently validated in ground tests, is used to demonstrate the sail analysis approach.

Wie, Bong↗

Active stabilization to prevent surge in centrifugal compression systems

This report documents an experimental and analytical study of the active stabilization of surge in a centrifugal engine. The aims of the research were to extend the operating range of a compressor as far as possible and to establish the theoretical framework for the active stabilization of surge from both an aerodynamic stability and a control theoretic perspective. In particular, much attention was paid to understanding the physical limitations of active stabilization and how they are influenced by control system design parameters. Previously developed linear models of actively stabilized compressors were extended to include such nonlinear phenomena as bounded actuation, bandwidth limits, and robustness criteria. This model was then used to systematically quantify the influence of sensor-actuator selection on system performance. Five different actuation schemes were considered along with four different sensors. Sensor-actuator choice was shown to have a profound effect on the performance of the stabilized compressor. The optimum choice was not unique, but rather shown to be a strong function of some of the non-dimensional parameters which characterize the compression system dynamics. Specifically, the utility of the concepts were shown to depend on the system compliance to inertia ratio ('B' parameter) and the local slope of the compressor speedline. In general, the most effective arrangements are ones in which the actuator is most closely coupled to the compressor, such as a close-coupled bleed valve inlet jet, rather than elsewhere in the flow train, such as a fuel flow modulator. The analytical model was used to explore the influence of control system bandwidth on control effectiveness. The relevant reference frequency was shown to be the compression system's Helmholtz frequency rather than the surge frequency. The analysis shows that control bandwidths of three to ten times the Helmholtz frequency are required for larger increases in the compressor flow range. This has important implications for implementation in gas turbine engines since the Helmholtz frequencies can be over 100 Hz, making actuator design extremely challenging.

Epstein, Alan H.↗

A positive and entropy-satisfying finite volume scheme for the Baer–Nunziato model

We present a relaxation scheme for approximating the entropy dissipating weak solutions of the Baer–Nunziato two-phase flow model. This relaxation scheme is straightforwardly obtained as an extension of the relaxation scheme designed in for the isentropic Baer–Nunziato model and consequently inherits its main properties. To our knowledge, this is the only existing scheme for which the approximated phase fractions, phase densities and phase internal energies are proven to remain positive without any restrictive condition other than a classical fully computable CFL condition. For ideal gas and stiffened gas equations of state, real values of the phasic speeds of sound are also proven to be maintained by the numerical scheme. It is also the only scheme for which a discrete entropy inequality is proven, under a CFL condition derived from the natural sub-characteristic condition associated with the relaxation approximation. This last property, which ensures the non-linear stability of the numerical method, is satisfied for any admissible equation of state. We provide a numerical study for the convergence of the approximate solutions towards some exact Riemann solutions. The numerical simulations show that the relaxation scheme compares well with two of the most popular existing schemes available for the Baer–Nunziato model, namely Schwendeman–Wahle–Kapila's Godunov-type scheme and Tokareva–Toro's HLLC scheme . The relaxation scheme also shows a higher precision and a lower computational cost (for comparable accuracy) than a standard numerical scheme used in the nuclear industry, namely Rusanov's scheme. Finally, we assess the good behavior of the scheme when approximating vanishing phase solutions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Amphiphilic Bottlebrush Block Copolymers: Analysis of Aqueous Self-Assembly by Small-Angle Neutron Scattering and Surface Tension Measurements

A systematic series of 16 amphiphilic bottlebrush block copolymers (BCPs) containing polystyrene and poly(N-acryloylmorpholine) (PACMO) side chains were prepared by a combination of atom-transfer radical polymerization (ATRP), photoiniferter polymerization, and ring-opening metathesis polymerization (ROMP). The grafting-through method used to prepare the polymers enabled a high degree of control over backbone and side-chain molar masses for each block. Surface tension measurements on the self-assembled amphiphilic bottlebrush BCPs in water revealed an ultralow critical micelle concentration (cmc), 1–2 orders of magnitude lower than linear BCP analogues on a molar basis, even for micelles with >90% PACMO content. Combined with coarse-grained molecular dynamics simulations, fitting of small-angle neutron scattering traces (SANS) allowed us to evaluate solution conformations for individual bottlebrush BCPs and micellar nanostructures for self-assembled macromolecules. Bottlebrush BCPs showed an increase in anisotropy with increasing PACMO content in toluene-d 8 , which is a good solvent for both blocks, reflecting an extended conformation for the PACMO block. SANS traces of bottlebrush BCPs assembled into micelles in D 2 O, a selective solvent for PACMO, were fitted to a core–shell–shell model, suggesting the presence of a partially hydrated inner shell. Here, results showed an average micelle diameter of 40 nm with combined shell diameters ranging from 16 to 18 nm. A general trend of increased stability of micelles (i.e., resistance to precipitation) was observed with increases in PACMO content. These results demonstrate the stability of bottlebrush polymer micelles, which self-assemble to form spherical micelles with ultralow (<70 nmol/L) cmc’s across a broad range of compositions.

36 MATERIALS SCIENCE↗

Solitary Waves of a $\mathcal {P}$ $\mathcal {T}$ -Symmetric Nonlinear Dirac Equation

In our study we consider we consider a prototypical example of a mathcalP mathcalT-symmetric Dirac model. We discuss the underlying linear limit of the model and identify the threshold of the mathcalP mathcalT -phase transition in an analytical form. We then focus on the examination of the nonlinear model. We consider the continuation in the mathcalP mathcalT -symmetric model of the solutions of the corresponding Hamiltonian model and find that the solutions can be continued robustly as stable ones all the way up to the mathcalP mathcalT-transition threshold. In the latter, they degenerate into linear waves. We also examine the dynamics of the model. Given the stability of the waveforms in the mathcalP mathcalT-exact phase, we consider them as initial conditions for parameters outside of that phase. We also find that both oscillatory dynamics and exponential growth may arise, depending on the size of the corresponding “quench”. The former can be characterized by an interesting form of bifrequency solutions that have been predicted on the basis of the SU symmetry. Finally, we explore some special, analytically tractable, but not mathcalP mathcalT-symmetric solutions in the massless limit of t- e model.

97 MATHEMATICS AND COMPUTING↗

Numerical Study of Steady and Unsteady Canard-Wing-Body Aerodynamics

The use of canards in advanced aircraft for control and improved aerodynamic performance is a topic of continued interest and research. In addition to providing maneuver control and trim, the influence of canards on wing aerodynamics can often result in increased maximum lift and decreased trim drag. In many canard-configured aircraft, the main benefits of canards are realized during maneuver or other dynamic conditions. Therefore, the detailed study and understanding of canards requires the accurate prediction of the non-linear unsteady aerodynamics of such configurations. For close-coupled canards, the unsteady aerodynamic performance associated with the canard-wing interaction is of particular interest. The presence of a canard in close proximity to the wing results in a highly coupled canard-wing aerodynamic flowfield which can include downwash/upwash effects, vortex-vortex interactions and vortex-surface interactions. For unsteady conditions, these complexities of the canard-wing flowfield are further increased. The development and integration of advanced computational technologies provide for the time-accurate Navier-Stokes simulations of the steady and unsteady canard-wing-body flox,fields. Simulation, are performed for non-linear flight regimes at transonic Mach numbers and for a wide range of angles of attack. For the static configurations, the effects of canard positioning and fixed deflection angles on aerodynamic performance and canard-wing vortex interaction are considered. For non-static configurations, the analyses of the canard-wing body flowfield includes the unsteady aerodynamics associated with pitch-up ramp and pitch oscillatory motions of the entire geometry. The unsteady flowfield associated with moving canards which are typically used as primary control surfaces are considered as well. The steady and unsteady effects of the canard on surface pressure integrated forces and moments, and canard-wing vortex interaction are presented in detail including the effects of the canard on the static and dynamic stability characteristics. The current study provides an understanding of the steady and unsteady canard-wing-body flowfield. Emphasis is placed on the effects of the canard on aerodynamic performance as well as the detailed flow physics of the canard-wing flowfield interactions. The computational tools developed to accurately predict the time-accurate flowfield of moving canards provides for the capability of coupled fluids-controls simulations desired in the detailed design and analysis of advanced aircraft.

Eugene, L. Tu↗

Modeling the Effects of a Backward-Facing Step on Boundary-Layer Transition

We model transition to turbulence in a two-dimensional boundary layer downstream of a backward-facing step (BFS) along a flat plate. With the goal of evaluating the available engineering models for predicting the effects of step excrescences on the transition characteristics, two separate methodologies are used to monitor the streamwise shift in the transition onset location as the step height and the flow speed are varied across the range of a previously reported experiment involving step-height-to-local-displacement-thickness ratios of 0 < h/δ* < 1.6. Unlike the variable N -factor method from the previous literature, both of these methods are general in scope and do not involve any empirical correlations that are specific to step excrescences. The first of these techniques involves an N -factor method that directly accounts for the change in boundary-layer instability characteristics due to the step. Stability computations using the harmonic linearized Navier-Stokes equations (HLNSE), which fully account for the nonparallel-mean-flow effects close to the BFS, indicate that the measured transition locations at nearly all test conditions ( h/δ* < 1.3) correlate well with a computed N -factor of N tr = 7.6, demonstrating a successful stability-based transition criterion related to step excrescences. Linear stability theory, which does not account for nonparallel effects, demonstrates reasonable agreement with the HLNSE results, yielding good predictions for the overall trends, but predicts a somewhat earlier onset of transition than HLNSE. The other methodology used in this work involves transport-equation-based transition models. We first show that the Langtry-Menter y - Re θt transition model cannot accurately predict the location of transition onset for moderate BFS heights because it is unable to accurately account for the flow history effects. Along with the Langtry-Menter transition model, we also show the amplification factor transport model does not produce accurate transition locations for subsonic flow over steps even though it accounts for some flow history effects.

Transition↗

Modeling the Effects of a Backward-Facing Step on Boundary-Layer Transition

We model transition to turbulence in a two-dimensional boundary layer downstream of a backward-facing step (BFS) along a flat plate. With the goal of evaluating the available engineering models for predicting the effects of step excrescences on the transition characteristics, two separate methodologies are used to monitor the streamwise shift in the transition onset location as the step height and the flow speed are varied across the range of a previously reported experiment involving step-height-to-local-displacement-thickness ratios of 0 < h/δ* < 1.6. Unlike the variable N -factor method from the previous literature, both of these methods are general in scope and do not involve any empirical correlations that are specific to step excrescences. The first of these techniques involves an N -factor method that directly accounts for the change in boundary-layer instability characteristics due to the step. Stability computations using the harmonic linearized Navier-Stokes equations (HLNSE), which fully account for the nonparallel-mean-flow effects close to the BFS, indicate that the measured transition locations at nearly all test conditions ( h/δ* < 1.3) correlate well with a computed N -factor of N tr = 7.6, demonstrating a successful stability-based transition criterion related to step excrescences. Linear stability theory, which does not account for nonparallel effects, demonstrates reasonable agreement with the HLNSE results, yielding good predictions for the overall trends, but predicts a somewhat earlier onset of transition than HLNSE. The other methodology used in this work involves transport-equation-based transition models. We first show that the Langtry-Menter y - Re θt transition model cannot accurately predict the location of transition onset for moderate BFS heights because it is unable to accurately account for the flow history effects. Along with the Langtry-Menter transition model, we also show the amplification factor transport model does not produce accurate transition locations for subsonic flow over steps even though it accounts for some flow history effects.

Transition↗

Feedback Control of Plasma Instabilities with Charged Particle Beams and Study of Plasma Turbulence

A new non-perturbing technique for feedback control of plasma instabilities has been developed in the Columbia Linear Machine (CLM). The feedback control scheme involves the injection of a feedback modulated ion beam as a remote suppressor. The ion beam was obtained from a compact ion beam source which was developed for this purpose. A Langmuir probe was used as the feedback sensor. The feedback controller consisted of a phase-shifter and amplifiers. This technique was demonstrated by stabilizing various plasma instabilities to the background noise level, like the trapped particle instability, the ExB instability and the ion-temperature-gradient (ITG) driven instability. An important feature of this scheme is that the injected ion beam is non-perturbing to the plasma equilibrium parameters. The robustness of this feedback stabilization scheme was also investigated. The principal result is that the scheme is fairly robust, tolerating about 100% variation about the nominal parameter values. Next, this scheme is extended to the unsolved general problem of controlling multimode plasma instabilities simultaneously with a single sensor-suppressor pair. A single sensor-suppressor pair of feedback probes is desirable to reduce the perturbation caused by the probes. Two plasma instabilities the ExB and the ITG modes, were simultaneously stabilized. A simple 'state' feedback type method was used where more state information was generated from the single sensor Langmuir probe by appropriate signal processing, in this case, by differentiation. This proof-of-principle experiment demonstrated for the first time that by designing a more sophisticated electronic feedback controller, many plasma instabilities may be simultaneously controlled. Simple theoretical models showed generally good agreement with the feedback experimental results. On a parallel research front, a better understanding of the saturated state of a plasma instability was sought partly with the help of feedback. A plasma instability is usually observed in its saturated state and appears as a single feature in the frequency spectrum with a single azimuthal and parallel wavenumbers. The physics of the non-zero spectral width was investigated in detail because the finite spectral width can cause "turbulent" transport. One aspect of the "turbulence" was investigated by obtaining the scaling of the linear growth rate of the instabilities with the fluctuation levels. The linear growth rates were measured with the established gated feedback technique. The research showed that the ExB instability evolves into a quasi-coherent state when the fluctuation level is high. The coherent aspects were studied with a bispectral analysis. Moreover, the single spectral feature was discovered to be actually composed of a few radial harmonics. The radial harmonics play a role in the nonlinear saturation of the instability via three-wave coupling.

Tham, Philip Kin-Wah↗

Determination of He and D permeability of neutron-irradiated SiC tubes to examine the potential for release due to micro-cracking

Driven by the need to enlarge the safety margins of light water reactors in both design-basis and beyond-design-basis accident scenarios, the research and development of accident-tolerant fuel (ATF) has become an importance topic in the nuclear engineering and materials community. Continuous SiC fiber-reinforced SiC matrix ceramic composites are under consideration as a replacement for traditional zirconium alloy cladding owing to their high-temperature stability, chemical inertness, and exceptional irradiation resistance. Among the key technical feasibility issues, potential failure of the fission product containment due to probabilistic penetrating cracking has been identified as one of the two most critical feasibility issues, together with the radiolysisassisted hydrothermal corrosion of SiC. The experimental capability to evaluate the hermeticity of SiC-based claddings is an urgent need. In this report, we present the development of a comprehensive permeation testing station established in the Low Activation Materials Development and Analysis laboratory at Oak Ridge National Laboratory. Preliminary results for the hermeticity evaluation of un-irradiated monolithic SiC tubes, uncoated and coated SiC/SiC composite tubes, and neutron-irradiated monolithic SiC tubes at room temperature are exhibited. The results indicate that this new permeation testing station is capable of evaluating the hermeticity of SiC-based tubes by determining the helium and deuterium permeation flux as a function of gas pressure at a high resolution of 8.07 x 10 -12 atm-cc/s for helium and 2.83 x 10 -12 atm-cc/s for deuterium, respectively. The detection limit of this system is sufficient to evaluate the maximum allowable helium leakage rate of lab-scale tubular samples, which is linearly extrapolated from the evaluation standard used for a commercial as-manufactured light water reactor fuel rod at room temperature. The un-irradiated monolithic SiC tube is hermetic, as is manifested by the un-detectable deuterium permeation flux at various feeding gas pressures. A large helium leakage rate was detected for the uncoated SiC/SiC composite tube exposed to atmosphere, indicating it is inherently not hermetic. The hermeticity of coated SiC/SiC composite tubes is strongly dependent on the coating materials and the preparation of the substrate SiC/SiC composite samples. To simulate the practical application environment, monolithic CVD SiC tubes were exposed to neutron irradiation at the High Flux Isotope Reactor under high heat flux from the internal surface to the external surface. Although finite element analysis and resonant ultrasound spectroscopy measurement indicated that the combined neutron irradiation and high heat flux gave rise to a high probability of cracking within the sample, the hermeticity evaluation of the tested sample still exhibited gas tightness, emphasizing that SiC cracking is inherently a statistical phenomenon. The developed permeation testing station is capable of measuring the gas permeation flux in the range of interest with full confidence based on the presented results. It is considered a critical pre- /post-irradiation examination technique to characterize SiC-based cladding materials in asreceived and irradiated states to aid the research and development of ATF.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Low-Power, Flexible Sensor Arrays with Solderless Board-to-Board Connectors for Monitoring Soil Deformation and Temperature

Landslides are a global and frequent natural hazard, affecting many communities and infrastructure networks. Technological solutions are needed for long-term, large-scale condition monitoring of infrastructure earthworks or natural slopes. However, current instruments for slope stability monitoring are often costly, require a complex installation process and/or data processing schemes, or have poor resolution. Wireless sensor networks comprising low-power, low-cost sensors have been shown to be a crucial part of landslide early warning systems. Here, we present the development of a novel sensing approach that uses linear arrays of three-axis accelerometers for monitoring changes in sensor inclination, and thus the surrounding soil’s deformation. By combining these deformation measurements with depth-resolved temperature measurements, we can link our data to subsurface thermal–hydrological regimes where relevant. In this research, we present a configuration of cascaded I2C sensors that (i) have ultra-low power consumption and (ii) enable an adjustable probe length. From an electromechanical perspective, we developed a novel board-to-board connection method that enables narrow, semi-flexible sensor arrays and a streamlined assembly process. The low-cost connection method relies on a specific FR4 printed circuit board design that allows board-to-board press fitting without using electromechanical components or solder connections. The sensor assembly is placed in a thin, semi-flexible tube (inner diameter 6.35 mm) that is filled with an epoxy compound. The resulting sensor probe is connected to an AA-battery-powered data logger with wireless connectivity. We characterize the system’s electromechanical properties and investigate the accuracy of deformation measurements. Our experiments, performed with probes up to 1.8 m long, demonstrate long-term connector stability, as well as probe mechanical flexibility. Furthermore, our accuracy analysis indicates that deformation measurements can be performed with a 0.390 mm resolution and a 95% confidence interval of ±0.73 mm per meter of probe length. This research shows the suitability of low-cost accelerometer arrays for distributed soil stability monitoring. In comparison with emerging low-cost measurements of surface displacement, our approach provides depth-resolved deformation, which can inform about shallow sliding surfaces.

deformation monitoring↗

Baseline Computational Fluid Dynamics Methodology for Longitudinal-Mode Liquid-Propellant Rocket Combustion Instability

A computational method for the analysis of longitudinal-mode liquid rocket combustion instability has been developed based on the unsteady, quasi-one-dimensional Euler equations where the combustion process source terms were introduced through the incorporation of a two-zone, linearized representation: (1) A two-parameter collapsed combustion zone at the injector face, and (2) a two-parameter distributed combustion zone based on a Lagrangian treatment of the propellant spray. The unsteady Euler equations in inhomogeneous form retain full hyperbolicity and are integrated implicitly in time using second-order, high-resolution, characteristic-based, flux-differencing spatial discretization with Roe-averaging of the Jacobian matrix. This method was initially validated against an analytical solution for nonreacting, isentropic duct acoustics with specified admittances at the inflow and outflow boundaries. For small amplitude perturbations, numerical predictions for the amplification coefficient and oscillation period were found to compare favorably with predictions from linearized small-disturbance theory as long as the grid exceeded a critical density (100 nodes/wavelength). The numerical methodology was then exercised on a generic combustor configuration using both collapsed and distributed combustion zone models with a short nozzle admittance approximation for the outflow boundary. In these cases, the response parameters were varied to determine stability limits defining resonant coupling onset.

Litchford, R. J.↗

Quasi-linear modeling of lower hybrid current drive in ITER and DEMO

First pass absorption of the Lower Hybrid waves in thermonuclear devices like ITER and DEMO is modeled by coupling the ray tracing equations with the quasi-linear evolution of the electron distribution function in 2D velocity space. As usually assumed, the Lower Hybrid Current Drive is not effective in a plasma of a tokamak fusion reactor, owing to the accessibility condition which, depending on the density, restricts the parallel wavenumber to values greater than n{sub ∥crit} and, at the same time, to the high electron temperature that would enhance the wave absorption and then restricts the RF power deposition to the very periphery of the plasma column (near the separatrix). In this work, by extensively using the “ray{sup star}” code, a parametric study of the propagation and absorption of the LH wave as function of the coupled wave spectrum (as its width, and peak value), has been performed very accurately. Such a careful investigation aims at controlling the power deposition layer possibly in the external half radius of the plasma, thus providing a valuable aid to the solution of how to control the plasma current profile in a toroidal magnetic configuration, and how to help the suppression of MHD mode that can develop in the outer part of the plasma. This analysis is useful not only for exploring the possibility of profile control of a pulsed operation reactor as well as the tearing mode stabilization, but also in order to reconsider the feasibility of steady state regime for DEMO.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Buoyancy instability of homologous implosions

With this study, I consider the hydrodynamic stability of imploding ideal gases as an idealized model for inertial confinement fusion capsules, sonoluminescent bubbles and the gravitational collapse of astrophysical gases. For oblate modes (short-wavelength incompressive modes elongated in the direction of the mean flow), a second-order ordinary differential equation is derived that can be used to assess the stability of any time-dependent flow with planar, cylindrical or spherical symmetry. Upon further restricting the analysis to homologous flows, it is shown that a monatomic gas is governed by the Schwarzschild criterion for buoyant stability. Under buoyantly unstable conditions, both entropy and vorticity fluctuations experience power-law growth in time, with a growth rate that depends upon mean flow gradients and, in the absence of dissipative effects, is independent of mode number. If the flow accelerates throughout the implosion, oblate modes amplify by a factor (2C) |N 0 |t i , where C is the convergence ratio of the implosion, N 0 is the initial buoyancy frequency and t i is the implosion time scale. If, instead, the implosion consists of a coasting phase followed by stagnation, oblate modes amplify by a factor exp(π|N 0 |t s ), where N 0 is the buoyancy frequency at stagnation and t s is the stagnation time scale. Even under stable conditions, vorticity fluctuations grow due to the conservation of angular momentum as the gas is compressed. For non-monatomic gases, this additional growth due to compression results in weak oscillatory growth under conditions that would otherwise be buoyantly stable; this over-stability is consistent with the conservation of wave action in the fluid frame. The above analytical results are verified by evolving the complete set of linear equations as an initial value problem, and it is demonstrated that oblate modes are the fastest-growing modes and that high mode numbers are required to reach this limit (Legendre mode ℓ ≳ 100 for spherical flows). Finally, comparisons are made with a Lagrangian hydrodynamics code, and it is found that a numerical resolution of ~30 zones per wavelength is required to capture these solutions accurately. This translates to an angular resolution of ~(12/ℓ)°, or ≲ 0.1° to resolve the fastest-growing modes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Detecting Permafrost Active Layer Thickness Change From Nonlinear Baseflow Recession

We report permafrost underlies about one fifth of the global land area and affects ground stability, freshwater runoff, soil chemistry, and surface-atmosphere gas exchange. The depth of thawed ground overlying permafrost (active layer thickness) has broadly increased across the Arctic in recent decades, coincident with a period of increased streamflow, especially the lowest flows (baseflow). Mechanistic links between active layer thickness and baseflow have recently been explored using linear reservoir theory, but most watersheds behave as nonlinear reservoirs. We derive theoretical nonlinear relationships between long-term average saturated soil thickness $\bar{η}$(proxy for active layer thickness) and long-term average baseflow. When applied to 38 years of daily streamflow data for the Kuparuk River basin on the North Slope of Alaska, the theory predicts $\bar{η}$ increased 0.17 ± 0.22[2σ] cm a -1 between 1983 and 2020 (6.4 ± cm total). The rate of increase nearly doubled to 0.29 ± 0.31 cm a -1 between 1990 and 2020, during which time local field measurements from Circumpolar Active Layer Monitoring sites indicate the active layer increased 0.31 ± 0.22 cm a -1 . The predicted rate of increase more than doubled again between 2002 and 2020, outpacing a near doubling of observed active layer thickening, consistent with trends in terrestrial water storage inferred from Gravity Recovery and Climate Experiment satellite gravimetry and Modern-Era Retrospective Analysis for Research and Applications climate reanalysis. Overall, hydrologic change is accelerating in the Kuparuk River basin, and we provide a theoretical framework for estimating basin-scale changes in active layer water storage from streamflow measurements.

54 ENVIRONMENTAL SCIENCES↗