Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Direct numerical simulation”

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,621 records · Page 90
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Creep anisotropy modeling and uncertainty quantification of an additively manufactured Ni-based superalloy

The advantages offered by additive manufacturing over traditional processes has driven a great deal of industrial and academic interest in recent years. However, the process is relatively new and requires additional investigation to become sufficiently mature for wide scale industrial adoption. Electron beam melting powder bed fusion is one technology that has shown promise for fabricating high temperature resistant materials such as nickel based superalloys. The resulting microstructures typically exhibit a strong fiber texture in the build direction giving rise to anisotropic time-dependent deformation behavior. In order to accelerate the qualification of these materials for industrial adoption accurate numerical models are needed for simulating their behavior. In this work a crystal plasticity model including non-Schmid effects is presented for capturing creep anisotropy observed in additively manufactured IN738LC. The model is calibrated via a probabilistic framework where model parameters are treated as random variables. An iterative sequential design strategy is utilized to efficiently identify the probability density of the unknown model parameters. As a case study the model is utilized to investigate the behavior of randomly oriented equiaxed grain clusters sometimes observed embedded in the additively manufactured columnar structure. A synthetic realization is simulated and uncertainty is propagated through to the full-field response. Results indicate that these features are the source of significant creep relaxation and strain accumulation which partially explains observed grain boundary decohesion at these locations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Physics–Informed Neural Networks of the Saint–Venant Equations for Downscaling a Large–Scale River Model

Large-scale river models are being refined over coastal regions to improve the scientific understanding of coastal processes, hazards and responses to climate change. However, coarse mesh resolutions and approximations in physical representations of tidal rivers limit the performance of such models at resolving the complex flow dynamics near the river-ocean interface, resulting in inaccurate simulations of flood inundation. In this research, we propose a machine learning (ML) framework based on the state-of-the-art physics-informed neural network (PINN) to simulate the downscaled flow at the subgrid scale. First, we demonstrate that PINN is able to assimilate observations of various types and solve the one-dimensional (1-D) Saint-Venant equations (SVE) directly. We perform the flow simulations over a floodplain and along an open channel in several synthetic case studies. The PINN performance is evaluated against analytical solutions and numerical models. Our results indicate that the PINN solutions of water depth have satisfactory accuracy with limited observations assimilated. In the case of flood wave propagation induced by storm surge and tide, a new neural network architecture is proposed based on Fourier feature embeddings that seamlessly encodes the periodic tidal boundary condition in the PINN's formulation. Furthermore, we show that the PINN-based downscaling can produce more reasonable subgrid solutions of the along-channel water depth by assimilating observational data. The PINN solution outperforms the simple linear interpolation in resolving the topography and dynamic flow regimes at the subgrid scale. This study provides a promising path towards improving emulation capabilities in large-scale models to characterize fine-scale coastal processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Material and component developments for the DEMO divertor using fibre reinforcement and additive manufacturing

Within the research along the European Fusion Roadmap, water-cooled divertor PFCs are foreseen in the design of a first fusion demonstration power plant (DEMO) in order to provide reliable heat removal capability. In the frame of this concerted attempt, the Max Planck Institute for Plasma Physics is concentrating on the development and testing of composite materials based on tungsten (W, preferred armour material) and copper (Cu, preferred heat sink material). W fibres (W f ) as monofilaments and yarns as reinforcement play a central role in these investigations due to their extraordinary properties concerning ductility already at room temperature and high tensile strength. Recent investigations on the impact of radiation damage suggest that the fibres retain their ductility upon irradiation. W reinforced with W fibres (W f /W) allows to overcome the intrinsic brittleness of W. Quantitative mechanical fracture tests of W f /W confirm the basic mechanisms of fibre reinforcement and the increased resistance to mechanical fatigue. The good wettability of W with liquid Cu and the absence of any metallurgical solubility make up an ideal material pairing for composite production. W fibre-reinforced Cu (W f /Cu) cooling tubes provide a rather high thermal conductivity (> 250 W mK -1 ) and at least twice the strength of CuCrZr in hoop direction in the temperature range up to at least 500 °C. Very recent neutron irradiation experiments confirm the sustainment of ductility of the W f /Cu composite. Numerical simulations suggest that thermal stresses in W-Cu PFCs could be strongly reduced by tailoring the local W and Cu volume fraction. This 'freely' distributed material composition can be achieved by means of additively manufactured W skeletons consecutively infiltrated by Cu. Investigations with W preforms produced by Laser Beam Powder Bed Fusion and infiltrated by Cu demonstrate the feasibility of this approach while testing of specifically prepared specimen is ongoing.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Magnetic dipole emission in resonant metal-dielectric-metal structures

We explore a possibility to control magnetic dipole emission with plasmonic cavities, placing Eu 3+ emitters inside profile-modulated metal-dielectric-metal structures. Significant variations in the branching ratio of the magnetic and electric dipole transitions are observed as the function of the thickness of the intermediate layer. The experimental results are confirmed with numerical simulations which account for cavity and gap plasmon resonances and predict modifications in the spontaneous emission spectrum as the function of the gap size and a strong directionality of the emission for small thicknesses of the intermediate layer. The implications of having a competition between electric and magnetic dipole relaxation channels in Eu 3+ are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1985

Middle atmosphere tides

The recent direction of theoretical and observational research in middle atmosphere tides is reviewed. The diurnal tide has received considerable emphasis, including observations of evanescent components at high latitudes, numerical simulations of nonmigrating propagating components due to longitudinally varying water vapor insolation absorption in the troposphere, and implications and origins of unsteadiness in tidal oscillations over periods of days. Recent analyses of Nimbus 7 LIMS temperatures and partial reflection drift measurements of neutral winds in the tropical mesosphere have also emphasized the diurnal propagating tide. A description of the Atmospheric Tides Middle Atmosphere Program (ATMAP) and an overview of ATMAP campaign results are also presented.

Forbes, J. M.↗
NASA NTRS1985

Samics on the IBM PC

An updated version of the Solar Array Manufacturing Industry Costing Standards (SAMICS) was presented. This version will run on the IBM PC-XT or compatibles, and embodies user friendly input screens and numerous help options. Use of this microcomputer version will still produce the old main frame direct cost driven analysis, although at some cost in turnaround time. A typical simulation will take about 4 hours. The program allows for unattended report printout. Also available is the IBM PC (or compatibles) version of the Improved Price Estimation Guidelines (IPEG). This program allows rapid analysis of process variables using SAMICS-generated coefficients.

Chamberlain, R.↗
NASA NTRS1988

Design for performance enhancement in feedback control systems with multiple saturating nonlinearities

A systematic control design methodology is introduced for multi-input/multi-output systems with multiple saturations. The methodology can be applied to stable and unstable open loop plants with magnitude and/or rate control saturations and to systems in which state limitations are desired. This new methodology is a substantial improvement over previous heuristic single-input/single-output approaches. The idea is to introduce a supervisor loop so that when the references and/or disturbances are sufficiently small, the control system operates linearly as designed. For signals large enough to cause saturations, the control law is modified in such a way to ensure stability and to preserve, to the extent possible, the behavior of the linear control design. Key benefits of this methodology are: the modified compensator never produces saturating control signals, integrators and/or slow dynamics in the compensator never windup, the directional properties of the controls are maintained, and the closed loop system has certain guaranteed stability properties. The advantages of the new design methodology are illustrated by numerous simulations, including the multivariable longitudinal control of modified models of the F-8 (stable) and F-16 (unstable) aircraft.

Kapasouris, Petros↗
NASA NTRS2008

Explaining Warm Coronal Loops

One of the great mysteries of coronal physics that has come to light in the last few years is the discovery that warn (- 1 INK) coronal loops are much denser than expected for quasi-static equilibrium. Both the excess densities and relatively long lifetimes of the loops can be explained with bundles of unresolved strands that are heated impulsively to very high temperatures. Since neighboring strands are at different stages of cooling, the composite loop bundle is multi-thermal, with the distribution of temperatures depending on the details of the "nanoflare storm." Emission hotter than 2 MK is predicted, but it is not clear that such emission is always observed. We consider two possible explanations for the existence of over-dense warm loops without corresponding hot emission: (1) loops are bundles of nanoflare heated strands, but a significant fraction of the nanoflare energy takes the form of a nonthermal electron beam rather then direct plasma heating; (2) loops are bundles of strands that undergo thermal nonequilibrium that results when steady heating is sufficiently concentrated near the footpoints. We present numerical hydro simulations of both of these possibilities and explore the observational consequences, including the production of hard X-ray emission and absorption by cool material in the corona.

Klimchuk, James A.↗
DOE OSTI2024

New systems in MOOSE

The Multiphysics Object-Oriented Simulation Environment (MOOSE) serves as a common library of classes between applications developed for advanced reactor analysis, fusion device engineering, spent fuel cask analysis, geochemistry studies, among other fields. These applications drive the development of the framework to meet their needs. Systems in MOOSE group capabilities that share a common purpose and generally common code. They can be leveraged by all downstream applications, providing extensive code re-use and shared maintenance. They facilitate the discovery by new users of the classes meeting at least partially their needs, and offer the same opportunities for customization as other systems. The addition of a new system to MOOSE opens new ways of solving or discretizing nonlinear problems, of performing distributed postprocessing, and a plethora of other needs. While new systems can be introduced in downstream applications rather than at the framework level, the framework team monitors common needs across the community and often triggers their addition. Documentation, training material, development needs can be centralized, limiting duplicated work across the community. The last three years have seen a large expansion in the capabilities of MOOSE. The supporting role of the framework in the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program has created numerous feature requests to support neutronics, thermal hydraulics, computational fluid dynamics and thermo-mechanics simulations in the Griffin, SAM, Pronghorn and Bison applications respectively. Similarly, laboratory-directed research and development (LDRD) projects in additive manufacturing, high-Reynolds flow simulations, structure optimization also necessitate an expansion of the framework capabilities. This summary reports on the new systems created in MOOSE, their design, their capabilities and some of the relevant interfaces.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRSDate not supplied

Comparing 3D and 2D CFD for Mars Helicopter Ingenuity Rotor Performance Prediction

Single and coaxial rotor performance simulations for the Mars Helicopter Ingenuity rotor are performed for representative Mars atmospheric conditions. Analyses are presented using both a high-fidelity 3D CFD model of the rotor and 2D CFD models of the airfoil sections for comprehensive analyses that use CAMRADII (Comprehensive Analytical Model of Rotorcraft Aerodynamics and Dynamics). When available, the airfoil performance calculations are generated using a numerical approach identical to that used in the high-fidelity 3D model, allowing for a direct comparison between the approaches. Experimental data from a validation campaign to explore higher thrust from an Ingenuity rotor is provided to substantiate a discussion on the simulation fidelity required for both coaxial and single rotor performance predictions. The data is in support of the Sample Recovery Helicopter (SRH) element that serves as the primary backup for tube retrieval as part of the Mars Sample Return (MSR) Campaign. Insights on modeling turbulence at low Reynolds numbers and its influence on the rotor figure of merit are discussed. Key rotor performance metrics are compared. A detailed investigation into differences between 2D and 3D rotor performance predictions, spanwise loading, and rotor stall behavior is included.

Mars Helicopter↗
DOE OSTI2021

A Segregated Approach for Modeling the Electrochemistry in the 3-D Microstructure of Li-Ion Batteries and Its Acceleration Using Block Preconditioners

Abstract Battery performance is strongly correlated with electrode microstructure. Electrode materials for lithium-ion batteries have complex microstructure geometries that require millions of degrees of freedom to solve the electrochemical system at the microstructure scale. A fast-iterative solver with an appropriate preconditioner is then required to simulate large representative volume in a reasonable time. In this work, a finite element electrochemical model is developed to resolve the concentration and potential within the electrode active materials and the electrolyte domains at the microstructure scale, with an emphasis on numerical stability and scaling performances. The block Gauss-Seidel (BGS) numerical method is implemented because the system of equations within the electrodes is coupled only through the nonlinear Butler–Volmer equation, which governs the electrochemical reaction at the interface between the domains. The best solution strategy found in this work consists of splitting the system into two blocks—one for the concentration and one for the potential field—and then performing block generalized minimal residual preconditioned with algebraic multigrid, using the FEniCS and the Portable, Extensible Toolkit for Scientific Computation libraries. Significant improvements in terms of time to solution (six times faster) and memory usage (halving) are achieved compared with the MUltifrontal Massively Parallel sparse direct Solver. Additionally, BGS experiences decent strong parallel scaling within the electrode domains. Last, the system of equations is modified to specifically address numerical instability induced by electrolyte depletion, which is particularly valuable for simulating fast-charge scenarios relevant for automotive application.

25 ENERGY STORAGE↗
NASA NTRS1993

Simulation of wave propagation in three-dimensional random media

Quantitative error analysis for simulation of wave propagation in three dimensional random media assuming narrow angular scattering are presented for the plane wave and spherical wave geometry. This includes the errors resulting from finite grid size, finite simulation dimensions, and the separation of the two-dimensional screens along the propagation direction. Simple error scalings are determined for power-law spectra of the random refractive index of the media. The effects of a finite inner scale are also considered. The spatial spectra of the intensity errors are calculated and compared to the spatial spectra of intensity. The numerical requirements for a simulation of given accuracy are determined for realizations of the field. The numerical requirements for accurate estimation of higher moments of the field are less stringent.

Coles, William A.↗
NASA NTRS1981

Numerical simulation of soil brightness temperatures at wavelength of 21 cm

A simulation model is applied to reproduce some observed brightness temperatures at a wavelength of 21 cm. The simulated results calculated with two different soil textures are compared directly with observations measured over fields in Arizona and South Dakota. It is found that good agreement is possible by properly adjusting the surface roughness parameter. Correlation analysis and linear regression of the brightness temperatures versus soil moistures are also carried out.

Mo, T.↗
DOE OSTI2022

Implement and Test 3D Mortar Contact in BISON

We leverage the extension of the generation of mortar segment meshes to three dimensions in MOOSE’s framework to extend thermomechanical modeling capabilities to problems with three dimensions. A modular approach to gap heat transfer physics using the mortar finite element method was created and documented, mechanical contact was extended to three dimensions—including frictional behavior, performance and ease of use were improved, and steps towards scalability of solid mechanics problems involving contact were taken. Many of these new developments are demonstrated in the simulation of 3D light-water reactor (LWR) problems, where the thermomechanical interface problem is solved using the mortar finite element method. Usage of the mortar framework has improved convergence in 2D problems and has enabled employing friction in 3D problems, of which we show results of a short, local stack of 3D pellets. Consequently, the benefits of mortar in terms of solution convergence and quality are extended to three dimensions. Section 2 discusses fundamental developments that enabled the simulation of practical mortar problems in three dimensions and other general improvements, including the reduction of the derivative container size, the modification of dual basis computations when edge dropping (lack of secondary element projection) takes place, the improvement of conditioning when employing the VCP in-edge dropping conditions, and code usability and quality improvements. These latter code enhancements include the migration of tests using “old” mortar contact constraints to using dual mortar with a semi-smooth Newton solution strategy and the reuse of lower dimensional domains for straightforwardly setting up a mortar thermomechanical LWR problem, i.e. the MOOSE action is employed for mechanical contact and the thermal LWR action is employed to capture the gas conductance, contact, and radiation components of gap heat transfer physics. Independently of the mortar LWR thermal action, we developed a modular approach to gap heat transfer that resides in MOOSE and can be leveraged, e.g., in metallic fuel problems. This approach, whose code design based on MOOSE’s user objects to model specific physics was proposed by the maintenance activity, is detailed in Section 3. Based on the dual mortar finite element method, the frictional contact constraints were extended to three dimensions. A block sheared in two directions in and out of contact with a rigid plane is employed in Section 4 to show the way the approach handles changes in frictional states (e.g. stick to slip) within a competitive number of Newton iterations. Equations and numerical results on the use of the VCP with Cartesian Lagrange multipliers, whose combination enables their direct condensation, are described in Section 5.3. Two-dimensional and three-dimensional BISON LWR simulations are discussed in Section 6. Particularly, a stack of five eccentric pellets with a surface defect is simulated and the effect of pellet-cladding friction is assessed. Finally, conclusions are outlined in Section 7.

42 ENGINEERING↗
DOE OSTI2020

Dynamic Response of a Semiactive Suspension System with Hysteretic Nonlinear Energy Sink Based on Random Excitation by means of Computer Simulation

This paper aims to investigate the property and behavior of the hysteretic nonlinear energy sink (HNES) coupled to a half vehicle system which is a nine-degree-of-freedom, nonlinear, and semiactive suspension system in order to improve the ride comfort and increase the stability in shock mitigation by using the computer simulation method. The HNES model is a semiactive suspension device, which comprises the famous Bouc–Wen (B-W) model employed to describe the force produced by both the purely hysteretic spring and linear elastic spring of potentially negative stiffness connected in parallel, for the half vehicle system. Nine nonlinear motion equations of the half vehicle system are derived in terms of the seven displacements and the two dimensionless hysteretic variables, which are integrated numerically by employing the direct time integration method for studying both the variables of vertical displacements, velocities, accelerations, chassis pitch angle, and the ride comfort and driver safety, respectively, based on the bump and random road inputs of the pseudoexcitation method as excitation signal. Simulation results show that, compared with the HNES model and the magnetorheological (MR) model coupled to the half vehicle system, the ride comfort and stability have been evidently improved. A successful validation process has been performed, which indicated that both the ride comfort and driver safety properties of the HNES model coupled to half vehicle significantly improved.

Chen, Hui↗
NASA NTRS1986

Simulating AOTV heating environment in an arc jet

The application of arc-jet facilities to the simulation of the combined radiative and convective heating environments calculated for AOTV return missions is presented. The AOTV mission requirement and design configuration determine the magnitude of the radiative and convective heating rates. This paper considers methods for obtaining radiation from the high temperature plasma in arc-jet flows. A computer code that predicts arc-heater performance was used to determine the effect of gas mixtures on arc-heater performance, with the intent of increasing the gas temperature and resultant radiative emission. Numerical radiation solutions were compared with direct measurements of the radiative intensity from the gas cap. Convective heating measurements were made as well. The arc-jet was operated at AOTV enthalpy and pressure levels for these measurements to assess the simulation capabilities for AOTV material development/testing. A parametric study is presented that defines the capability of the Interaction Heating Facility arc jet at Ames to simulate both the radiative and the convective heating for current AOTV configurations.

Balboni, J.↗
NASA NTRS1999

Highlights of Recent Developments in the International GPS Service and Perspectives for Future Directions

In December 1998, the IGS Governing Board officially changed the name of this IAG service from 'International GPS Service for Geodynamics to simply the 'International GPS Service'. This change of name reflects the fact that today the IGS supports numerous scientific projects outside the traditional geodetic and geodynamic disciplines. A number of IGS projects and working groups have been established, each concentrating on a particular science application, such as the ionosphere, atmosphere, reference frame, precise time transfer, etc. These activities are enabled and simulated by the IGS and directly contribute to the continuing development of the service. The IGS is currently poised to respond to evolving user requirements that focus on timeliness and reliability of data and products, particularly in support of a slate of Low Earth Orbiter missions over the next decade. Perspectives on the future of the IGS will be developed based on current directions as well as anticipated external influences, such as GPS satellite modernization, GLONASS, availability of global communications, and plans for the European GALILEO (Global Navigation Satellite System - GNSS). We will address development of user friendly interfaces and IGS product tutorials.

Neilan, R.↗
NASA NTRS1988

A diagnostic technique used to obtain cross range radiation centers from antenna patterns

A diagnostic technique to obtain cross range radiation centers based on antenna radiation patterns is presented. This method is similar to the synthetic aperture processing of scattered fields in the radar application. Coherent processing of the radiated fields is used to determine the various radiation centers associated with the far-zone pattern of an antenna for a given radiation direction. This technique can be used to identify an unexpected radiation center that creates an undesired effect in a pattern; on the other hand, it can improve a numerical simulation of the pattern by identifying other significant mechanisms. Cross range results for two 8' reflector antennas are presented to illustrate as well as validate that technique.

Lee, T. H.↗