Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Displacement”

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,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Investigation of main bearing operating conditions in a three-Point mount wind turbine drivetrain

Wear-related failures of spherical roller bearings in the main bearing position of three-point mount wind turbines have been higher than expected and can contribute to higher than anticipated operation-and-maintenance costs. In this paper, the operational conditions of such a main bearing—including measured axial displacement and velocity subject to the estimated axial loads—are described for an instrumented, commercial wind turbine. The field measurements suggest a maximum axial speed between the bearing rings that is less than 2 millimeters per second. It is estimated that the axial speed between the ring and rollers is approximately 25% of this value. When compared to a speed of rolling from approximately 200–352 millimeters per second, the measured axial sliding is therefore significantly less than 1% of the speed of rolling. Previous numerical studies of lubricant film formation in rolling contacts have shown that the effect of axial sliding starts to be noticeable only when this ratio exceeds 10%; therefore, the axial velocity represents only a small disturbance to the nominal pure rolling case, and the influence on oil film building can be neglected. A simple analytic model of the main bearing motion was also developed and demonstrated similar displacement and velocity characteristics.

17 WIND ENERGY↗
DOE OSTI2023

Towards a unified nonlocal, peridynamics framework for the coarse-graining of molecular dynamics data with fractures

Molecular dynamics (MD) has served as a powerful tool for designing materials with reduced reliance on laboratory testing. However, the use of MD directly to treat the deformation and failure of materials at the mesoscale is still largely beyond reach. In this work, we propose a learning framework to extract a peridynamics model as a mesoscale continuum surrogate from MD simulated material fracture data sets. Firstly, we develop a novel coarse-graining method, to automatically handle the material fracture and its corresponding discontinuities in the MD displacement data sets. Inspired by the weighted essentially non-oscillatory (WENO) scheme, the key idea lies at an adaptive procedure to automatically choose the locally smoothest stencil, then reconstruct the coarse-grained material displacement field as the piecewise smooth solutions containing discontinuities. Then, based on the coarse-grained MD data, a two-phase optimization-based learning approach is proposed to infer the optimal peridynamics model with damage criterion. In the first phase, we identify the optimal nonlocal kernel function from the data sets without material damage to capture the material stiffness properties. Then, in the second phase, the material damage criterion is learnt as a smoothed step function from the data with fractures. As a result, a peridynamics surrogate is obtained. As a continuum model, our peridynamics surrogate model can be employed in further prediction tasks with different grid resolutions from training, and hence allows for substantial reductions in computational cost compared with MD. We illustrate the efficacy of the proposed approach with several numerical tests for the dynamic crack propagation problem in a single-layer graphene. Our tests show that the proposed data-driven model is robust and generalizable, in the sense that it is capable of modeling the initialization and growth of fractures under discretization and loading settings that are different from the ones used during training.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Behavioral responses of fish to a current-based hydrokinetic turbine under mutliple operational conditions

There is significant international interest in developing current-based marine and hydrokinetic (MHK) technologies to capture the power of tidal energy. However, concerns have been raised regarding the ecological effects of these projects on fish, including the risk of blade collision and behavioral impacts such as the disruption of migratory behavior and food acquisition and displacement from preferred habitats. In this study, we conducted mobile hydroacoustic surveys to track fish as they approached a tidal turbine deployed in Cobscook Bay, Maine. There was a significant decline in fish numbers with decreasing distance to the turbine, beginning approximately 140 m from the turbine. Similar declines were not observed at control transects or when the turbine was not spinning. The decline in fish numbers appeared to be the result of horizontal displacement, not vertical, movements to avoid the turbine. Noise rather than visual cues or flow field disturbance seemed to be a likely explanation for the reduced number of fish near the turbine. This finding, combined with near-field blade collision studies indicating a low probability of encounter, suggests that a single turbine poses a low collision risk to pelagic fish and that a single turbine is likely to result in minimal behavioral responses by fish. However, the risk may be different with additional devices, which will become more relevant as commercial-scale MHK arrays come under consideration. Therefore, the risks associated with commercial-scale operations will ultimately have to be evaluated to fully understand the ecological impacts of MHK devices.

42 ENGINEERING↗
DOE OSTI2022

A multi-resolution approach to hydraulic fracture simulation

Abstract We present a multi-resolution approach for constructing model-based simulations of hydraulic fracturing, wherein flow through porous media is coupled with fluid-driven fracture. The approach consists of a hybrid scheme that couples a discrete crack representation in a global domain to a phase-field representation in a local subdomain near the crack tip. The multi-resolution approach addresses issues such as the computational expense of accurate hydraulic fracture simulations and the difficulties associated with reconstructing crack apertures from diffuse fracture representations. In the global domain, a coupled system of equations for displacements and pressures is considered. The crack geometry is assumed to be fixed and the displacement field is enriched with discontinuous functions. Around the crack tips in the local subdomains, phase-field sub-problems are instantiated on the fly to propagate fractures in arbitrary, mesh independent directions. The governing equations and fields in the global and local domains are approximated using a combination of finite-volume and finite element discretizations. The efficacy of the method is illustrated through various benchmark problems in hydraulic fracturing, as well as a new study of fluid-driven crack growth around a stiff inclusion.

58 GEOSCIENCES↗
DOE OSTI2022

A smeared crack modeling framework accommodating multi-directional fracture at finite strains

A generic smeared crack modeling framework predicated on the deformation gradient decomposition (DGD) approach is proposed for use in dynamic fracture problems at finite strains, accommodating failure along multiple mutually orthogonal fracture planes embedded within an independently defined bulk material model. Within this constitutive framework, the traction equilibrium conditions imposed at each failure surface are used to determine the associated crack displacements stored as internal state variables. In general, the enforcement of interfacial equilibrium entails the implicit solution of a non-linear system of equations within the constitutive update procedure. However, if inertial effects arising due to the relative motion of the fractured material are incorporated within the model, the traction equilibrium conditions are shown to give rise to corresponding dynamic equations of motion governing the time-evolution of the crack opening displacements. For dynamic problems, an explicit time-integration procedure is devised to efficiently update the material state, subject to a set of internal frictionless contact constraints to prevent material inversion. Finally, the efficacy of the proposed modeling framework is investigated through several benchmark dynamic fracture problems run within the explicit finite element code DYNA3D.

42 ENGINEERING↗
DOE OSTI2020

A Mixed Discontinuous Galerkin Method for Linear Elasticity with Strongly Imposed Symmetry

In this paper, we study a mixed discontinuous Galerkin (MDG) method to solve linear elasticity problem with arbitrary order discontinuous finite element spaces in d -dimension ( \(d=2,3\) ). This method uses polynomials of degree \(k+1\) for the stress and of degree k for the displacement ( \(k\ge 0\) ). The mixed DG scheme is proved to be well-posed under proper norms. Specifically, we prove that, for any \(k \ge 0\) , the \(H(\mathrm{div})\) -like error estimate for the stress and \(L^2\) error estimate for the displacement are optimal. We further establish the optimal \(L^2\) error estimate for the stress provided that the \({\mathcal {P}}_{k+2}-{\mathcal {P}}_{k+1}^{-1}\) Stokes pair is stable and \(k \ge d\) . Finally, we also provide numerical results of MDG showing that the orders of convergence are actually sharp.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Influence of Wetting on Viscous Fingering Via 2D Lattice Boltzmann Simulations

We present simulations of two-phase flow using the Rothman and Keller colour gradient Lattice Boltzmann method to study viscous fingering when a “red fluid” invades a porous model initially filled with a “blue” fluid with different viscosity. We conducted eleven suites of 81 numerical experiments totalling 891 simulations, where each suite had a different random realization of the porous model and spanned viscosity ratios in the range $$M\in [0.01,100]$$ M ϵ [ 0.01 , 100 ] and wetting angles in the range $$\theta _w\in [180^\circ ,0^\circ ]$$ θ w ϵ [ 180 ° , 0 ° ] to allow us to study the effect of these parameters on the fluid-displacement morphology and saturation at breakthrough (sweep). Although sweep often increased with wettability, this was not always so and the sweep phase space landscape, defined as the difference in saturation at a given wetting angle relative to saturation for the non-wetting case, had hills, ridges and valleys. At low viscosity ratios, flow at breakthrough is localized through narrow fingers that span the model. After breakthrough, the flow field continues to evolve and the saturation continues to increase albeit at a reduced rate, and eventually exceeds 90% for both non-wetting and wetting cases. The existence of a complicated sweep phase space at breakthrough, and continued post-breakthrough evolution suggests the hydrodynamics and sweep is a complicated function of wetting angle, viscosity ratio and time, which has major potential implications to Enhanced Oil Recovery by water flooding, and hence, on estimates of global oil reserves. Validation of these results via experiments is required to ensure they translate to field studies.

Engineering↗
DOE OSTI2023

Path-Integrated X-Ray Images for Multi-Surface Digital Image Correlation (PI-DIC)

X-ray imaging offers unique possibilities for Digital Image Correlation (DIC), opening the door for full-field deformation measurements of a test article in complex environments where optical DIC suffers severe biases or is impossible. While X-ray DIC has been performed in the past with standard DIC codes designed for optical images, the path-integrated nature of X-ray images places constraints on the experimental setup, predominantly that only a single surface of interest moves/deforms. These requirements are difficult to realize for many practical situations and limit the amount of information that can be garnered in a single test. Other X-ray based diagnostics such as Digital Volume Correlation (DVC) and Projection DVC (P-DVC) overcome these obstacles, but DVC is limited to quasi-static tests, and both DVC and P-DVC necessitate high-resolution computed tomography (CT) scan(s) and often require a potentially invasive pattern throughout the volume of the specimen. Here this work presents a novel approach to measure time-resolved displacements and strains on multiple surfaces from a single series of 2D, path-integrated (PI) X-ray images, called PI-DIC. The principle of optical flow or conservation of intensity—the foundation of DIC—was reframed for path-integrated images, for an exemplar setup comprised of two plates moving and deforming independently. Synthetic images were generated for rigid translations, rigid rotations, and uniform stretches, where each plate underwent a unique motion/deformation. Experimental specimens were fabricated (either an aluminum plate with tantalum features or a plastic plate with steel features) and the two specimens were independently translated. PI-DIC was successfully demonstrated with the synthetic images and validated with the experimental images. Prescribed displacements were recovered for each plate from the single set of path-integrated, deformed images. Errors were approximately 0.02 px for the synthetic images with 1.5% image noise, and 0.05 px for the experimental images. These results provide the foundation for PI-DIC to measure motion and deformation of multiple, independent surfaces with subpixel accuracy from a single series of path-integrated X-ray images.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

In situ microstructural evolution in face-centered and body-centered cubic complex concentrated solid-solution alloys under heavy ion irradiation

This study characterizes the microstructural evolution of single-phase complex concentrated solid-solution alloy (CSA) compositions under heavy ion irradiation with the goal of evaluating mechanisms for CSA radiation tolerance in advanced fission systems. Three such alloys, Cr 18 Fe 27 Mn 27 Ni 28 , Cr 15 Fe 35 Mn 15 Ni 35 , and equimolar NbTaTiV, along with reference materials (pure Ni and E90 for the Cr-FeMnNi family and pure V for NbTaTiV) were irradiated at 50 K and 773 K with 1 MeV Kr ++ ions to various levels of displacements per atom (dpa) using in-situ transmission electron microscopy. Cryogenic irradiation resulted in small defect clusters and faulted dislocation loops as large as 12 nm in face-centered cubic (FCC) CSAs. With thermal diffusion suppressed at cryogenic temperatures, defect densities were lower in all CSAs than in their less compositionally complex reference materials indicating that point defect production is reduced during the displacement cascade stage. High temperature irradiation of the two FCC CSA resulted in the formation of interstitial dislocation loops which by 2 dpa grew to an average size of 27 nm in Cr 18 Fe 27 Mn 27 Ni 28 and 10 nm in Cr 15 Fe 35 Mn 15 Ni 35 . This difference in loop growth kinetics was attributed to the difference in Mn-content due to its effect on the nucleation rate by increasing vacancy mobility or reducing the stacking-fault energy.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Residual stresses and microstructure within Allvac 718Plus laser powder bed fusion bars

The residual stresses within Allvac 718Plus bars built using LBPF with different laser powers, speeds and table displacements were measured using diffraction and mechanical methods and modelled with the baseplate attached and then removed. The residual stress profiles within the bars from the top surface down through the bulk were all quite similar, becoming less tensile due to a strain hardening mechanism. Table displacement has the greatest impact on residual stress, decreasing with increasing displacement/powder layer thickness. There was good agreement amongst the modeling and measurements. The microstructures were examined and varied slightly with energy density with higher densities having larger grains and enhanced post solidification diffusion. We found energy density had minimal impact on the residuals stresses within the parameters to produce dense material.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Insights into palladium poisoning of Cu/SSZ-13 selective catalytic reduction catalysts

Here, the impacts of Pd poisoning to the activity, selectivity, and hydrothermal stability of Cu/SSZ-13 selective catalytic reduction (SCR) catalysts are reported. Pd lowers DeNOx efficiency of Cu/SSZ-13 via two mechanisms: (1) displacing SCR active sites in the form of isolated Pd-ions, and (2) catalyzing non-selective NH 3 oxidation in the form of PdO. The first mechanism works by the diffusion of isolated Pd-ions into chabazite cages to displace ZCu II OH SCR active species, and it occurs on Cu/SSZ-13 catalysts with and without external surface CuO clusters in similar fashion. In contrast, the second mechanism works differently with and without external surface CuO clusters. For Cu/SSZ-13 catalyst without external CuO clusters (i.e., primarily isolated Cu-ions), PdO leads to significant decrease of DeNOx efficiency at reaction temperatures above ~400°C due to its non-selective NH 3 oxidation activity at high temperatures. However, this poisoning effect becomes much less impactful on Cu/SSZ-13 catalyst containing CuO clusters. This is due to the formation of CuO-PdO solid solution via interactions between CuO and PdO, which reduces non-selective NH 3 oxidation potential. Furthermore, this solid solution formation even mitigates adverse effects caused by hydrothermal aging. Hence, the poisoning effects of Pd are closely related to Cu speciation and spatial distribution of a Cu/SSZ-13 catalyst. Finally, this study suggests a strategy in eliminating Pd poisoning, that is, the introduction of an oxide phase that effectively traps PdO but does not adversely influence SCR.

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

Perspectives on novel refractory amorphous high-entropy alloys in extreme environments

Two new refractory amorphous high-entropy alloys (RAHEAs) within the W–Ta–Cr–V and W–Ta–Cr–V–Hf systems were herein synthesized using magnetron-sputtering and tested under high-temperature annealing and displacing irradiation using in situ Transmission Electron Microscopy. While the 14W-41Ta-26Cr-19V in at.% RAHEA (defined as WTaCrV RAHEA) was found to be unstable under such tests, additions of Hf in this system composing a new quinary 24W-40Ta-18Cr-5V-13Hf in at.% RAHEA (defined as WTaCrVHf RAHEA) was found to be a route to achieve stability both under annealing and irradiation. A new effect of nanoprecipitate reassembling observed to take place within the WTaCrVHf RAHEA under irradiation indicates that a duplex microstructure composed of an amorphous matrix with crystalline nanometer-sized precipitates enhances the radiation response of the system. Further, it is demonstrated that tunable chemical complexity arises as a new alloy design strategy to foster the use of novel RAHEAs within extreme environments. New perspectives for the alloy design and application of chemically-complex amorphous metallic alloys in extreme environments are presented with focus on their thermodynamic phase stability when subjected to high-temperature annealing and displacing irradiation.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Overcoming the sensitivity vs. throughput tradeoff in Coulter counters: A novel side counter design

We report microfabricated Coulter counters are attractive for point of care (POC) applications since they are label free and compact. However, these approaches inherently suffer from a trade off between sample throughput and sensitivity. The counter measures a change in impedance due to displaced fluid volume by passing cells, and thus the counter's signal increases with the fraction of the sensing volume displaced. Reducing the size of the sensing region requires reductions in volumetric throughput in the absence of increased hydraulic pressure and sensor bandwidth. The risk of mechanical clog formation, rendering the counter inoperable, increases markedly with reductions in the size of the constriction aperture. We present here a microfluidic coplanar Coulter counter device design that overcomes the problem of constriction clogging while capable of operating in microfluidic channels filled entirely with highly conductive sample. The device utilizes microfabricated planar electrodes projecting into one side of the microfluidic channel and is easily integrated with upstream electronic, hydrodynamic, or other focusing units to produce efficient counting which could allow for dramatically increased volumetric and sample throughput. The design lends itself to simple, cost effective POC applications.

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

Asymptotically compatible reproducing kernel collocation and meshfree integration for the peridynamic Navier equation

Here, we study reproducing kernel (RK) collocation method for peridynamic Navier equation. In the first part, we apply a linear RK approximation to both displacement and dilatation, and then back-substitute dilatation and solve the peridynamic Navier equation in a pure displacement form. The RK collocation scheme converges to the nonlocal limit for a fixed nonlocal interaction length and also to the local limit as nonlocal interactions vanish. The stability is shown by comparing the collocation scheme with the standard Galerkin scheme using Fourier analysis. In the second part, we apply the RK collocation to the quasi-discrete peridynamic Navier equation and show its convergence to the correct local limit when the ratio between the nonlocal length scale and the discretization parameter is fixed. The analysis is carried out on a special family of rectilinear Cartesian grids for the RK collocation method with a designated kernel with finite support. We assume the Lamé parameters satisfy λ ≥ μ to avoid extra assumptions on the nonlocal kernel. Finally, numerical experiments are conducted to validate the theoretical results.

42 ENGINEERING↗
DOE OSTI2021

Nonlinear convergence in contact mechanics: Immersed boundary finite volume

In this report we present an immersed boundary finite volume (IBM) method for simulating quasistatic contact mechanics of linearly elastic domains at small strains. In IBM, all external boundaries and internal contacts of an object are represented by embedded surfaces inside a Cartesian mesh, which need not conform to the grid lines. The contact constraints consist of the non-penetrability condition and Coulomb’s friction law, which are discretized using special interpolation stencils and enforced via penalty parameters. The resulting nonlinear system depends on displacement unknowns only. To solve it, we use the Newton method but find that it diverges frequently. To understand the divergence pattern, we analyze a simplified 2-cell problem and show that the global convergence of Newton cannot be ensured for any choice of penalty parameters. We thus propose a modified Newton solver, which guarantees convergence for the 2-cell problem and is numerically verified to converge for all the challenging simulations considered herein. While both 1 st - and 2 nd -order variants of IBM, in displacement unknowns, are proposed, the modified Newton solver applies only to the 1 st -order variant.

42 ENGINEERING↗
DOE OSTI2024

Three-dimensional continuum point cloud method for large deformation and its verification

This study presents a strong form based meshfree collocation method, which is named Continuum Point Cloud Method, to solve nonlinear field equations derived from classical mechanics for deformed bodies in three-dimensional Euclidean space. The method and its implementation are benchmarked against a nonlinear vector field using manufactured solutions. The analysis of mechanical fields firstly focuses on the study of St. Venant Kirchhoff and compressible neo-Hookean materials. Results for various initial boundary value problems are presented, including benchmark cases involving unidirectional tension and simple shear. Subsequently, the study concludes with an analysis of a displacement-controlled simulation of a compressible neo-Hookean material, specifically a bar that is pulled to 50% of its original length and rotated 90°. The pure tension case yields a 1.5% error in displacement between computed and expected values and a combined tension and torsion loading case provides further insight into material behavior under complex loading conditions. The resulting normal axial and transverse stress-strain curves are also presented. Lastly, the consistency and robustness of the proposed nonlinear numerical schemes are successfully demonstrated through various numerical experiments.

Compressible neo-Hookean materials↗
DOE OSTI2020

A direct numerical simulation of Jet A flame kernel quenching

The safe operation of aeronautical engines requires an understanding of flame ignition, propagation and extinction. In this study, direct numerical simulations are performed using a 29 species reduced chemical mechanism for jet fuel surrogate Jet A to understand the flame quenching process. Here, initially laminar spherical flames of varying sizes and equivalence ratios are subject to an identical periodic domain of decaying and isotropic high intensity turbulence with a turbulent Reynolds number of 2400. All cases become quenched, except for the larger kernel with lower Karlovitz number. An analysis of the flame structure shows broadened preheat zone, flame shortening on the product side, differential species diffusion and partial fuel pyrolysis in the fresh mixture. Two extinction mechanisms are identified arising from flame shortening and high flame stretch. Flame shortening occurs due to turbulence-chemistry interactions that resemble the flame–flame interaction in a laminar counterflow reactant-to-reactant configuration, which contorts and breaks up the ignition kernel. Flame stretch is a local effect that attenuates the heat release rate and causes the flame to retreat towards the product mixtures, similar to what has been observed for reactant-to-product laminar counterflow flames. Chemical explosive mode analysis was also performed to quantify the flame structure and local combustion mode. The diffusion–reaction balance in pinched-off flame islands favors extinction of these smaller structures, while auto-ignition modes are observed within the flame kernel after fresh mixture is engulfed and preheated in the product kernel. Statistics of the density-weighted displacement speed conditional on local combustion mode indicates strong correlation between the local extinction mode and negative displacement speed. The local balance between diffusion and reaction ultimately determines the propensity for local extinction in both laminar and turbulent flames, the extent of which has an impact on global flame propagation.

42 ENGINEERING↗
DOE OSTI2022

Non-premixed counterflow methane flames in DC/AC/NS electric fields

The response of counterflow diffusion flames to sub-breakdown DC and AC electric fields, as well as their superposition with ns pulse discharge waveforms, is studied in the plane-to-plane electrode geometry. Sub-breakdown DC and low-frequency AC electric fields cause the flame displacement toward the grounded electrode, in the direction of the applied field, indicating that the body force on the positive ions exceeds that on the electrons and negative ions. As the AC frequency increases, the flame response becomes less pronounced, due to the reduction of the electrohydrodynamic (EHD) body force impulse over the AC half-period. The electric field in the electrode gap is determined by ps Electric Field Induced Second Harmonic (E-FISH) generation, with absolute calibration using sub-breakdown ns pulses overlapped with the measured electric field waveform. The results show that the electric field distribution across the flame in the current saturation regime follows the Laplacian field. This indicates that the space charge density in the gap is too low to distort the applied DC or AC field, consistent with the kinetic modeling predictions. Combining a nanosecond pulse discharge with a sub-breakdown DC field generates a diffuse plasma across the entire gap. Time-resolved and spatially resolved measurements of the electric field in the discharge indicate the ionization wave propagation between the electrodes. The present results do not exhibit a detectable flame displacement enhancement by ns discharge pulses combined with a sub-breakdown field, observed previously. Kinetic modeling calculations show that the absence of this effect in the plane-to-plane geometry is due to the rapid plasma self-shielding. This indicates that alternative electrode geometries limiting the self-shielding would be more effective for the plasma / electric field enhanced flameholding and flame stabilization applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗