Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Theoretical description and modelling”

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 55 records · Page 3

Motion of Drops on Surfaces with Wettability Gradients

A liquid drop present on a solid surface can move because of a gradient in wettability along the surface, as manifested by a gradient in the contact angle. The contact angle at a given point on the contact line between a solid and a liquid in a gaseous medium is the angle between the tangent planes to the liquid and the solid surfaces at that point and is measured within the liquid side, by convention. The motion of the drop occurs in the direction of increasing wettability. The cause of the motion is the net force exerted on the drop by the solid surface because of the variation of the contact angle around the periphery. This force causes acceleration of an initially stationary drop, and leads to its motion in the direction of decreasing contact angle. The nature of the motion is determined by the balance between the motivating force and the resisting hydrodynamic force from the solid surface and the surrounding gaseous medium. A wettability gradient can be chemically induced as shown by Chaudhury and Whitesides who provided unambiguous experimental evidence that drops can move in such gradients. The phenomenon can be important in heat transfer applications in low gravity, such as when condensation occurs on a surface. Daniel et al have demonstrated that the velocity of a drop on a surface due to a wettability gradient in the presence of condensation can be more than two orders of magnitude larger than that observed in the absence of condensation. In the present research program, we have begun to study the motion of a drop in a wettability gradient systematically using a model system. Our initial efforts will be restricted to a system in which no condensation occurs. The experiments are performed as follows. First, a rectangular strip of approximate dimensions 10 x 20 mm is cut out of a silicon wafer. The strip is cleaned thoroughly and its surface is exposed to the vapor from an alkylchlorosilane for a period lasting between one and two minutes inside a desiccator. This is done using an approximate line source of the vapor in the form of a string soaked in the alkylchlorosilane. Ordinarily, many fluids, including water, wet the surface of silicon quite well. This means that the contact angle is small. But the silanized surface resists wetting, with contact angles that are as large as 100 degs. Therefore, a gradient of wettability is formed on the silicon surface. The region near the string is highly hydrophobic, and the contact angle decreases gradually toward a small value at the hydrophilic end away from this region. The change in wettability occurs over a distance of several mm. The strip is placed on a platform within a Plexiglas cell. Drops of a suitable liquid are introduced on top of the strip near the hydrophobic end. An optical system attached to a video camera is trained on the drop so that images of the moving drop can be captured on videotape for subsequent analysis. We have performed preliminary experiments with water as well as ethylene glycol drops. Results from these experiments will be presented in the poster. Future plans include the refinement of the experimental system so as to permit images to be recorded from the side as well as the top, and the conduct of a systematic study in which the drop size is varied over a good range. Experiments will be conducted with different fluids so as to obtain the largest possible range of suitably defined Reynolds and Capillary numbers. Also, an effort will be initiated on theoretical modeling of this motion. The challenges in the development of the theoretical description lie in the proper analysis of the region in the vicinity of the contact line, as well as in the free boundary nature of the problem. It is known that continuum models assuming the no slip condition all the way to the contact line fail by predicting that the stress on the solid surface becomes singular as the contact line is approached. One approach for dealing with this issue has been to relax the no-slip boundary condition using the Navier model. Molecular dynamics simulations of the contact line region show that for a non-polar liquid on a solid surface, the no-slip boundary condition is in fact incorrect near the contact line. Furthermore, the same simulations also show that the usual relationship between stress and the rate of deformation breaks down in the vicinity of the contact line. In developing continuum theoretical models of the system, we shall accommodate this knowledge to the extent possible.

Subramanian, R. Shankar↗

Interfacing electron and neutrino quasielastic scattering cross sections with the spectral function in GENIE

Progress in neutrino-nucleus cross section models is being driven by the need for highly accurate predictions for the neutrino oscillation community. These sophisticated models are being developed within a microscopic description of the nucleus with the goal of encompassing all reaction modes relevant for the accelerator neutrino program. The disconnect between these microscopic models and the event generators that will be used in the next generation of experiments represents a critical obstacle that must be overcome in order to precisely measure the neutrino oscillation parameters. To this end we have developed a hadron tensor interface for lepton-nucleus quasielastic (QE) scattering within the GENIE event generator as a proof of principle, with the broader goal of creating an efficient pipeline for incorporating advanced theoretical models in event generators. As a demonstration of this interface we have implemented the spectral function model into GENIE by connecting theorist provided fortran code through the hadron tensor interface. The spectral function model offers a more complete description of the nuclear ground state, as well as the ability to provide quantifiable theoretical uncertainties. Finally, we validate this implementation and compare its predictions against data and against QE models already available in GENIE.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

SINDA-NASTRAN interfacing program theoretical description and user's manual

The task of converting SINDA finite difference thermal model temperature results into NASTRAN finite element model thermal loads can be very labor intensive if there is not one node-to-one element, or systematic node-to-element. correlation between models. This paper describes the SINDA-NASTRAN Interfacing Program (SNIP), a FORTRAN computer code that generates NASTRAN structural model thermal load cards given by SINDA (or similar thermal model) temperature results and thermal model geometric data. SNIP generates NASTRAN thermal load cards for NASTRAN plate, shell, bar, and beam elements. The paper describes the interfacing procedures used by SNIP, and discusses set-up and operation of the program. Sample cases are included to demonstrate use of the program and show its performance under a variety of conditions. SNIP can provide structural model thermal loads that accurately reflect thermal model results while reducing the time required to interface thermal and structural models when compared to other methods.

Winegar, Steven R.↗

User’s Manual for RESRAD-RDD&IND Code Version 2: Vol. 2—User’s Guide for RESRAD-RDD&IND Code

Version 2.0 of the RESRAD-RDD&IND computer code is designed to support the implementation of protective action guides (PAGs) after a nuclear emergency incident including a radiological dispersal device (RDD) and/or an improvised nuclear device (IND) incident (EPA 2017). Eight different group types, addressing various decisions, are available for selection. The RESRAD-RDD&IND code calculates radiological doses, stay times, etc., for the selected group that the user wishes to focus on. (That is, the results for all the groups are not calculated simultaneously, and the input for those other groups do not matter, although some parameter values are shared between groups.) Version 2.0 has a user-friendly interface so that the RESRAD-RDD&IND code can be used with minimal training. For example, the user can select the major characteristics of the problem-event type, source term, and decision type from the left side of the interface and then calculate the results with the default assumptions for the exposure scenarios. More in-depth analysis would include specifying site-specific exposure scenario characteristics in the right side of the interface. The procedures for data entry and results viewing are self-explanatory. This is because common window maneuvering features and text instructions were incorporated in the interface design. General and context-specific help are available to aid users entering parameter values, as well. The RESRAD-RDD&IND computer code gives the user the option to select either an RDD or IND incident for analysis. For an RDD event analysis, 11 radionuclides (Am-241, Cf-252, Cm-244, Co-60, Cs-137, Ir-192, Po-210, Pu-238, Pu-239, Ra-226, and Sr-90) are included. These 11 radionuclides are the radionuclides most likely used for an RDD. More than 90 radionuclides can be selected for an IND event analysis. Initial default concentrations are provided for 44 radionuclides for a uranium-fueled IND event. These 44 radionuclides are those that would contribute significantly to the radiation dose associated with a uranium-fueled bomb detonation. The radionuclides generated from ingrowth of these 44 initial radionuclides are also automatically included in the analysis. Pu-239, Cs-134m, Ru-105, and Rb-89 and their progeny can be selected for analysis if they are detected and their concentrations are determined. This user’s guide, which is Volume 2 of the User’s Manual for RESRAD-RDD&IND Code Version 2, provides instructions to users on how to install the RESRAD-RDD&IND code, navigate the interface, and use the various features, including those discussed above, to set up an analysis and view/print the results in text outputs. Volume 1 of the User’s Manual for RESRAD-RDD&IND Code Version 2 (Yu et al. 2026), which contains descriptions of the methodology and theoretical basis for dose modeling and the mathematical equations implemented in the code, can be accessed and viewed through the Help menu in the code or can be downloaded from the RESRAD website (https://resrad.evs.anl.gov).

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Laser amplification of incoherent radiation

The amplification of noise in a laser amplifier is treated theoretically. The model for the active medium and its description using density-matrix techniques, are taken from the theory of laser operation. The spectral behavior of the radiation in the nonlinear regime is studied and the formalism is written from the onset in the frequency domain. The statistics of the light are gradually modified by the nonlinear amplification process, and expressions are derived for the rate of change of fluctuations in intensity as a measure of statistical changes. In addition, the range of validity of Litvak's Gaussian-statistics approximation is discussed. In the homogeneous-broadening case, the evolution of initially broadband Gaussian radiation toward quasimonochromatic oscillations with laserlike statistics is explored in several numerical examples. The connections of this study with the time-domain work on self-pulsing in a ring-laser configuration, are established. Finally, spectral-narrowing and -rebroadening effects in Doppler-broadened media are discussed both analytically and with numerical examples. These examples show the distinct contribution of pulsations in the population ('Raman-type terms'), and saturation phenomena.

Menegozzi, L. N.↗

Turbulent heating of colliding streams in the solar wind.

Turbulent heating of colliding plasma streams has previously been observed in the solar wind. The original data were interpreted in terms of a fluid model. It is contended that a plasma-kinetic description is the more appropriate theoretical approach and is necessary in order to better understand the microscopic physical phenomena that underlie all fluid models. Microscopic solar-wind parameters characteristic of conditions during the observations were used, together with the quasi-linear plasma-kinetic theory, to compute the expected magnetic field and temperature enhancements in the interaction region between two counterstreaming plasma beams. The physical mechanism of excitation is the electromagnetic two-stream instability in which Alfven waves are unstable. A total field in the interaction region of about 8 gamma and a change in temperature of about 100,000 K are obtained.

Goldstein, M. L.↗

The overshoot region at the bottom of the solar convection zone

The extent and thermal stratification of the region of convective overshoot underneath the convection zone of the sun are investigated. The phenomenon of convective overshoot in general is discussed, and some of the modal and model approaches to studying it are briefly reviewed. A detailed theoretical description of the motion of plumes in a stably stratified medium is given, leading to a 'derivation' of the plume equations from the hydrodynamic equations. Entrainment is discussed, and it is shown how the plume equations can be used to compute convective overshoot in the sun. The limitations of the plume model are addressed, arguing that a thin boundary layer must exist which separates convective and radiative regions. The results of numerical integrations of the plume equations, as applied to the region of convective overshoot underneath the solar convective zone, are discussed.

Schmitt, J. H. M. M.↗

Theoretical description of proton-deuteron interactions using exact two-body dynamics of the femtoscopic correlation method

Modeling proton-deuteron interactions is particularly challenging. Due the deuteron's large size, the interaction can extend over several femtometers. The degree to which it can be modeled as a two-body problem might also be questioned. One way to study these interactions is through femtoscopic correlation measurements of particle pairs, extracting information using available theoretical models. In this work, we examine two approaches for describing proton-deuteron correlations: the Lednický-Lyuboshits formalism and full numerical solutions of the Schrödinger equation. Here, our results show that the differences between these methods are significant. Furthermore, we demonstrate that incorporating higher-order partial waves—particularly the p wave—is essential for accurately capturing the dynamics of proton–deuteron interactions and the full potential of the strong force.

Nucleon induced nuclear reactions↗

Development of Parametric Mass and Volume Models for an Aerospace SOFC/Gas Turbine Hybrid System

In aerospace power systems, mass and volume are key considerations to produce a viable design. The utilization of fuel cells is being studied for a commercial aircraft electrical power unit. Based on preliminary analyses, a SOFC/gas turbine system may be a potential solution. This paper describes the parametric mass and volume models that are used to assess an aerospace hybrid system design. The design tool utilizes input from the thermodynamic system model and produces component sizing, performance, and mass estimates. The software is designed such that the thermodynamic model is linked to the mass and volume model to provide immediate feedback during the design process. It allows for automating an optimization process that accounts for mass and volume in its figure of merit. Each component in the system is modeled with a combination of theoretical and empirical approaches. A description of the assumptions and design analyses is presented.

Tornabene, Robert↗

Majorana modes with side features in magnet-superconductor hybrid systems

Magnet-superconductor hybrid (MSH) systems represent promising platforms to host Majorana zero modes (MZMs), the elemental building blocks for fault-tolerant quantum computers. Theoretical description of such MSH structures is mostly based on simplified models, not accounting for the complexity of real materials. Here, based on density functional theory, we derive a superconducting 80-band model to study an MSH system consisting of a magnetic manganese chain on the s wave superconductor niobium. For a wide range of values of the superconducting order parameter, the system is a topological superconductor, with MZMs exhibiting non-universal spatial patterns and a drastic accumulation of spectral weight on both sides along the magnetic chain. These side feature states can be explained by an effective model which is guided by the ab initio results. Performing scanning tunneling spectroscopy experiments on the same system, we observe a spatial structure in the low-energy local density of states that is consistent with the theoretical findings. Our results open a first-principle approach to the discovery of topological superconductors.

36 MATERIALS SCIENCE↗

Unified theoretical framework for black carbon mixing state allows greater accuracy of climate effect estimation

Black carbon (BC) plays an important role in the climate system because of its strong warming effect, yet the magnitude of this effect is highly uncertain owing to the complex mixing state of aerosols. Here we build a unified theoretical framework to describe BC’s mixing states, linking dynamic processes to BC coating thickness distribution, and show its self-similarity for sites in diverse environments. The size distribution of BC-containing particles is found to follow a universal law and is independent of BC core size. A new mixing state module is established based on this finding and successfully applied in global and regional models, which increases the accuracy of aerosol climate effect estimations. Our theoretical framework links observations with model simulations in both mixing state description and light absorption quantification.

54 ENVIRONMENTAL SCIENCES↗

Empirical investigation of nuclear correlation function distributions in lattice QCD

Two-point correlation functions of systems with baryon number 𝐵 ∈ {1,2,3,4} are investigated using lattice quantum chromodynamics (QCD). In particular, the empirical distributions of importance-sampling Monte-Carlo samples of these correlation functions are examined as a function of the spacetime separation between the two points and the baryon number. While the exact forms of these distributions are not known for QCD, recent work has determined asymptotic expressions for analogous correlation function distributions in simpler theories such as scalar field theory and the disordered phase of the 𝑂⁡(𝑁) model. The theoretical 𝑂⁡(𝑁) model distributions are found to provide an accurate description of the empirical QCD distributions at zero momentum over a wide range of temporal separations for each baryon number when assessed with a range of different statistical tests. In particular, the behavior of the baryon number 𝐵 QCD correlation function at large temporal separation is well reproduced by the 𝑂⁡(𝑁 ∼ 2/𝐵) model distribution.

Lattice field theory↗

Reentrant Correlated Insulators in Twisted Bilayer Graphene at 25 T ( 2 π Flux)

Twisted bilayer graphene (TBG) is remarkable for its topological flat bands, which drive strongly interacting physics at integer fillings, and its simple theoretical description facilitated by the Bistritzer-MacDonald Hamiltonian, a continuum model coupling two Dirac fermions. Because of the large moiré unit cell, TBG offers the unprecedented opportunity to observe reentrant Hofstadter phases in laboratory-strength magnetic fields near 25 T. This Letter is devoted to magic angle TBG at 2π flux where the magnetic translation group commutes. Here we use a newly developed gauge-invariant formalism to determine the exact single-particle band structure and topology. We find that the characteristic TBG flat bands reemerge at 2π flux, but, due to the magnetic field breaking C 2z T, they split and acquire Chern number ±1. We show that reentrant correlated insulating states appear at 2π flux driven by the Coulomb interaction at integer fillings, and we predict the characteristic Landau fans from their excitation spectrum.

36 MATERIALS SCIENCE↗

The Role of Flow Diagnostic Techniques in Fan and Open Rotor Noise Modeling

A principal source of turbomachinery noise is the interaction of the rotating and stationary blade rows with the perturbations in the airstream through the engine. As such, a lot of research has been devoted to the study of the turbomachinery noise generation mechanisms. This is particularly true of fan and open rotors, both of which are the major contributors to the overall noise output of modern aircraft engines. Much of the research in fan and open rotor noise has been focused on developing theoretical models for predicting their noise characteristics. These models, which run the gamut from the semi-empirical to fully computational ones, are, in one form or another, informed by the description of the unsteady flow-field in which the propulsors (i.e., the fan and open rotors) operate. Not surprisingly, the fidelity of the theoretical models is dependent, to a large extent, on capturing the nuances of the unsteady flowfield that have a direct role in the noise generation process. As such, flow diagnostic techniques have proven to be indispensible in identifying the shortcoming of theoretical models and in helping to improve them. This presentation will provide a few examples of the role of flow diagnostic techniques in assessing the fidelity and robustness of the fan and open rotor noise prediction models.

Noise Prediction↗

Observation of Coherent 𝜙⁡(1020) Meson Photoproduction in Ultraperipheral PbPb Collisions at $\sqrt{𝑠_{NN}}$ = 5.36 TeV

The first observation of coherent 𝜙⁡(1020) meson photoproduction off heavy nuclei is presented using ultraperipheral lead-lead collisions at a center-of-mass energy per nucleon pair of 5.36 TeV. The data were collected by the CMS experiment and correspond to an integrated luminosity of 1.62 μ⁢𝑏 −1 . The 𝜙⁡(1020) meson signals are reconstructed via the 𝐾 + ⁢𝐾 − decay channel. The production cross section is presented as a function of the 𝜙⁡(1020) meson rapidity in the range 0.3 < |𝑦| < 1.0, probing gluons that carry a fraction of the nucleon momentum (𝑥) around 10 −4 . The observed cross section exhibits little dependence on rapidity and is significantly suppressed, by a factor of ∼5, compared to a baseline model that treats a nucleus as a collection of free nucleons. Theoretical models that incorporate the nuclear shadowing effect generally provide a better description of the 𝜙⁡(1020) data than those incorporating gluon saturation. This study establishes a powerful new tool for exploring nuclear effects and nuclear gluonic structure in the small-𝑥 regime at a unique energy scale bridging the perturbative and nonperturbative quantum chromodynamics domains.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Computational interferometric description of nested flow fields

Computer graphics and theoretical descriptions of density are used to obtain computer generated flow visualizations called computational interferograms. Computational interferograms are pictorially analogous to optical interferograms, and examples showing the fringe pattern for the flow about a sharp tip cone in a supersonic air stream are presented. To ascertain the effect of unsteady behavior, local density disturbances are added to the steady state flow field. This introduces irregularities to the computational interferogram like those seen in the optical interferograms. These theoretical disturbances can be varied in geometry, density description, translated with time, and strengthened or dissipated. The accuracy of computational interferometry relies on the accuracy of the theoretical density descriptions and therefore, it provides a way of verifying existing models of flow fields, especially those containing unsteady or turbulent behavior. In addition to being a unique method of flow visualization, computational interferometry can be used to develop and modify theories or numerical solutions to both simple and complex flow fields. The presented research is a general description of this process.

Havener, A. George↗

Hydraulic Actuator System for Rotor Control

In the last ten years, several different types of actuators were developed and fabricated for active control of rotors. A special hydraulic actuator system capable of generating high forces to rotating shafts via conventional bearings is addressed. The actively controlled hydraulic force actuator features an electrohydraulic servo valve which can produce amplitudes and forces at high frequencies necessary for influencing rotor vibrations. The mathematical description will be given in detail. The experimental results verify the theoretical model. Simulations already indicate the usefulness of this compact device for application to a real rotor system.

Ulbrich, Heinz↗

Landing impact of seaplanes

The theory of landing impact is briefly stated and the applicability of a previously suggested formula is extended. Theoretical considerations regarding impact measurements on models and actual seaplanes are followed by a brief description of the instruments used in actual flight tests. The report contains a description of the strength conditions and deals exhaustively with force measurements on the float gear of an "HE 9a" with flat-bottom and with V-bottom floats. The experimental data are given and compared with the theoretical results.

Pabst, Wilhelm↗