Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Boundary condition”

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 577 records · Page 32

Considerations for testing full-scale wind turbine nacelles with hardware-in-the-loop

In this paper, pole placement control techniques are utilized to compensate for variations in the torsional dynamics resulting from the different boundary conditions experienced by a wind turbine nacelle when mounted on a test bench. Simulation case studies are used to understand the baseline torsional characteristics of the nacelle coupled to the test bench versus that coupled to the rotor. It was found that the frequency shifts by up to 200% and 400% for the first and second torsional modes, respectively. The feasibility of utilizing pole placement technique to match the torsional characteristics of the turbine on the test bench is investigated using PI and PID controllers. Here, the performance of the tuned controller is then verified under two test scenarios: Low Voltage Ride Through (LVRT), and a highly dynamic turbulent wind input. It is demonstrated that tuning the test bench speed controller can effectively shift the poles of the closed-loop system to match the desired first and second modes of the full turbine and replicate transient field events on the test bench.

17 WIND ENERGY↗

Comparing Transformation Pathways Across Major Economies

This paper aims at providing an overview of the transformation pathways for the seven largest greenhouse-gas (GHG) emitting regions and countries in the world, in line with their contributions to limiting global average temperature increase to well below 2oC as compared to pre-industrial levels. We introduce the methodology for developing these pathways by initially discussing the iterative process by which global integrated assessment models (IAMs) initiated discussions by suggesting boundary conditions in the form of carbon budgets for different regions and countries. The carbon budgets subsequently agreed upon were then used in eleven different regional and national energy-economy models and IAMs for deriving low-carbon pathways in line with a well-below-2oC world. We present a comparative assessment of the resulting regional and national pathways, and of the challenges and opportunities associated with them.

Schaeffer, Roberto↗

High fidelity simulations of contaminant dispersion in an urban environment with comparison to magnetic resonance imaging measurements

The dispersion of a contaminant in an urban environment has the potential to impact a large population of people. In this work, a complex urban canopy flow based on the Oklahoma City downtown business district circa 2003 is studied using Magnetic Resonance Imaging (MRI) and high-fidelity Large Eddy Simulations (LES). MRI is a novel experimental technique that can provide high-resolution measurements in four dimensions (three spatial and temporal) for lab scale models. The experiments and simulations use the same geometry and boundary conditions providing a one-to-one comparison of the two methods. Results are presented on the time-averaged velocity and concentration fields, the temporal dynamics of the concentration plumes for a transient release, and a novel Cloud Identification Algorithm that can separate plumes produced by periodic contaminant releases used for ensemble averaging over many releases. The MRI and LES datasets both include millions of measurement voxels and the comparisons highlight the complex 3D nature of the flow including strong vertical velocities in spanwise street canyons and flow acceleration in streamwise street canyons. The concentration fields are qualitatively similar albeit the LES shows larger dispersion. A quantitative analysis with performance measures compares the datasets pointwise and demonstrates that the two 3D datasets are similar with respect to many measures including a fractional bias of 0.02 (ideal=0.0), correlation coefficient of 0.87 (ideal = 1.0), and the fraction points within a factor of 2 is 0.98 (ideal = 1.0). Plume analysis compares the arrival and residence time of contaminant and is found to vary significantly with location within the urban environment with arrival times between 0 and 1.25 and differences within the contaminant cloud less than 10% at most locations.

54 ENVIRONMENTAL SCIENCES↗

Field–potential finite-difference time-domain (FiPo FDTD) technique for computational electromagnetics

Modeling light–matter interactions at the nanoscale requires accurate handling of coupled quantum and electromagnetic systems. This coupling requires information about the electric scalar potential Φ and the magnetic vector potential A, which are not typically calculated in standard computational electromagnetics implementations. To that end, we have developed a field–potential finite-difference time-domain (FiPo FDTD) algorithm, which solves a set of first-order equations for Φ and A alongside equations for the electric and magnetic fields E and H. The FiPo Basic code is essentially conventional FDTD, but with an added module that calculates the potentials. The FiPo Hybrid code self-consistently calculates both fields and potentials and is particularly suitable for coupling with quantum electronic transport solvers because it can be sourced by the potentials themselves. To terminate the domain and mimic infinite space, we have derived and implemented a convolutional perfectly matched layer (CPML) absorbing boundary condition for FiPo FDTD whose performance is on par with state-of-the-art CPMLs for standard FDTD. We present FiPo simulation results on several example systems.

Avazpour, L. [University of Wisconsin-Madison, WI ↗

Nonoverlapping Localized Exponential Time Differencing Methods for Diffusion Problems

In this paper, we propose nonoverlapping localized exponential time differencing (ETD) methods for diffusion problems. The model time-dependent diffusion equation is first reformulated on subdomains based on the nonoverlapping domain decomposition, in which Neumann boundary conditions are imposed on the interfaces for the subdomain problems and Dirichlet type conditions are enforced to form a space-time interface problem. After spatial discretization by standard central finite differences and temporal integration with the first or second order ETD methods, the fully discrete interface problem is obtained. Such an interface problem is then solved iteratively either at each time step or over the whole time interval: the former involves the solution of stationary problems in each subdomain at each iteration while the latter involves the solution of time-dependent subdomain problems at each iteration. For both approaches, we prove that localized ETD solutions conserve mass exactly and converge in time to the exact space semidiscrete solution. Numerical experiments in two dimensions are also presented to illustrate the performance of the proposed methods.

Mass conservation↗

Nonlinear Impact-Echo Test for Quantitative Evaluation of ASR Damage in Concrete

Nonlinear acoustic methods demonstrate high sensitivity for concrete damage evaluation. Among these methods, nonlinear resonance acoustic methods have been widely used in laboratory tests on small scale test samples. In this study, a nonlinear impact-echo (IE) method for concrete damage evaluation was introduced, in which the IE frequency shift was measured at different impact force levels. The nonlinear IE test and the nonlinear impact resonance acoustic spectroscopy (NIRAS) test share similarities in the experimental setup and analysis method. However, the nonlinear IE test excites a local thickness resonance mode of a member instead of the global resonance vibration. Therefore, the analyzed mode is unaffected by member boundary conditions and is applicable to large concrete structural members. To identify the fundamental IE frequencies, multiple impacts were applied along the members. Once the fundamental IE frequencies were determined, the effectiveness of the nonlinear IE method was evaluated by testing seven concrete beam specimens with varying levels of alkali-silica reaction (ASR) damage. In conclusion, the nonlinear IE test demonstrates high sensitivity to concrete damage and allows for quantitative assessment of the damage state of concrete without a baseline measurement.

36 MATERIALS SCIENCE↗

Mutual Diffusion Coefficients and Refractive Index Increments of K 2 SO 4 (aq) at 298.15 K from Rayleigh Interferometry

Here, isothermal mutual diffusion coefficients (interdiffusion coefficients) were measured for K 2 SO 4 (aq) at 298.15 ± 0.005 K, at numerous concentrations ranging from dilute solutions to near saturation (0.59648 mol∙dm -3 ; 0.61349 mol∙kg -1 ) under free diffusion boundary conditions, using high precision Rayleigh interferometry. Under the experimental conditions these diffusion coefficients are on the volume-fixed reference frame D v . Two series of experiments were performed, the first using the traditional experimental approach with a mercury lamp source and with the interference patterns being recorded on glass photographic plates, and the second with a He-Ne laser light source and computerized data acquisition using a photodiode array. The results from both series of experiments are in excellent agreement, and generally yield diffusion coefficients precise to at least 0.002 x 10 -9 m 2 ∙s 1 (0.15% to 0.19%) and in most cases to 0.001 x 10 -9 m 2 ∙s -1 (0.07% to 0.1%). These experiments also yield accurate values of the refractive index differences for the solution pairs used in the diffusion experiments. The new diffusion coefficients are compared to two sets of published values of diffusion coefficients for K 2 SO 4 (aq) which are somewhat discrepant from each other. This study extends and complements our earlier work on the diffusion coefficients of the most common brine salts: NaCl(aq), KCl(aq), MgCl 2 (aq), CaCl 2 (aq), Na 2 SO 4 (aq), ] MgSO 4 (aq), NaHCO 3 (aq), and KHCO 3 (aq) at 298.15 K.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Uncoupling Electrokinetic Flow Solutions

The continuum-scale electrokinetic porous-media flow and excess charge redistribution equations are uncoupled using eigenvalue decomposition. The uncoupling results in a pair of independent diffusion equations for “intermediate” potentials subject to modified material properties and boundary conditions. The fluid pressure and electrostatic potential are then found by recombining the solutions to the two intermediate uncoupled problems in a matrix-vector multiplication. Expressions for the material properties or source terms in the intermediate uncoupled problem may require extended precision or careful rewriting to avoid numerical cancellation, but the solutions themselves can typically be computed in double precision. The approach works with analytical or gridded numerical solutions and is illustrated through two examples. The solution for flow to a pumping well is manipulated to predict streaming potential and electroosmosis, and a periodic one-dimensional analytical solution is derived and used to predict electroosmosis and streaming potential in a laboratory flow cell subjected to low frequency alternating current and pressure excitation. The examples illustrate the utility of the eigenvalue decoupling approach, repurposing existing analytical solutions or numerical models and leveraging solutions that are simpler to derive for coupled physics.

58 GEOSCIENCES↗

Multiscale Modeling Meets Machine Learning: What Can We Learn?

Machine learning is increasingly recognized as a promising technology in the biological, biomedical, and behavioral sciences. There can be no argument that this technique is incredibly successful in image recognition with immediate applications in diagnostics including electrophysiology, radiology, or pathology, where we have access to massive amounts of annotated data. However, machine learning often performs poorly in prognosis, especially when dealing with sparse data. This is a field where classical physics-based simulation seems to remain irreplaceable. In this review, we identify areas in the biomedical sciences where machine learning and multiscale modeling can mutually benefit from one another: Machine learning can integrate physics-based knowledge in the form of governing equations, boundary conditions, or constraints to manage ill-posted problems and robustly handle sparse and noisy data; multiscale modeling can integrate machine learn- ing to create surrogate models, identify system dynamics and parameters, analyze sensitivities, and quantify uncertainty to bridge the scales and understand the emergence of function. With a view towards applications in the life sciences, we discuss the state of the art of combining machine learning and multiscale modeling, identify applications and opportunities, raise open questions, and address potential challenges and limitations. We anticipate that it will stimulate discussion within the community of computational mechanics and reach out to other disciplines including mathematics, statistics, computer science, artificial intelligence, biomedicine, systems biology, and precision medicine to join forces towards creating robust and efficient models for biological systems.

machine learning, multiscale modeling, physics-bas↗

On the Impact of Mechanics on Electrochemistry of Lithium-Ion Battery Anodes

Abstract Models exploring electrochemistry-mechanics coupling in liquid electrolyte lithium-ion battery anodes have traditionally incorporated stress impact on thermodynamics, bulk diffusive transport, and fracture, while stress-kinetics coupling is more explored in the context of all solid-state batteries. Here, we showcase the existence of strong link between active particle surface pressure and reaction kinetics affecting performance even in liquid electrolyte systems. Traction-free and immobile particle surface mechanical boundary conditions are used to delineate the varying pressure magnitudes in graphite host during cycling. Both tensile and compressive stresses are generated in traction-free case, while a fixed surface subjects the entire particle to a compression state. Pressure magnitudes are nearly two to three orders of magnitude higher for the latter resulting in significant depression of open circuit potential and improvement of exchange current densities compared to stress-free state. The results demonstrate the need for incorporating stress-kinetics linkage in models and provide a rationale for putting battery electrodes under compression to improve kinetics.

25 ENERGY STORAGE↗

An innovative heterogeneous transfer learning framework to enhance the scalability of deep reinforcement learning controllers in buildings with integrated energy systems

Deep Reinforcement Learning (DRL)-based control shows enhanced performance in the management of integrated energy systems when compared with Rule-Based Controllers (RBCs), but it still lacks scalability and generalisation due to the necessity of using tailored models for the training process. Transfer Learning (TL) is a potential solution to address this limitation. However, existing TL applications in building control have been mostly tested among buildings with similar features, not addressing the need to scale up advanced control in real-world scenarios with diverse energy systems. This paper assesses the performance of an online heterogeneous TL strategy, comparing it with RBC and offline and online DRL controllers in a simulation setup using EnergyPlus and Python. The study tests the transfer in both transductive and inductive settings of a DRL policy designed to manage a chiller coupled with a Thermal Energy Storage (TES). The control policy is pre-trained on a source building and transferred to various target buildings characterised by an integrated energy system including photovoltaic and battery energy storage systems, different building envelope features, occupancy schedule and boundary conditions (e.g., weather and price signal). The TL approach incorporates model slicing, imitation learning and fine-tuning to handle diverse state spaces and reward functions between source and target buildings. Results show that the proposed methodology leads to a reduction of 10% in electricity cost and between 10% and 40% in the mean value of the daily average temperature violation rate compared to RBC and online DRL controllers. Moreover, online TL maximises self-sufficiency and self-consumption by 9% and 11% with respect to RBC. Conversely, online TL achieves worse performance compared to offline DRL in either transductive or inductive settings. However, offline Deep Reinforcement Learning (DRL) agents should be trained at least for 15 episodes to reach the same level of performance as the online TL. Therefore, the proposed online TL methodology is effective, completely model-free and it can be directly implemented in real buildings with satisfying performance.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

On the Order Parameter of the Continuous Phase Transition in the Classical and Quantum Mechanical Limits

The mean field theory is revisited in the classical and quantum mechanical limits. Taking into account the boundary conditions at the phase transition and the third law of the thermodynamics the physical properties of the ordered and disordered phases were reported. The equation for the order parameter predicts the occurrence of a saturation of Ψ 2 = 1 near Θ S , the temperature below the quantum mechanical ground state is reached. The theoretical predictions are also compared with high resolution thermal expansion data of SrTiO 3 monocrystalline samples and other some previous results. An excellent agreement has been found suggesting a universal behavior of the theoretical model to describe continuous structural phase transitions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Balancing Impedance and Controllability in Response Reconstruction with $\mathrm{TS-IMMAT}$

In a prior work the authors proposed a variant on the Impedance Matched Multi-Axis Test (IMMAT) in which a fixture is defined, called the Transmission Simulator (TS), and the desired environment is matched at a set of sensors on the TS. If the motion of the TS is matched then the response of the rest of the component will also match, provided that the attached component has the same dynamics as it did when the environment was measured. Hence, one would like the TS to be flexible so that it reproduces the boundary conditions that the component of interest experiences during flight, but the more flexible the TS, the more shakers might be needed to control its response. This work presents a derivation that gives expressions for these two potential error sources in TS-IMMAT. Then, various case studies are presented, both on simulated and real hardware, to understand the importance of each error term in practical testing. The theory explains the phenomena that were observed when using measurements from a component that flew on a sounding rocket. The environmental response was measured and then various fixtures were attached, each comprising more of the next assembly, or the hardware to which the component was attached in flight. MIMO testing was repeated with each fixture and the results were compared to seek to understand the role of the impedance match in this type of testing. Here the results show that the number of modes that are active in the transmission simulator is also very important, and so the best solution balances these two considerations. An improved method of simulating the MIMO test is then proposed, so simulations can be used to predict what fixture, or transmission simulator, will give the best results in a TS-IMMAT test.

42 ENGINEERING↗

A Simplified Method for Predicting Shaker Voltage in IMMATs

Impedance Matched Multi-Axis Tests (IMMATs) can replicate in-service vibration induced stress more accurately than single axis shaker table tests as they can better match a part’s operational boundary conditions and excite it in multiple degrees of freedom simultaneously. Here, the shakers used in IMMATs are less powerful than shaker tables, so shaker force limits can be exceeded during tests if they are not placed adequately for the desired environment. The ability to predict shaker voltage and force before performing a test is, therefore, helpful in selecting shaker locations so that their limits are not exceeded. In this study, electrodynamic shakers were modeled as discrete electromechanical systems, and the shaker parameters were chosen to match experimentally obtained acceleration/voltage frequency response functions (FRFs). These models were coupled to a finite element model of the device under test (DUT) via dynamic substructuring, and the substructured model was demonstrated to accurately predict shaker voltage as well as the error in reproducing the environment at multiple accelerometer locations. A simple method called the FRF Multiplication method, in which the FRF of the substructured system is approximated as the product of two separate FRFs of the shaker and DUT respectively, was proposed and applied to the same system, yielding similar voltage and error predictions to those obtained using substructuring. Simple case studies were presented to explore the applicability of the proposed method, and it was demonstrated to have similar accuracy to the substructuring method in a range of cases. Additionally, we showed that while it was not possible to derive a unique model of the shakers from acceleration/voltage FRFs alone, the models that could be obtained were sufficient to predict test error almost perfectly and shaker voltage with less than 40 percent error.

42 ENGINEERING↗

Optimization of an Impedance-Matched Test Fixture with the Modal Projection Error

Across many industries and engineering disciplines, components and systems are designed and deployed into their operational environment of intended use. It is the desire of the design agency to be able to predict whether their component or system will function in its shock and vibration environments or if it will fail due to mechanical stresses. One method to determine if the component will survive the shock and vibration environments is to expose the component to the operational environment in a laboratory. One difficulty in executing a representative laboratory test is that the component may not have the same boundary condition in the laboratory as in the operational configuration. This paper examines the use of parameterized optimization to design the test fixture in order to better match the operational configuration. Several frequency and modal-based objective functions are examined for the optimization. The study shows that the Modal Projection Error objective function performs the best of the functions studied. The efficacy of the Modal Projection Error is demonstrated with respect to dynamic test fixture design on several analytical and experimental exemplars.

dynamic↗

Impedance Based Fixture Neutralization (FINE) Method to Replicate Field Environments in Laboratory Settings

Test article vibration test responses will vary from one test laboratory to another due to differences in shaker and fixture dynamic characteristics. This is also seen between the field and laboratory due to the differing dynamic characteristics between field assembly loading and boundary conditions compared to the laboratory configuration. This work introduces a technique called Impedance Fixture Neutralization which customizes input forces to cause consistent responses for the Device Under Test across different vibration testing conditions. The customized force neutralizes the dynamic variations between configurations for the Device Under Test. The device responses can be replicated in several situations: when using a different fixture with the same attachment points, when the test article mounting location on the fixture changes, when the force location changes, or any combination of these situations. Impedance Fixture Neutralization uses the uncoupled dynamic characteristics of the test article and excitation fixtures to customize the input thereby causing the same test article responses between two mounting configurations (either field to laboratory or two different laboratory configurations). The application of the technique is shown using an analytical model of a two-beam system and analytically using experimental FRFs from plate and frame component hardware characterization tests. In both cases a device under test is attached to two test fixtures wherein the dynamic differences in the configurations are neutralized.

Laboratory↗

Numerical and experimental investigation of the geometry dependent layer-wise evolution of temperature during laser powder bed fusion of Ti–6Al–4V

Abstract Laser powder bed fusion (L-PBF) is currently the additive manufacturing process with the widest industrial use for metal parts. Yet some hurdles persist on the way to a widespread industrial serial production, with reproducibility of the process and the resulting part properties being a major concern. As the geometry changes, so do the local boundary conditions for heat dissipation. Consequently, the use of global, geometry-independent processing parameters, which are today’s state of the art, may result in varying part properties or even defects. This paper presents a numerical simulation as a method to predict the geometry-dependent temperature evolution during the build. For demonstration, an overhang structure with varying angles towards the build platform was manufactured using Ti–6Al–4V. A calibrated infrared camera was integrated into a commercial L-PBF system to measure the temperature evolution over time for a total build height of 10 mm, and the results are used for validation of the simulation. It is shown that the simulation is capable of predicting the temperature between layers. The deviations between simulation and measurement remain in single digit range for smaller overhang structures (90°, 60° and 45°). For large overhang structures (30°), the simulation tends to over-predict the temperatures up to 15 °C. Experiments with varying process parameters showed the feasibility of energy reduction as compensation of the heat accumulation produced by overhang structures.

Li, Gefei (ORCID:0000000174837952)↗

Exascale granular microstructure reconstruction in 3D volumes of arbitrary geometries with generative learning

Reconstructing 3D granular microstructures within volumes of arbitrary geometries from limited 2D image data is crucial for predicting the material properties, as well as performances of structural components accounting for material microstructural effects. We present a novel generative learning framework that enables exascale reconstruction of granular microstructures within complex 3D geometric volumes. Building upon existing transfer learning techniques using pre-trained convolutional neural networks (CNN), we introduce several key innovations to overcome the difficulties inherent in arbitrary geometries. Our framework incorporates periodic boundary conditions using circular padding techniques, ensuring continuity and representativeness of the reconstructed microstructures. We also introduce a novel seamless transition reconstruction (STR) method that creates statistically equivalent transition zones to integrate multiple pre-existing 3D microstructure volumes. Based on STR, we propose a cost-effective strategy for reconstructing microstructures within complex geometric volumes, minimizing computational waste. Validation through numerical experiments using kinetic Monte Carlo simulations demonstrates accurate reproduction of grain statistics, including grain size distributions and morphology. A case study involving the reconstruction of a 4-blade propeller microstructure illustrates the method’s capability to efficiently handle complex geometries. In conclusion, the proposed framework significantly reduces computational demands while maintaining high reconstruction quality, paving the way for scalable microstructure reconstruction in materials design and analysis.

36 MATERIALS SCIENCE↗