Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Development of a C-ELS Specimen-Based Numerical Benchmark for Mode II Delamination and Assessment of Two VCCT-Based Propagation Strategies

A finite element (FE) benchmark example inspired by the calibrated end-loaded split (C-ELS) specimen is developed and used to assess the performance of delamination propagation capabilities based on linear elastic fracture mechanics (LEFM). The C-ELS specimen has the advantage of a longer region of stable delamination propagation compared to the existing mode II benchmark case. The new benchmark example may therefore provide a better assessment tool by enabling more stable crack growth in regions further away from the boundary conditions or load application. First, a benchmark result is created manually using two-dimensional finite element models of the C-ELS specimen with different delamination lengths. Second, the performance of the virtual crack closure technique (VCCT) delamination propagation capabilities in the Abaqus/Standard®1 FE code and the recently developed Progressive Release eXplicit-VCCT (PRX-VCCT) method are assessed by comparing the results to the benchmark case. Two examples with different starter delamination lengths are studied. A shorter starter length is chosen to create a scenario with unstable delamination propagation. A longer delamination encourages stable delamination propagation. Detailed results from three-dimensional analyses with aligned and misaligned meshes and two levels of mesh refinement are provided. In general, good agreement can be achieved between the results obtained from the quasi-static propagation analysis and the benchmark analysis. Numerical artifacts including anomalous unreleased nodes in the crack wake and zig-zag crack fronts occur for propagation analyses using Abaqus/Standard VCCT. In comparison, continuous, smooth, delamination fronts are observed for PRX-VCCT. The use of the benchmark case to assess different VCCT-based propagation strategies illustrates the value of establishing benchmark cases.

Composite Materials↗
NASA NTRS1986

Sediment-transport experiments in zero-gravity

One of the important parameters in the analysis of sediment entrainment and transport is gravitational attraction. The availability of a laboratory in Earth orbit would afford an opportunity to conduct experiments in zero and variable gravity environments. Elimination of gravitational attraction as a factor in such experiments would enable other critical parameters (such as particle cohesion and aerodynamic forces) to be evaluated much more accurately. A Carousel Wind Tunnel (CWT) is proposed for use in conducting experiments concerning sediment particle entrainment and transport in a space station. In order to test the concept of this wind tunnel design a one third scale model CWT was constructed and calibrated. Experiments were conducted in the prototype to determine the feasibility of studying various aeolian processes and the results were compared with various numerical analysis. Several types of experiments appear to be feasible utilizing the proposed apparatus.

Iversen, J. D.↗
NASA NTRS1987

Sediment-transport experiments in zero-gravity

One of the important parameters in the analysis of sediment entrainment and transport is gravitational attraction. The availability of a laboratory in earth orbit would afford an opportunity to conduct experiments in zero and variable gravity environments. Elimination of gravitational attraction as a factor in such experiments would enable other critical parameters (such as particle cohesion and aerodynamic forces) to be evaluated much more accurately. A Carousel Wind Tunnel (CWT) is proposed for use in conducting experiments concerning sediment particle entrainment and transport in a space station. In order to test the concept of this wind tunnel design a one third scale model CWT was constructed and calibrated. Experiments were conducted in the prototype to determine the feasibility of studying various aeolian processes and the results were compared with various numerical analysis. Several types of experiments appear to be feasible utilizing the proposed apparatus.

Iversen, James D.↗
NASA NTRS1976

On a numerical solution of the supersonic panel flutter eigenproblem

An automated digital computer procedure is presented in this paper which enables efficient solution of the eigenvalue problem associated with the supersonic panel flutter phenomena. The step-by-step incremental solution procedure is based on an inverse iteration technique which effectively utilizes solution results from the previous step in determining such results during the current solution step. Also, the computations are limited to the determination of a few specific roots only, which are expected to contain the flutter mode, and this is achieved at each step without having to compute any other root. The structural discretization achieved by the finite-element method yields highly banded stiffness, mass, and aerodynamic matrices; the aerodynamic matrix evaluated by the linearized piston theory is real but unsymmetric in nature. The solution algorithm presented in this paper fully exploits the banded form of the associated matrices, and the resulting computer program written in FORTRAN V for the JPL UNIVAC 1108 computer proves to be most efficient and economical when compared to existing procedures of such analysis. Numerical results are presented for a two-dimensional panel flutter problem.

Gupta, K. K.↗
DOE OSTI2025

Time correlations from steady-state expectation values

Recovering properties of correlation functions is typically challenging. On one hand, experimentally, it requires measurements with a temporal resolution finer than the system's dynamics. On the other hand, analytical or numerical analysis requires solving the system evolution. Here, we use recent results of quantum metrology with continuous measurements to derive general lower bounds on the relaxation and second-order correlation times that are both easy to calculate and measure. These bounds are based solely on steady-state expectation values and their derivatives with respect to a control parameter, and can be readily extended to the autocorrelation of arbitrary observables. We validate our method on two examples of critical quantum systems: a critical driven-dissipative resonator, where the bound matches analytical results for the dynamics, and the infinite-range Ising model, where only the steady state is solvable and thus the bound provides information beyond the reach of existing analytical approaches. Our results can be applied to experimentally characterize ultrafast systems, and to theoretically analyze many-body models with dynamics that are analytically or numerically hard.

Górecki, Wojciech [INFN, Pavia]↗
DOE OSTI2025

Time correlations from steady-state expectation values

Recovering properties of correlation functions is typically challenging. On one hand, experimentally, it requires measurements with a temporal resolution finer than the system's dynamics. On the other hand, analytical or numerical analysis requires solving the system evolution. Here, we use recent results of quantum metrology with continuous measurements to derive general lower bounds on the relaxation and second-order correlation times that are both easy to calculate and measure. These bounds are based solely on steady-state expectation values and their derivatives with respect to a control parameter, and can be readily extended to the autocorrelation of arbitrary observables. We validate our method on two examples of critical quantum systems: a critical driven-dissipative resonator, where the bound matches analytical results for the dynamics, and the infinite-range Ising model, where only the steady state is solvable and thus the bound provides information beyond the reach of existing analytical approaches. Our results can be applied to experimentally characterize ultrafast systems, and to theoretically analyze many-body models with dynamics that are analytically or numerically hard.

Górecki, Wojciech [INFN, Pavia]↗
NASA NTRS1999

Computation of Steady and Unsteady Laminar Flames: Theory

In this paper we describe the numerical analysis underlying our efforts to develop an accurate and reliable code for simulating flame propagation using complex physical and chemical models. We discuss our spatial and temporal discretization schemes, which in our current implementations range in order from two to six. In space we use staggered meshes to define discrete divergence and gradient operators, allowing us to approximate complex diffusion operators while maintaining ellipticity. Our temporal discretization is based on the use of preconditioning to produce a highly efficient linearly implicit method with good stability properties. High order for time accurate simulations is obtained through the use of extrapolation or deferred correction procedures. We also discuss our techniques for computing stationary flames. The primary issue here is the automatic generation of initial approximations for the application of Newton's method. We use a novel time-stepping procedure, which allows the dynamic updating of the flame speed and forces the flame front towards a specified location. Numerical experiments are presented, primarily for the stationary flame problem. These illustrate the reliability of our techniques, and the dependence of the results on various code parameters.

Hagstrom, Thomas↗
NASA NTRS2017

MPCV Exercise Operational Volume Analysis

In order to minimize the loss of bone and muscle mass during spaceflight, the Multi-purpose Crew Vehicle (MPCV) will include an exercise device and enough free space within the cabin for astronauts to use the device effectively. The NASA Digital Astronaut Project (DAP) has been tasked with using computational modeling to aid in determining whether or not the available operational volume is sufficient for in-flight exercise.Motion capture data was acquired using a 12-camera Smart DX system (BTS Bioengineering, Brooklyn, NY), while exercisers performed 9 resistive exercises without volume restrictions in a 1g environment. Data were collected from two male subjects, one being in the 99th percentile of height and the other in the 50th percentile of height, using between 25 and 60 motion capture markers. Motion capture data was also recorded as a third subject, also near the 50th percentile in height, performed aerobic rowing during a parabolic flight. A motion capture system and algorithms developed previously and presented at last years HRP-IWS were utilized to collect and process the data from the parabolic flight [1]. These motions were applied to a scaled version of a biomechanical model within the biomechanical modeling software OpenSim [2], and the volume sweeps of the motions were visually assessed against an imported CAD model of the operational volume. Further numerical analysis was performed using Matlab (Mathworks, Natick, MA) and the OpenSim API. This analysis determined the location of every marker in space over the duration of the exercise motion, and the distance of each marker to the nearest surface of the volume. Containment of the exercise motions within the operational volume was determined on a per-exercise and per-subject basis. The orientation of the exerciser and the angle of the footplate were two important factors upon which containment was dependent. Regions where the exercise motion exceeds the bounds of the operational volume have been identified by determining which markers from the motion capture exceed the operational volume and by how much. A credibility assessment of this analysis was performed in accordance with NASA-STD-7009 prior to delivery to the MPCV program.

mathematical models↗
NASA NTRS2008

A Modeling Technique and Representation of Failure in the Analysis of Triaxial Braided Carbon Fiber Composites

Quasi-static tests have been performed on triaxially braided carbon fiber composite materials with large unit cell sizes. The effects of different fibers and matrix materials on the failure mode were investigated. Simulations of the tests have been performed using the transient dynamic finite element code, LS-DYNA. However, the wide range of failure modes observed for the triaxial braided carbon fiber composites during tests could not be simulated using composite material models currently available within LS-DYNA. A macroscopic approach has been developed that provides better simulation of the material response in these materials. This approach uses full-field optical measurement techniques to measure local failures during quasi-static testing. Information from these experiments is then used along with the current material models available in LS-DYNA to simulate the influence of the braided architecture on the failure process. This method uses two-dimensional shell elements with integration points through the thickness of the elements to represent the different layers of braid along with a new analytical method for the import of material stiffness and failure data directly. The present method is being used to examine the effect of material properties on the failure process. The experimental approaches used to obtain the required data will be described, and preliminary results of the numerical analysis will be presented.

Littell, Justin D.↗
NASA NTRS2008

A Modeling Technique and Representation of Failure in the Analysis of Triaxial Braided Carbon Fiber Composites

Quasi-static tests have been performed on triaxially braided carbon fiber composite materials with large unit cell sizes. The effects of different fibers and matrix materials on the failure mode were investigated. Simulations of the tests have been performed using the transient dynamic finite element code, LS-DYNA. However, the wide range of failure modes observed for the triaxial braided carbon fiber composites during tests could not be simulated using composite material models currently available within LS-DYNA. A macroscopic approach has been developed that provides better simulation of the material response in these materials. This approach uses full-field optical measurement techniques to measure local failures during quasi-static testing. Information from these experiments is then used along with the current material models available in LS-DYNA to simulate the influence of the braided architecture on the failure process. This method uses two-dimensional shell elements with integration points through the thickness of the elements to represent the different layers of braid along with a new analytical method for the import of material stiffness and failure data directly. The present method is being used to examine the effect of material properties on the failure process. The experimental approaches used to obtain the required data will be described, and preliminary results of the numerical analysis will be presented.

Littell, Justin D.↗
NASA NTRS2008

Investigation, Analysis, and Testing of Self-contained Oxygen Generators

Self Contained Oxygen Generators (SCOGs) have widespread use in providing emergency breathing oxygen in a variety of environments including mines, submarines, spacecraft, and aircraft. These devices have definite advantages over storing of gaseous or liquid oxygen. The oxygen is not generated until a chemical briquette containing a chlorate or perchlorate oxidizer and a solid metallic fuel such as iron is ignited starting a thermal decomposition process allowing gaseous oxygen to be produced. These devices are typically very safe to store, easy to operate, and have primarily only a thermal hazard to the operator that can be controlled by barriers or furnaces. Tens of thousands of these devices are operated worldwide every year without major incident. This report examines the rare case of a SCOG whose behavior was both abnormal and lethal. This particular type of SCOG reviewed is nearly identical to a flight qualified version of SCOG slated for use on manned space vehicles. This Investigative Report is a compilation of a NASA effort in conjunction with other interested parties including military and aerospace to understand the causes of the particular SCOG accident and what preventative measures can be taken to ensure this incident is not repeated. This report details the incident and examines the root causes of the observed SCOG behavior from forensic evidence. A summary of chemical and numerical analysis is provided as a background to physical testing of identical SCOG devices. The results and findings of both small scale and full scale testing are documented on a test-by-test basis along with observations and summaries. Finally, conclusions are presented on the findings of this investigation, analysis, and testing along with suggestions on preventative measures for any entity interested in the safe use of these devices.

Keddy, Christopher P.↗
DOE OSTI2025

A detailed study of pre-heating effects in electron beam melting powder bed fusion process

Metal-based additive manufacturing processes, such as powder bed fusion with electron beam (PBF-EB) process, also referred to as electron beam melting (EBM), can produce high-density parts with minimal residual stresses due to the uniform and coherent preheating of the powder bed. However, understanding and controlling the multiple stages of preheating is required to enable the production of high-quality, consistent parts of various materials. This work presents a large-scale, multi-layer, three-dimensional numerical analysis focused on studying the preheating stages for predicting thermal history during the PBF-EB process. The model follows a continuous multi-stage cyclic process, that incorporates all the main stages of the PBF-EB process for 316 L stainless steel. This includes the gradual deposition of a new powder layer, the first and second preheating levels of the powder bed, and the energy deposition during melting (excluding the actual melt-pool behavior simulation). The model employs an adaptive time-scaling approach that automatically adjusts the energy deposition for each solution time-increment. This allows for localized changes in time-resolution over an otherwise computationally expensive multi-layer procedure. The material property variations are also taken into account, with an emphasis on the subtle irreversible changes in powder effective thermal conductivity after the two requisite preheating stages of the powder bed. This effect is studied using simplified conductivity models from the literature for partially sintered powder, validated by a dedicated experiment and numerical simulation. The large-scale model is then used to estimate the actual temperatures during first and second preheating levels for 316 L steel, which is not yet fully supported commercially for PBF-EB. Model predictions are corroborated by experiments, using and analyzing IR images, taken at the completion of each layer by the machine’s built-in infrared camera. The current model also incorporates a qualitative assessment for the effects of conductivity change during pre-heating, as well as evaluates the applicability of the time-scaling approach.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Digital Rock Physics and 3D Printing for Fractured Porous Media

Imaging techniques for the analysis of porous structures have revolutionized our ability to quantitatively characterize geomaterials. Digital representations of rock from CT images and physics modeling based on these pore structures provide the opportunity to further advance our quantitative understanding of fluid flow, geomechanics, and geochemistry, and the emergence of coupled behaviors. Additive manufacturing, commonly known as 3D printing, has revolutionized production of custom parts with complex internal geometries. For the geosciences, recent advances in 3D printing technology may be co-opted to print reproducible porous structures derived from CT-imaging of actual rocks for experimental testing. The use of 3D printed microstructure allows us to surmount typical problems associated with sample-to-sample heterogeneity that plague rock physics testing and to test material response independent from pore-structure variability. Together, imaging, digital rocks and 3D printing potentially enables a new workflow for understanding coupled geophysical processes in a real, but well-defined setting circumventing typical issues associated with reproducibility, enabling full characterization and thus connection of physical phenomena to structure. Here we report on our research exploring the possibilities that these technologies can bring to geosciences for coupled multiscale experimental and numerical analysis using 3D printed fractured rock specimens.

42 ENGINEERING↗
DOE OSTI2021

Exsolution Synthesis of Nanocomposite Perovskites with Tunable Electrical and Magnetic Properties

Nanostructured functional oxides play an important role in enabling clean energy technologies and novel memory and processor devices. Using thin-film La 0.6 Sr 0.4 FeO 3 (LSF) as a model system, the novel utility of exsolution in fabricating self-assembled metal oxide nanocomposites with tunable functionalities is shown. Exsolution triggers the formation of metallic iron (Fe 0 ) nanoparticles, Ruddlesden–Popper domains, and nm-scale percolated Fe-deficient channels in LSF. Combining multimodal characterization with numerical modeling, the chemical, magnetic, and electrical properties of the exsolution-synthesized nanocomposite at different stages of Fe 0 exsolution as well as during redox cycling are assessed. After exsolution, the electronic conductivity of the nanocomposite LSF increased by more than two orders of magnitude. Based on numerical analysis representing all the constituents, it is expected that this increase in conductivity originates mainly from the Fe-deficient percolating channels formed during exsolution. Moreover, the exsolved nanocomposite is redox-active even at moderate temperatures. Such redox capabilities can enable dynamic control of the nanocomposite functionality by tailoring the oxygen non-stoichiometry. This concept is demonstrated with a continuous modulation of magnetization between 0 and 110 emu cm -3 . These findings point out that exsolution may serve as a platform for scalable fabrication of complex metal oxide nanocomposites for electrochemical and electronic applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Hybrid Plasmonic Coupled‐Mode Metasurface Design for Shot‐Noise Limited Ppb‐Level Hydrogen Detection

Plasmonic hydrogen sensors have enabled hydrogen detection below parts-per-million (ppm) range by boosting the sensitivity using localized surface plasmonic resonant (LSPR) structures. However, the intrinsic optical losses of Palladium (Pd), the primary plasmonic metal used for hydrogen detection, result in a low quality (Q) factor LSPR, which fundamentally hinders further improvement. In this work, a hybrid plasmonic metasurface is proposed that couples Pd-based LSPR structure with an Au film supporting surface plasmon polariton mode (Au-SPP). The coupled near-perfect absorber resonance yields a spectrally narrow, high Q response that retains strong sensitivity to hydrogen while improving resonance localization. Numerical analysis shows that, under shot-noise-limited conditions, the limit of detection (LoD) can be improved by over threefold compared to the state-of-the-art designs. Furthermore, this hybrid plasmonic coupled-mode metasurface thus presents a promising pathway to achieve parts-per-billion-level (ppb-level) hydrogen detection with enhanced spectral precision and robustness.

coupled modes↗
DOE OSTI2020

Cascade decays of heavy Higgs bosons through vectorlike quarks in two Higgs doublet models

We study cascade decays of heavy neutral Higgs bosons through vectorlike quarks. We focus on scenarios where decay modes into pairs of vectorlike quarks are not kinematically open which extends the sensitivity of the LHC to larger masses. Assuming only mixing with the third family of standard model quarks the new decay modes of heavy Higgs bosons are: H → t 4 t → Wbt, Ztt, htt and H → b 4 b → Wtb, Zbb, hbb, where t 4 (b 4 ) is the new up-type (down-type) quark mass eigenstate. In the numerical analysis we assume the CP even Higgs boson in the two Higgs doublet model type-II but the signatures are relevant for many other scenarios. We identify the region of the parameter space where these decay modes are significant or can even dominate, and thus they provide the best opportunities for the simultaneous discovery of a new Higgs boson and vectorlike quarks. We further explore the reach of the High Luminosity LHC for two representative decay modes,t 4 → Zt → ℓℓt and b 4 → Zb → ℓℓb, and found that cross sections at a 0.1 fb level can be probed with simple cut based analyses. We also find that the rates for Higgs cascade decays can be much larger than the rates for a single production of vectorlike quarks. Furthermore, the reach for vectorlike quarks in Higgs cascade decays and pair production extends to comparable masses.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

One-loop charge-breaking minima in the two-Higgs doublet model

We analyze the vacuum structure of the one-loop effective potential in the two Higgs doublet model. We find that electroweak-breaking vacuua can coexist with charge breaking ones, contradicting a theorem valid at tree-level. We perform a numerical analysis of the model and supply explicit parameter values for which charge-breaking vacuua can be the global minimum of the theory, and deeper than charge-preserving ones.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

A flippon related singlet at the LHC II

Here, we consider the 750 GeV diphoton resonance at the 13 TeV LHC in the ℱ-SU(5) model with a Standard Model (SM) singlet field which couples to TeV-scale vector-like particles, dubbed flippons . This singlet field assumes the role of the 750 GeV resonance, with production via gluon fusion and subsequent decay to a diphoton via the vector-like particle loops. We present a numerical analysis showing that the observed 8 TeV and 13 TeV diphoton production cross-sections can be generated in the model space with realistic electric charges and Yukawa couplings for light vector-like masses. We further discuss the experimental viability of light vector-like masses in a General No-Scale ℱ-SU(5) model, offering a few benchmark scenarios in this consistent GUT that can satisfy all experimental constraints imposed by the LHC and other essential experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗