SEARCH · Engineering Papers

Results for “numerical 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 901 records · Page 50

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Hydro-chemo-mechanical phase field formulation for corrosion induced cracking in reinforced concrete

Steel corrosion induced cracking in reinforced concrete structures is often caused by volumetric straining at the steel-concrete interface. This volumetric strain is mainly controlled by the rate of generation of corrosion product, which depends on the diffusion of Cl{sup −} and dissolved O{sub 2} as well as the chemical reaction at the steel-concrete interface. While many numerical modeling techniques have been proposed to tackle this complex phenomenon, most have been limited to using staggered solution schemes, where the diffusion, chemical reaction, mechanical strain, and crack fields are solved for separately. Unfortunately, such approaches introduce numerical errors that can lead to unrealistic predictions. Therefore, a coupled solution, is needed to remove the formulation induced errors and understand this complex phenomenon. To provide a unifying solution framework for corrosion induced cracking in reinforced concrete, a phase field formulation is presented and implemented with FEM algorithms. The governing equations are derived from the system's free energy and Fick's law of diffusion. The proposed formulation was first validated in 1D by comparing the modeling results with the semi-analytical solution approximated via Laurent series to avoid singularities in the complex domain. 2D plane strain solutions to the embedded rebar undergoing corrosion induced cracking with pre-cracks are also presented and compared with the approximated solutions from linear elastic fracture mechanics theory (LEFM). The modeling of crack propagation enabled by the proposed formulation is also presented along with parametric studies to reveal the roles of the fracture toughness and permeability of the concrete on the corrosion induced cracking. A case study was also conducted on H-pile steel with concrete jackets to demonstrate the feasibility of modeling corrosion induced fracture in reinforced concrete structures with complex steel-concrete interfaces.

36 MATERIALS SCIENCE↗

NASA NTRS1997

A Three-Dimensional Linearized Unsteady Euler Analysis for Turbomachinery Blade Rows

A three-dimensional, linearized, Euler analysis is being developed to provide an efficient unsteady aerodynamic analysis that can be used to predict the aeroelastic and aeroacoustic responses of axial-flow turbo-machinery blading.The field equations and boundary conditions needed to describe nonlinear and linearized inviscid unsteady flows through a blade row operating within a cylindrical annular duct are presented. A numerical model for linearized inviscid unsteady flows, which couples a near-field, implicit, wave-split, finite volume analysis to a far-field eigenanalysis, is also described. The linearized aerodynamic and numerical models have been implemented into a three-dimensional linearized unsteady flow code, called LINFLUX. This code has been applied to selected, benchmark, unsteady, subsonic flows to establish its accuracy and to demonstrate its current capabilities. The unsteady flows considered, have been chosen to allow convenient comparisons between the LINFLUX results and those of well-known, two-dimensional, unsteady flow codes. Detailed numerical results for a helical fan and a three-dimensional version of the 10th Standard Cascade indicate that important progress has been made towards the development of a reliable and useful, three-dimensional, prediction capability that can be used in aeroelastic and aeroacoustic design studies.

Montgomery, Matthew D.↗

NASA NTRS2000

Constituent-Based Life Models Being Developed for SiC/SiC Composites

For the successful utilization of ceramic matrix composites (CMC) as hot-section components in advanced aeropropulsion engines, the CMC constituents will need to be tailored and optimized to meet all the critical property demands of each component. Under the High-Speed Research (HSR) and Advanced High-Temperature Engine Materials Technology (HITEMP) Programs, the NASA Glenn Research Center at Lewis Field initiated research to develop mechanistic models for key CMC thermostructural properties. These models would describe the effects of different constituent factors (composition, geometry, and volume fraction) and of potential application conditions (stress, time, temperature, and environment) on these properties. Particular focus was placed on both analytical and numerical modeling of state-of-the-art SiC/SiC composites where the primary load-bearing constituents are stoichiometric SiC fibers in a complex multiphase SiC matrix produced by chemical vapor infiltration and melt infiltration. Recent studies have resulted in computer-based numerical models for the elastic modulus, thermal expansion, and thermal conductivity properties of the SiC/SiC system. Additional studies have generated analytical and empirical models for the time dependence of composite rupture strength at temperatures above 2200 F (1200 C), where CMC's have an important thermostructural advantage over current nickel-based superalloys. These life models utilize thermal activation theory and fiber stress-rupture results measured at Glenn to generate Larson-Miller (L M) plots of fiber rupture strength versus q, a single time- and temperature-dependent parameter. Assuming a worse case in which the SiC matrix is cracked, rupture is then controlled by the time-dependent fracture characteristics of the fiber bundles bridging the matrix cracks. With this as the controlling mechanism, one can then use simple composite theory and the fiber L M plots to predict CMC rupture strength versus the q parameter. The dashed line shows the predicted rupture strength of a SiC/SiC composite that is reinforced by a state-of-the-art stoichiometric SiC fiber. For the q parameter, time is in hours and temperature in degrees kelvin. To generate these predictions, a two-dimensional 0 /90 composite with approximately 16 percent fiber in the applied stress direction and an air test environment were assumed. As such, it is possible to compare the model predictions against limited stress-rupture data for this CMC as shown by the data points. The good agreement confirms the rupture model at least for the selected CMC and test conditions. Thus for this particular SiC/SiC composite, one can estimate a 1000-hr rupture strength of approximately 12 ksi at 2400 F (T = 1588 K and q = 39 700 K). At lower CMC application stresses, the SiC matrix is typically uncracked, so both the fiber and matrix constituents share the composite load. In this case, CMC rupture is controlled by the constituent with the longest rupture time based on the creep rate of the composite. Measured Monkman-Grant plots of rupture time versus creep rate for the two SiC constituents have been used to develop CMC life models for this important application condition. NASA and DOD are currently using this information to establish application and material goals for more advanced CMC's that can be used at even higher temperatures.

DiCarlo, James A.↗

DOE OSTI2022

Three‐Dimensional Permeability Inversion Using Convolutional Neural Networks and Positron Emission Tomography

Abstract Quantification of heterogeneous multiscale permeability in geologic porous media is key for understanding and predicting flow and transport processes in the subsurface. Recent utilization of in situ imaging, specifically positron emission tomography (PET), enables the measurement of three‐dimensional (3‐D) time‐lapse radiotracer solute transport in geologic media. However, accurate and computationally efficient characterization of the permeability distribution that controls the solute transport process remains challenging. Leveraging the relationship between local permeability variation and solute advection behavior, an encoder‐decoder based convolutional neural network (CNN) is implemented as a permeability inversion scheme using a single PET scan of a radiotracer pulse injection experiment as input. The CNN can accurately capture the 3‐D spatial correlation between the permeability and the radiotracer solute arrival time difference maps in geologic cores. We first test the inversion accuracy using synthetic test datasets and then test the accuracy on a suite of experimental PET imaging datasets acquired on four different geologic cores. The network‐predicted permeability maps from the geologic cores are used to parameterize forward numerical models that are directly compared with the experimental PET imaging data. The results indicate that a single trained network can generate robust 3‐D permeability inversion maps in seconds. Numerical models parameterized with these permeability maps closely capture the experimentally observed solute arrival time behavior. This work provides an unprecedented approach for efficiently characterizing multiscale permeability heterogeneity in complex geologic samples.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1995

Large HAWT wake measurement and analysis

From the theoretical fluid dynamics point of view, the wake region of a large horizontal-axis wind turbine has been defined and described, and numerical models of wake behavior have been developed. Wind tunnel studies of single turbine wakes and turbine array wakes have been used to verify the theory and further refine the numerical models. However, the effects of scaling, rotor solidity, and topography on wake behavior are questions that remain unanswered. In the wind tunnel studies, turbines were represented by anything from scaled models to tea strainers or wire mesh disks whose solidity was equivalent to that of a typical wind turbine. The scale factor compensation for the difference in Reynolds number between the scale model and an actual turbine is complex, and not typically accounted for. Though it is wise to study the simpler case of wakes in flat topography, which can be easily duplicated in the wind tunnel, current indications are that wind turbine farm development is actually occurring in somewhat more complex terrain. Empirical wake studies using large horizontal-axis wind turbines have not been thoroughly composited, and, therefore, the results have not been applied to the well-developed theory of wake structure. The measurement programs have made use of both in situ sensor systems, such as instrumented towers, and remote sensors, such as kites and tethered, balloonborne anemometers. We present a concise overview of the work that has been performed, including our own, which is based on the philosophy that the MOD-2 turbines are probably their own best detector of both the momentum deficit and the induced turbulence effect downwind. Only the momentum deficit aspects of the wake/machine interactions have been addressed. Both turbine power output deficits and wind energy deficits as measured by the onsite meteorological towers have been analyzed from a composite data set. The analysis has also evidenced certain topographic influences on the operation of spatially diverse wind turbines.

Miller, A. H.↗

NASA NTRS1981

Analytical and experimental investigation of the propagation and attenuation of sound in extended reaction lined ducts

Results are presented of an analytical and experimental study of the attenuation and propagation of harmonically excited sound waves in an extended reaction lined cylindrical duct. The duct geometry considered consisted of an annular outer region of bulk material surrounding an inner cylinder of air. The coupled wave equations governing the motion of the sound in both the inner and annular regions were solved numerically. The numerically predicted attenuation and propagation constants were in excellent agreement with measured values using Kevlar as the liner material for plane-wave mode (0,0) excitation over the frequency from 100 to 7,000 Hz. Although the numerical model was verified using Kevlar, it can be used with any fibrous constructed bulk liner. The results of this study demonstrate that a good start has been made on the numerical modelling of the acoustic performance of extended reaction liners.

Hersh, A. S.↗

DOE OSTI2022

A Review of Computational Capabilities and Requirements in High-Resolution Simulation of Nonferrous Pyrometallurgical Furnaces

We report that tremendous progress has been made over the last several decades in improving the pyrometallurgical processing routes of nonferrous metals, such as aluminum, copper, and lead. Advances in the numerical modeling of pyrometallurgical processes has aided in these improvements by providing a better understanding of the complex transport phenomena occurring in modern furnaces. However, there is a need for a comprehensive discussion of the numerical modeling of a primary and a secondary nonferrous pyrometallurgical furnaces. This review provides such a discussion by surveying recent attempts at capturing the physico-chemical phenomena occurring within these furnaces, including gas-phase combustion, melting/smelting, and multiphase heat- and mass-transfer with the molten phases. This work then identifies a complete set of approaches for simulating these types of furnaces and provides recommendations for applying high-resolution numerical tools towards full-furnace simulation. By identifying gaps in the current state of the art, this review offers suggestions for future developments in commercial codebases, outlining the path to fulfilling the promise of CFD-enabled pyrometallurgy.

36 MATERIALS SCIENCE↗

NASA NTRS2008

Measurement and Prediction of the Thermomechanical Response of Shape Memory Alloy Hybrid Composite Beams

An experimental and numerical investigation into the static and dynamic responses of shape memory alloy hybrid composite (SMAHC) beams is performed to provide quantitative validation of a recently commercialized numerical analysis/design tool for SMAHC structures. The SMAHC beam specimens consist of a composite matrix with embedded pre-strained SMA actuators, which act against the mechanical boundaries of the structure when thermally activated to adaptively stiffen the structure. Numerical results are produced from the numerical model as implemented into the commercial finite element code ABAQUS. A rigorous experimental investigation is undertaken to acquire high fidelity measurements including infrared thermography and projection moire interferometry for full-field temperature and displacement measurements, respectively. High fidelity numerical results are also obtained from the numerical model and include measured parameters, such as geometric imperfection and thermal load. Excellent agreement is achieved between the predicted and measured results of the static and dynamic thermomechanical response, thereby providing quantitative validation of the numerical tool.

Davis, Brian↗

NASA NTRS2017

The Latest on the Venus Thermospheric General Circulation Model: Capabilities and Simulations

Venus has a complex and dynamic upper atmosphere. This has been observed many times by ground-based, orbiters, probes, and fly-by missions going to other planets. Two over-arching questions are generally asked when examining the Venus upper atmosphere: (1) what creates the complex structure in the atmosphere, and (2) what drives the varying dynamics. A great way to interpret and connect observations to address these questions utilizes numerical modeling; and in the case of the middle and upper atmosphere (above the cloud tops), a 3D hydrodynamic numerical model called the Venus Thermospheric General Circulation Model (VTGCM) can be used. The VTGCM can produce climatological averages of key features in comparison to observations (i.e. nightside temperature, O2 IR nightglow emission). More recently, the VTGCM has been expanded to include new chemical constituents and airglow emissions, as well as new parameterizations to address waves and their impact on the varying global circulation and corresponding airglow distributions.

chemistry↗

DOE OSTI2023

Electrochemical Response of Alkaline Batteries Subject to Quasi-Static and Dynamic Loading

We report electrochemical changes were successfully monitored through the coupling of external loading, electrochemical impedance spectroscopy (EIS), and distribution of relaxation times (DRT) analysis. The development of the testing methods allows for detailed observation of changes due to mechanical loading and distinguishes responses between different cell geometries. Comparison of the force vs displacement, voltage vs displacement, and distribution of relaxation times plots to the different cell geometries display failure modes specific to geometry as well as the impact of different loading profiles. Numerical modeling confirmed the movement within individual cells and predicted locations with the most deformation based on the external loading condition applied. The results gathered from combining mechanical loading, electrochemical response, and numerical modeling, yield a viable approach to establishing an improved understanding of the effects of mechanical loading on the electrochemical response of multiple battery geometries and the methods herein may be extensible to additional battery chemistries.

25 ENERGY STORAGE↗

DOE OSTI2020

Extreme Scale Infrasound Inversion and Prediction for Weather Characterization and Acute Event Detection

Accurate and timely weather predictions are critical to many aspects of society with a profound impact on our economy, general well-being, and national security. In particular, our ability to forecast severe weather systems is necessary to avoid injuries and fatalities, but also important to minimize infrastructure damage and maximize mitigation strategies. The weather community has developed a range of sophisticated numerical models that are executed at various spatial and temporal scales in an attempt to issue global, regional, and local forecasts in pseudo real time. The accuracy however depends on the time period of the forecast, the nonlinearities of the dynamics, and the target spatial resolution. Significant uncertainties plague these predictions including errors in initial conditions, material properties, data, and model approximations. To address these shortcomings, a continuous data collection occurs at an effort level that is even larger than the modeling process. It has been demonstrated that the accuracy of the predictions depends on the quality of the data and is independent to a certain extent on the sophistication of the numerical models. Data assimilation has become one of the more critical steps in the overall weather prediction business and consequently substantial improvements in the quality of the data would have transformational benefits. This paper describes the use of infrasound inversion technology, enabled through exascale computing, that could potentially achieve orders of magnitude improvement in data quality and therefore transform weather predictions with significant impact on many aspects of our society.

58 GEOSCIENCES↗

NASA NTRS1997

Finite Element Analysis of Magnetic Damping Effects on G-Jitter Induced Fluid Flow

This paper reports some interim results on numerical modeling and analyses of magnetic damping of g-jitter driven fluid flow in microgravity. A finite element model is developed to represent the fluid flow, thermal and solute transport phenomena in a 2-D cavity under g-jitter conditions with and without an applied magnetic field. The numerical model is checked by comparing with analytical solutions obtained for a simple parallel plate channel flow driven by g-jitter in a transverse magnetic field. The model is then applied to study the effect of steady state g-jitter induced oscillation and on the solute redistribution in the liquid that bears direct relevance to the Bridgman-Stockbarger single crystal growth processes. A selection of computed results is presented and the results indicate that an applied magnetic field can effectively damp the velocity caused by g-jitter and help to reduce the time variation of solute redistribution.

Pan, Bo↗

NASA NTRS2010

Waterhammer Modeling for the Ares I Upper Stage Reaction Control System Cold Flow Development Test Article

The Upper Stage Reaction Control System provides three-axis attitude control for the Ares I launch vehicle during active Upper Stage flight. The system design must accommodate rapid thruster firing to maintain the proper launch trajectory and thus allow for the possibility to pulse multiple thrusters simultaneously. Rapid thruster valve closure creates an increase in static pressure, known as waterhammer, which propagates throughout the propellant system at pressures exceeding nominal design values. A series of development tests conducted in the fall of 2009 at Marshall Space Flight Center were performed using a water-flow test article to better understand fluid performance characteristics of the Upper Stage Reaction Control System. A subset of the tests examined waterhammer along with the subsequent pressure and frequency response in the flight-representative system and provided data to anchor numerical models. This thesis presents a comparison of waterhammer test results with numerical model and analytical results. An overview of the flight system, test article, modeling and analysis are also provided.

Williams, Jonathan H.↗

NASA NTRS1973

Derivation of a three dimensional numerical water quality model for estuary and continental shelf application

A derivation is given for a three dimensional mass transport equation which is appropriate for numerical modeling of estuary and continental shelf water quality variations for both the time dependent and steady state cases. A finite difference approximation to the derived equation is presented and a solution scheme for the resulting equations outlined. Preliminary results are obtained using the model for the extremely simple problems which have analytical solutions. The numerical model, as presented, will provide a scheme to study water quality problems in coastal waters for both steady state and time dependent cases.

Spaulding, M.↗

DOE OSTI2024

Quantifying the Influence of Size, Shape, and Density of Microplastics on Their Transport Modes: A Modeling Approach

Microplastics (MPs) pose significant risks to marine ecosystems and human health, necessitating accurate predictions of their distributions in aquatic environments for effective risk mitigation. However, understanding MP transport dynamics is challenging because of the inadequate representation of MP characteristics such as size, shape, and density in numerical models. Further, the accuracy of the MP vertical profiles in existing models has not been thoroughly validated. Thus, we developed an MP transport model within the Finite Volume Community Ocean Model framework (FVCOM-MP) by integrating MP characteristics. We validated FVCOM-MP against experimental and analytical data, focusing on various MP transport modes and transitions. FVCOM-MP successfully replicates MP profiles in different transport modes, including the bedload, surface-load, suspended-load, and mixed-load modes. Additionally, we introduce phase diagrams for classifying MP transport modes based on particle characteristics, enhancing our understanding of MP dynamics in aquatic systems. The transport modes for a number of real-world MP particles, including fishing line, plastic bag/bottle fragments, synthetic fibers, tire wear particles, polyvinyl chloride and expanded polystyrene pellets, were analyzed with these phase diagrams.

Microplastic transport, Settling velocity, Rising ↗

DOE OSTI2020

Spatial Mapping of Riverbed Grain-Size Distribution Using Machine Learning

Recent alluvial sediments in riverbeds play a significant role in controlling hydrologic exchange flows (HEFs) in river systems. The alluvial layer is usually associated with strong heterogeneity in physical properties (e.g., permeability and hydraulic conductivity), which affects local HEFs and therefore biogeochemical processes. The spatial distribution of these physical properties needs to be determined to inform the numerical models used to reveal the realistic hydro-biogeochemical behaviors. Such information can be obtained based on the intrinsic link between sediment grain-size distribution and hydraulic properties where sediment texture information is available. However, grain-size measurements are usually spatially sparse and do not have adequate coverage and resolution, particularly for a relatively large domain such as the Hanford Reach of the Columbia River. In this paper, we adopted machine learning (ML) approaches for categorizing and mapping the spatial distributions of riverbed substrate grain size and filling in missing areas of substrate data using the ML models along the reach. Such ML models for substrate size mapping were trained at 13,372 locations using measured substrate sizes along with observed and simulated attributes, including bathymetric attributes (e.g., elevation, slope, and aspect ratio) from LIDAR and bathymetric surveys, and hydrodynamic properties (e.g., water depth, velocity, shear stress, and their statistical moments). An ensemble bagging-based ML technique, Random Forest, was adopted to identify the most influential factors as predictors to develop the predictive models with over-fitting issues addressed. The models were evaluated with respect to each individual substrate size class and the lumped group, and then used to generate the final substrate size maps covering all the grid cells in the numerical modeling domain.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2000

The fluid dynamic and shear environment in the NASA/JSC rotating-wall perfused-vessel bioreactor

The rotating-wall perfused-vessel (RWPV) bioreactor, used for both microgravity and Earth-based cell science experiments, is characterized in terms of the fluid dynamic and fluid shear stress environment. A numerical model of the flow field is developed and verified with laser Doppler velocimeter measurements. The effects of changes in operating conditions, including rotation rates and fluid perfusion rates, are investigated with the numerical model. The operating conditions typically used for ground-based experiments (equal rotation of the inner and outer cylinders) leads to flow patterns with relatively poor mass distribution characteristics. Approximately 50% of the inlet-perfused fluid bypasses the bulk of the fluid volume and flows to the perfusion exit. For operating conditions typical in microgravity, small differential rotation rates between the inner and outer cylinders lead to greatly improved flow distribution patterns and very low fluid shear stress levels over a large percentage of the fluid volume. Differences in flow patterns for the different operating conditions are explored. Large differences in the hydrodynamic environments for operating conditions typical of true microgravity and ground-based "microgravity simulations" are demonstrated.

Bioreactors↗

NASA NTRS1996

Computer Modeling of Non-Isothermal Crystallization

A realistic computer model for simulating isothermal and non-isothermal phase transformations proceeding by homogeneous and heterogeneous nucleation and interface-limited growth is presented. A new treatment for particle size effects on the crystallization kinetics is developed and is incorporated into the numerical model. Time-dependent nucleation rates, size-dependent growth rates, and surface crystallization are also included. Model predictions are compared with experimental measurements of DSC/DTA peak parameters for the crystallization of lithium disilicate glass as a function of particle size, Pt doping levels, and water content. The quantitative agreement that is demonstrated indicates that the numerical model can be used to extract key kinetic data from easily obtained calorimetric data. The model can also be used to probe nucleation and growth behavior in regimes that are otherwise inaccessible. Based on a fit to data, an earlier prediction that the time-dependent nucleation rate in a DSC/DTA scan can rise above the steady-state value at a temperature higher than the peak in the steady-state rate is demonstrated.

Kelton, K. F.↗