Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “analytical model”

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,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Thin shell formation in radiative shocks. 1: Supernova remnants in low-density media

This paper explores the onset of thin-shell formation in interstellar shocks associated with supernova explosions. We outline a simple but useful scheme that indicates the time at which thin shell formation begins for supernova remnants (SNRs) evolving in a range of interstellar environments, extending the previous analytical models to arbitrary power-law density media. The result depends on the gas cooling properties and the shock velocity and radius. This is then applied to the specific case of SNRs in low-density media. The procedure for defining the time for the onset of shell formation, t(sub sf), equates the value of the adiabat, kappa = p/rho(exp gamma), to zero using the known time dependence of the shock radius and velocity. For the case of a power-law density ambient medium of the form rho(r) = Br(exp -omega), it is found that shell formation can be prevented when the ambient density drops faster than a critical rate. For a cooling function of the form Lambda = Lambda(sub 0) tau(exp beta), with beta = -0.5 (appropriate for line cooling), shell formation never occurs for omega greater than or equal to 9/5. The shell formation time is then computed for spherical shocks in a power-law density medium. For omega = 0, the onset of shell formation is found to be at t(sub sf) approx. equal to 2.87 x 10(exp 4) E51(exp 3/14) n(sub 0 exp -4.7) yr, which agrees well with previous estimates derived by other means. We compare the analytical shell formation time with the results of detailed numerical models for omega = 0 and three different ambient densities and find good agreement. The extension of the criterion for the onset of thin shell formation using the ratio of cooling to swept-up column density is also described. This method provides a useful approximation for cases when the exact solution is not known.

Franco, Jose↗
NASA NTRS1981

A groundwater convection model for Rio Grande rift geothermal resources

It has been proposed that forced convection, driven by normal groundwater flow through the interconnected basins of the Rio Grande rift is the primary source mechanism for the numerous geothermal anomalies along the rift. A test of this concept using an analytical model indicates that significant forced convection must occur in the basins even if permeabilities are as low as 50-200 millidarcies at a depth of 2 km. Where groundwater flow is constricted at the discharge areas of the basins forced convection can locally increase the gradient to a level where free convection also occurs, generating surface heat flow anomalies 5-15 times background. A compilation of groundwater data for the rift basins shows a strong correlation between constrictions in groundwater flow and hot springs and geothermal anomalies, giving strong circumstantial support to the convection model.

Morgan, P.↗
NASA NTRS1999

Modeling and Validation of a Navy A6-Intruder Actively Controlled Landing Gear System

Concepts for long-range air travel are characterized by airframe designs with long, slender, relatively flexible fuselages. One aspect often overlooked is ground-induced vibration of these aircraft. This paper presents an analytical and experimental study of reducing ground-induced aircraft vibration loads by using actively controlled landing gear. A facility has been developed to test various active landing gear control concepts and their performance, The facility uses a Navy A6 Intruder landing gear fitted with an auxiliary hydraulic supply electronically controlled by servo valves. An analytical model of the gear is presented, including modifications to actuate the gear externally, and test data are used to validate the model. The control design is described and closed-loop test and analysis comparisons are presented.

Horta, Lucas G.↗
DOE OSTI2020

Integrated modelling: Coupling of surface evolution and plasma-impurity transport

During the interaction of the scrape off layer (SOL) plasma with the first wall the evolution of both wall and plasma are tightly coupled: The erosion of the first wall leads to an impurity concentration in the plasma which affects the particle and power balance in the plasma. In turn the impurities, when leaving the plasma via transport, can form deposits and mixed materials, far away from their initial source location. These deposits can be eroded, allowing the impurities to stepwise migrate through the fusion device until they end up at a location where the plasma at the wall is cold enough and no further erosion occurs. To describe these processes an integrated model of surface evolution and plasma transport of impurities is needed. The WallDYN code achieves this required coupling of processes by parameterising the output of surface evolution- and plasma impurity-migration-codes by analytical models. For a given fixed background plasma it evolves the surface composition, derives impurity flux into and from the plasma and can from this derive the impurity densities in the plasma. This paper will show the importance of including this recycling of impurities at the wall in impurity migration modelling: The 13 CH 4 seeding experiment performed in the Wendelstein 7-X Stellerator is modelled using the recent extension of the WallDYN code to 3D plasma and wall geometries. A comparison with post mortem analysis of the 13 C deposition shows both qualitative and quantitative agreement with the WallDYN calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Fast and Universal Kohn-Sham Density Functional Theory Algorithm for Warm Dense Matter to Hot Dense Plasma

Understanding many processes, e.g., fusion experiments, planetary interiors, and dwarf stars, depends strongly on microscopic physics modeling of warm dense matter and hot dense plasma. This complex state of matter consists of a transient mixture of degenerate and nearly free electrons, molecules, and ions. This regime challenges both experiment and analytical modeling, necessitating predictive ab initio atomistic computation, typically based on quantum mechanical Kohn-Sham density functional theory (KS-DFT). However, cubic computational scaling with temperature and system size prohibits the use of DFT through much of the warm dense matter regime. A recently developed stochastic approach to KS-DFT can be used at high temperatures, with the exact same accuracy as the deterministic approach, but the stochastic error can converge slowly and it remains expensive for intermediate temperatures (< 50 eV). Here we have developed a universal mixed stochastic-deterministic algorithm for DFT at any temperature. This approach leverages the physics of KS-DFT to seamlessly integrate the best aspects of these different approaches. We demonstrate that this method significantly accelerated self-consistent field calculations for temperatures from 3 to 50 eV, while producing stable molecular dynamics and accurate diffusion coefficients.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1984

A study of aeroelastic and structural dynamic effects in multi-rotor systems with application to hybrid heavy lift vehicles

An aeroelastic model suitable for the study of aeroelastic and structural dynamic effects in multirotor vehicles simulating a hybrid heavy lift vehicle was developed and applied to the study of a number of diverse problems. The analytical model developed proved capable of modeling a number of aeroelastic problems, namely: (1) isolated blade aeroelastic stability in hover and forward flight, (2) coupled rotor/fuselage aeromechanical problem in air or ground resonance, (3) tandem rotor coupled rotor/fuselage problems, and (4) the aeromechanical stability of a multirotor vehicle model representing a hybrid heavy lift airship (HHLA). The model was used to simulate the ground resonance boundaries of a three bladed hingeless rotor model, including the effect of aerodynamic loads, and the theoretical predictions compared well with experimental results. Subsequently the model was used to study the aeromechanical stability of a vehicle representing a hybrid heavy lift airship, and potential instabilities which could occur for this type of vehicle were identified. The coupling between various blade, supporting structure and rigid body modes was identified.

Friedmann, P. P.↗
NASA NTRS1990

Error analysis for pulsed maneuvers of a dual-spin spacecraft

The Galileo is a dual-spin interplanetary spacecraft which is scheduled for launch toward Jupiter in October of 1989. The spacecraft has twelve thrusters, located on the rotor, which can be operated in continuous or pulsed modes. In this paper probabilistic error models are presented for the pulsed maneuvers which are used to change the spacecraft's velocity. These models take into account all error sources which might affect the maneuver accuracies, such as wobble, nutation, plume impingement, thrust level variation and misalignment, gyro drift, burn timing errors and mass property uncertainties. The analytic models are sufficiently general so that they may be applied to all spin and dual-spin stabilized spacecraft. When these models are applied to the case of the Galileo spacecraft, the numerical results demonstrate that the spacecraft can achieve the maneuver accuracies required in order to successfully reach Jupiter and perform its scientific mission.

Longuski, James M.↗
NASA NTRS1999

Multi-Body Analysis of the 1/5 Scale Wind Tunnel Model of the V-22 Tiltrotor

The paper presents a multi-body analysis of the 1/5 scale wind tunnel model of the V-22 tiltrotor, the Wing and Rotor Aeroelastic Testing System (WRATS), currently tested at NASA Langley Research Center. An original multi-body formulation has been developed at the Dipartimento di Ingegneria Aerospaziale of the Politecnico di Milano, Italy. It is based on the direct writing of the equilibrium equations of independent rigid bodies, connected by kinematic constraints that result in the addition of algebraic constraint equations, and by dynamic constraints, that directly contribute to the equilibrium equations. The formulation has been extended to the simultaneous solution of interdisciplinary problems by modeling electric and hydraulic networks, for aeroservoelastic problems. The code has been tailored to the modeling of rotorcrafts while preserving a complete generality. A family of aerodynamic elements has been introduced to model high aspect aerodynamic surfaces, based on the strip theory, with quasi-steady aerodynamic coefficients, compressibility, post-stall interpolation of experimental data, dynamic stall modeling, and radial flow drag. Different models for the induced velocity of the rotor can be used, from uniform velocity to dynamic in flow. A complete dynamic and aeroelastic analysis of the model of the V-22 tiltrotor has been performed, to assess the validity of the formulation and to exploit the unique features of multi-body analysis with respect to conventional comprehensive rotorcraft codes; These are the ability to model the exact kinematics of mechanical systems, and the possibility to simulate unusual maneuvers and unusual flight conditions, that are particular to the tiltrotor, e.g. the conversion maneuver. A complete modal validation of the analytical model has been performed, to assess the ability to reproduce the correct dynamics of the system with a relatively coarse beam model of the semispan wing, pylon and rotor. Particular care has been used to model the kinematics of the gimbal joint, that characterizes the rotor hub, and of the control system, consisting in the entire swashplate mechanism. The kinematics of the fixed and the rotating plates have been modeled, with variable length control links used to input the controls, the rotating flexible links, the pitch horns and the pitch bearings. The investigations took advantage of concurring wind tunnel test runs, that were performed in August 1998, and allowed the acquisition of data specific to the multi-body analysis.

Ghiringhelli, G. L.↗
DOE OSTI2020

Nanostructured Tube Wakefield Accelerator

Unprecedented $\rm TeVm^{-1}$ acceleration gradients are modeled to be realizable using a nonlinear surface crunch-in mode in nanostructured tubes. This mode is realizable using advances in nanofabrication and solid energy density attosecond bunch compression. Three dimensional computational and analytical modeling demonstrates GeV energy gain in sub-millimeter long tubes with effective wall densities $n_{\rm t}\sim10^{22-24}\rm cm^{-3}$ and hundreds of nanometer core radius when driven by submicron near solid electron beams, $n_{\rm b}\sim0.05n_{\rm t}$. Besides the many $\rm TVm^{-1}$ average gradients, strong self-focusing and nanomodulation of the beam which increases its peak density and the wakefield strength also opens up controlled high-energy photon production.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Simulations of vertical displacement oscillatory modes and global Alfvén Eigenmodes in JET geometry

Vertical Displacement Oscillatory Modes (VDOM), with frequency in the Alfvén range, are natural modes of oscillation of magnetically confined laboratory plasmas with elongated cross-section. These axisymmetric modes arise from the interaction between the plasma current, which is in equilibrium with currents flowing in external coils, and perturbed currents induced on a nearby conducting wall. The restoring force exerted by these perturbed currents on the vertical motion of the plasma column leads to its oscillatory behavior. An analytic model for VDOM was proposed based on an idealized 'straight tokamak' equilibrium with uniform equilibrium current density. This article introduces the first numerical simulations of VDOM in a realistic JET tokamak configuration, using the extended-MHD code NIMROD and drawing comparisons with Global Alfvén Eigenmodes (GAE). The results show qualitative agreement with analytic predictions regarding mode frequency and radial structure, supporting the identification of VDOM as a fundamental oscillation mode in tokamak plasmas. VDOM and GAE are modeled in a representative JET discharge, where axisymmetric perturbations with toroidal mode number n = 0 driven unstable by fast ions were observed. The two modes are examined separately using a forced oscillator within the NIMROD code, which enables a comparison of their characteristics and helps identify the experimentally observed mode possibly as a GAE.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Experimental validation of multiphysics model simulations of the thermal response of a cement clinker rotary kiln at laboratory scale

Abstract An increasing demand for buildings, transportation systems and civil infrastructure development has driven expansion of cement consumption world‐wide, producing a significant increase in related global energy demand. With approximately 7% of the world‐wide industrial energy consumption (10.7 exajoules [EJ]), the cement industry is the third most energy intensive industrial processes and a key component for concrete, the most consumed composite material in the global construction industry. In cement manufacturing, the cement kiln accounts for most of the energy consumption in the production process. As the heart of a cement plant, the cement kiln is where the kiln feed primarily containing calcium oxide (CaO), silica (SiO 2 ), alumina (Al 2 O 3 ), and iron (Fe 2 O 3 ) are thermally and chemically transformed into clinker minerals. The presented work developed a multiphysics model, designed and built a laboratory‐scale rotary cement clinker kiln, and produced cement clinker at laboratory‐scale. The model was developed to study the interaction between the various thermal, fluid dynamic and chemical interactions involved in the sintering process used to form Portland cement clinker in an effort to reduce energy use. The analytical model was validated through experimental testing using a unique laboratory‐scale rotary cement kiln developed during the investigation. Also demonstrated was the feasibility of producing clinker at laboratory scale. This modeling and lab scale tests were designed to better understand the clinker sintering process so that operational and quality decisions can be made to optimize energy consumption without compromising cement clinker quality. The computational fluid dynamics modeling was developed in COMSOL Multiphysics 6.0. The characteristics of the combustion fluid flow, concentration of species, temperature and heat transfer were studied for a turbulent flow of methane (CH 4 ) gas and oxygen (O 2 ). Theory suggests that heat transfer impacts the cement production process but the multiphysics model more accurately describes the convection, conduction, and radiant heat transfer in the kilning process and thus allows for a better understanding of the energy exchange driving the chemical reactions that produce Portland cement. Clinker minerals were formed because of appropriate burning conditions implemented during experimental model validation.

Tabares, Juan David↗
DOE OSTI2021

Robustness of cement-based materials: From dosage variations to yield stress fluctuations

The idea behind the present work is to assess robustness of fresh cement-based materials (i.e. their ability to display the same properties when submitted to variations in the way they are produced) from a generic and analytical point of view. First, analytical models from literature relating components proportions and yield stress are combined into a multi-scale and physical approach and compared to our experimental yield stress measurements. In the second part, the derivation of these analytical relationships is used to assess robustness as a function of components proportions, which allows for the most sensitive variations to be mapped onto component dosage from a yield stress point of view. In the case of fluid concretes, this approach showcases the interest of working at slightly higher water content and slightly lower aggregates fraction along with benefits of viscosity modifying agents / plasticizers combination.

36 MATERIALS SCIENCE↗
NASA NTRS1980

Acoustoelasticity - Review and extension

Research on the interaction between the sound pressure field and the flexible wall of an enclosure is reviewed. The topics discussed include: a mathematical model, simplified solutions for sinusoidal excitation, damping, random excitation, numerical results, including comparisons with representative experimental data, and design analysis methodology. Results of various authors using finite element techniques and modal analysis suggest that effective analytical models are available.

Dowell, E. H.↗
NASA NTRS1995

Testing higher-order Lagrangian perturbation theory against numerical simulations. 2: Hierarchical models

We present results showing an improvement of the accuracy of perturbation theory as applied to cosmological structure formation for a useful range of scales. The Lagrangian theory of gravitational instability of Friedmann-Lemaitre cosmogonies is compared with numerical simulations. We study the dynamics of hierarchical models as a second step. In the first step we analyzed the performance of the Lagrangian schemes for pancake models, the difference being that in the latter models the initial power spectrum is truncated. This work probed the quasi-linear and weakly non-linear regimes. We here explore whether the results found for pancake models carry over to hierarchical models which are evolved deeply into the non-linear regime. We smooth the initial data by using a variety of filter types and filter scales in order to determine the optimal performance of the analytical models, as has been done for the 'Zel'dovich-approximation' - hereafter TZA - in previous work. We find that for spectra with negative power-index the second-order scheme performs considerably better than TZA in terms of statistics which probe the dynamics, and slightly better in terms of low-order statistics like the power-spectrum. However, in contrast to the results found for pancake models, where the higher-order schemes get worse than TZA at late non-linear stages and on small scales, we here find that the second-order model is as robust as TZA, retaining the improvement at later stages and on smaller scales. In view of these results we expect that the second-order truncated Lagrangian model is especially useful for the modelling of standard dark matter models such as Hot-, Cold-, and Mixed-Dark-Matter.

Melott, A. L.↗
DOE OSTI2022

Polariton relaxation under vibrational strong coupling: Comparing cavity molecular dynamics simulations against Fermi’s golden rule rate

Under vibrational strong coupling (VSC), the formation of molecular polaritons may significantly modify the photo-induced or thermal properties of molecules. In an effort to understand these intriguing modifications, both experimental and theoretical studies have focused on the ultrafast dynamics of vibrational polaritons. Here, following our recent work [Li et al., J. Chem. Phys. 154, 094124 (2021)], we systematically study the mechanism of polariton relaxation for liquid CO 2 under a weak external pumping. Classical cavity molecular dynamics (CavMD) simulations confirm that polariton relaxation results from the combined effects of (i) cavity loss through the photonic component and (ii) dephasing of the bright-mode component to vibrational dark modes as mediated by intermolecular interactions. The latter polaritonic dephasing rate is proportional to the product of the weight of the bright mode in the polariton wave function and the spectral overlap between the polariton and dark modes. Both these factors are sensitive to parameters such as the Rabi splitting and cavity mode detuning. Compared to a Fermi’s golden rule calculation based on a tight-binding harmonic model, CavMD yields a similar parameter dependence for the upper polariton relaxation lifetime but sometimes a modest disagreement for the lower polariton. We suggest that this disagreement results from polariton-enhanced molecular nonlinear absorption due to molecular anharmonicity, which is not included in our analytical model. Finally, we also summarize recent progress on probing nonreactive VSC dynamics with CavMD.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Revisiting a Drag Partition Model For Canopy-Like Roughness Elements

Turbulent flows over a large surface area (S) covered by n obstacles experience an overall drag due to the presence of the ground and the protruding obstacles into the flow. The drag partition between the roughness obstacles and the ground is analyzed using an analytical model proposed by Raupach (Boundary-Layer Meteorol 60:375-395, 1992) and is hereafter referred to as R92. The R92 is based on the premise that the wake behind an isolated roughness element can be described by a shelter area A and a shelter volume V. The individual sizes of A and V without any interference from other obstacles can be determined from scaling analysis for the spread of wakes. To upscale from an individual roughness element to n/S elements where wakes may interact, R92 adopted a background stress re-normalizing instead of reducing A or V with each element addition. This work demonstrates that R92’s approach results in a linear background stress reduction in A and V only when the ratio of n/S is small, due to a low probability of wake interactions. This probabilistic nature suggests that up-scaling from individual to multiple roughness elements can be re-formulated using stochastic averaging methods proposed here. The two approaches are shown to recover R92 under plausible conditions. An alternative scaling for the shelter volume is also proposed here using thermodynamic arguments of work and dissipation though the final outcome remains similar to R92. Comparisons between R92 and available data spanning more than two decades after R92 on blocks and vegetation-like roughness elements confirm the practical utility of R92. The agreement between R92 and this updated databases of experiments and simulations confirm the potential use of R92 in large-scale models provided that the relevant parameters accommodate certain features of the roughness element type (cube versus vegetation-like) and, to a lesser extent, their configuration throughout S. Last, a comparison between R92 and models based on first-order closure principles with constant mixing length suggests that R92 can outperform such models when evaluated across a wide range of roughness densities.

canopy turbulence↗
DOE OSTI2022

Machine-learning enabled thermodynamic model for the design of new rare-earth compounds

We employ a descriptor based machine-learning approach to assess the effect of chemical alloying on formation-enthalpy of rare-earth intermetallics. Application of machine-learning approaches in rare-earth intermetallic design have been sparse due to limited availability of reliable datasets. In this work, we developed an ‘in-house’ rare-earth database with more than 600 + compounds, each entry was populated with formation enthalpy and related atomic features using high-throughput density-functional theory (DFT). The SISSO (sure independence screening and sparsifying operator) based machine-learning method with meaningful atomic features was used for training and testing the formation enthalpies of rare earth compounds. The complex lattice function coupled with the machine-learning model was used to explore the effect of transition metal alloying on the energy stability of Ce based cubic Laves phases (MgCu 2 type). The SISSO predictions show good agreement with high-fidelity DFT calculations and X-ray powder diffraction measurements. Our study provides quantitative guidance for compositional considerations within a machine-learning model and discovering new metastable materials. The electronic-structure of Ce-Fe-Cu based compound was also analyzed to get an in-depth understanding of the electronic origin of phase stability. The interpretable analytical models in combination with density-functional theory and experiments provide a fast and reliable design guide for discovering technologically useful materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Towards a quantitative understanding of bonding in supersonic single particle impacts: A three-dimensional FIB-SEM exploration

Particle bonding is crucial to coating quality in cold spray, but it has been a challenge to accurately quantify bonding even in single particle impacts. This paper uses FIB-SEM to explicitly map the particle-substrate interface for Cu-on-Cu single microparticle impacts in a full 3D rendering that spans a wide range of impact velocities. This approach permits a detailed quantification of the total bonding area and all of its associated components. In addition to revealing why prior 2D characterization efforts have missed important details about impact bonding, these data quantitatively reveal the evolution of bonding from its onset at the “critical velocity” V cr (where bonding is generally poor, ∼6 %) to its peak at around 1.3‧V cr (where almost 90 % of the particle bonds). Further increase in the velocity to 1.5‧V cr and beyond finds the onset of hydrodynamic penetration and a decrease in bonding. These data then support the development of a simple analytical model based on oxide rarefication and extrusion of bare metal through gaps in the oxide layer as driving the development of bonding. As a result, the model reproduces the experiments and provides guidance on optimization of bonding as a function of material and process parameters.

Analytical methods↗