SEARCH · Engineering Papers

Results for “nonlinear phenomena”

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 721 records · Page 40

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Understanding plasma turbulence through exact coherent structures

Plasma turbulence is a key challenge in understanding transport phenomena in magnetically confined plasmas. This work presents a generalized framework to analyze plasma turbulence that utilizes periodic orbit theory. In periodic orbit theory, doubly periodic solutions (coherent structures) of the governing equation(s) serve as building blocks of the considered turbulent dynamics. To illustrate the concept and method, the particularly simple Kuramoto–Sivashinsky (referred to here as LMRT for the original authors: LaQuey, Mahajan, Rutherford, and Tang) trapped-ion mode toy model is used. By applying numerical optimization techniques to the LMRT equation, we extract coherent spacetime patterns that represent the library of allowable fundamental structures of the equation. These structures provide a framework to systematically describe turbulence as a composition of recurrent solutions, revealing an underlying order within chaotic plasma motion. Although illustrated here using the simplified LMRT model for clarity, this framework provides a general strategy that can be extended to more complex and realistic models of plasma turbulence, including gyrokinetic systems. This offers a new method for predicting and potentially controlling transport processes in fusion plasmas by providing a bridge between nonlinear dynamical systems theory and plasma physics in the form of a generalized framework with which to analyze and understand spatially extended nonlinear partial differential equations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1999

Sonoluminescence in Space: The Critical Role of Buoyancy in Stability and Emission Mechanisms

Sonoluminescence is the term used to describe the emission of light from a violently collapsing bubble. Sonoluminescence ("light from sound") is the result of extremely nonlinear pulsations of gas/vapor bubbles in liquids when subject to sufficiently high amplitude acoustic pressures. In a single collapse, a bubble's volume can be compressed more than a thousand-fold in the span of less than a microsecond. Even the simplest consideration of the thermodynamics yields pressures on the order of 10,000 ATM. and temperatures of at least 10,000 K. On the face of things, it is not surprising that light should be emitted from such an extreme process. Since 1990 (the year that Gaitan discovered light from a single bubble) there has been a tremendous amount of experimental and theoretical research in stable, single-bubble sonoluminescence. Yet there remain four fundamental mysteries associated with this phenomenon: 1) the light emission mechanism itself; 2) the mechanism for anomalous mass flux stability; 3) the disappearance of the bubble at some critical acoustic pressure; and 4) the appearance of quasiperiodic and chaotic oscillations in the flash timing. Gravity, in the context of the buoyant force, is implicated in all four of these unexplained phenomena. We are developing microgravity experiments probing the effect of gravity on single bubble sonoluminescence. By determining the stability boundaries experimentally in microgravity, and measuring not only light emission but mechanical bubble response, we will be able to directly test the unambiguous predictions of existing theories. By exploiting the microgravity environment we will gain new knowledge impossible to obtain in earth-based labs which will enable explanations for the above mysteries. We will also be in a position to make new discoveries about bubbles which emit light.

Holt, R. Glynn↗

NASA NTRS1989

Simulation studies using multibody dynamics code DART

DART is a multibody dynamics code developed by Photon Research Associates for the Air Force Astronautics Laboratory (AFAL). The code is intended primarily to simulate the dynamics of large space structures, particularly during the deployment phase of their missions. DART integrates nonlinear equations of motion numerically. The number of bodies in the system being simulated is arbitrary. The bodies' interconnection joints can have an arbitrary number of degrees of freedom between 0 and 6. Motions across the joints can be large. Provision for simulating on-board control systems is provided. Conservation of energy and momentum, when applicable, are used to evaluate DART's performance. After a brief description of DART, studies made to test the program prior to its delivery to AFAL are described. The first is a large angle reorientating of a flexible spacecraft consisting of a rigid central hub and four flexible booms. Reorientation was accomplished by a single-cycle sine wave shape torque input. In the second study, an appendage, mounted on a spacecraft, was slewed through a large angle. Four closed-loop control systems provided control of this appendage and of the spacecraft's attitude. The third study simulated the deployment of the rim of a bicycle wheel configuration large space structure. This system contained 18 bodies. An interesting and unexpected feature of the dynamics was a pulsing phenomena experienced by the stays whole playout was used to control the deployment. A short description of the current status of DART is given.

Keat, James E.↗

NASA NTRS2003

Thruster Limitation Consideration for Formation Flight Control

Physical constraints of any real system can have a drastic effect on its performance. Some of the more recognized constraints are actuator and sensor saturation and bandwidth, power consumption, sampling rate (sensor and control-loop) and computation limits. These constraints can degrade system s performance, such as settling time, overshoot, rising time, and stability margins. In order to address these issues, researchers have investigated the use of robust and nonlinear controllers that can incorporate uncertainty and constraints into a controller design. For instance, uncertainties can be addressed in the synthesis model used in such algorithms as H(sub infinity), or mu. There is a significant amount of literature addressing this type of problem. However, there is one constraint that has not often been considered; that is, actuator authority resolution. In this work, thruster resolution and controller schemes to compensate for this effect are investigated for position and attitude control of a Low Earth Orbit formation flight system. In many academic problems, actuators are assumed to have infinite resolution. In real system applications, such as formation flight systems, the system actuators will not have infinite resolution. High-precision formation flying requires the relative position and the relative attitude to be controlled on the order of millimeters and arc-seconds, respectively. Therefore, the minimum force resolution is a significant concern in this application. Without the sufficient actuator resolution, the system may be unable to attain the required pointing and position precision control. Furthermore, fuel may be wasted due to high-frequency chattering phenomena when attempting to provide a fine control with inadequate actuators. To address this issue, a Sliding Mode Controller is developed along with the boundary Layer Control to provide the best control resolution constraints. A Genetic algorithm is used to optimize the controller parameters according to the states error and fuel consumption criterion. The tradeoffs and effects of the minimum force limitation on performance are studied and compared to the case without the limitation. Furthermore, two methods are proposed to reduce chattering and improve precision.

Yunjun Xu↗

DOE OSTI2021

Driven-dissipative phase transition in a Kerr oscillator: From semiclassical PT symmetry to quantum fluctuations

We study a minimal model that has a driven-dissipative quantum phase transition, namely a Kerr non-linear oscillator subject to driving and dissipation. Using mean-field theory, exact diagonalization, and the Keldysh formalism, we analyze the critical phenomena in this system, showing which aspects can be captured by each approach and how the approaches complement each other. Then critical scaling and finite-size scaling are calculated analytically using the quantum Langevin equation. Furthermore, the physics contained in this simple model is surprisingly rich: it includes a continuous phase transition, Z 2 symmetry breaking, PT symmetry, state squeezing, and critical fluctuations. Due to its simplicity and solvability, this model can serve as a paradigm for exploration of open quantum many-body physics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Realizability-preserving discontinuous Galerkin method for spectral two-moment radiation transport in special relativity

Here we present a realizability-preserving numerical method for solving a spectral two-moment model to simulate the transport of massless, neutral particles interacting with a steady background material moving with relativistic velocities. The model is obtained as the special relativistic limit of a four-momentum-conservative general relativistic two-moment model. Using a maximum-entropy closure, we solve for the Eulerian-frame energy and momentum. The proposed numerical method is designed to preserve moment realizability, which corresponds to moments defined by a nonnegative phase-space density. The realizability-preserving method is achieved with the following key components: (i) a discontinuous Galerkin phase-space discretization with specially constructed numerical fluxes in the spatial and energy dimensions; (ii) a strong stability-preserving implicit-explicit time-integration method; (iii) a realizability-preserving conserved to primitive moment solver; (iv) a realizability-preserving implicit collision solver; and (v) a realizability-enforcing limiter. Component (iii) is necessitated by the closure procedure, which closes higher order moments nonlinearly in terms of primitive moments. The nonlinear conserved to primitive and the implicit collision solves are formulated as fixed-point problems, which are solved with custom iterative solvers designed to preserve the realizability of each iterate. With a series of numerical tests, we demonstrate the accuracy and robustness of this discontinuous-Galerkin-implicit-explicit method.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Retrieval of the molecular orientation distribution from atom-pair angular distributions

Imaging laser-induced rotational dynamics is an important and active field due to its applications in capturing reactions in the molecular frame and in molecular imaging. Experimental measurement of the molecular orientation distribution, as a function of the Euler angles, has only been demonstrated for special cases when the detectable signal is generated along the molecular symmetry axis. Here we developed the general theory that maps the probability density distribution of the molecular orientation to the atom-pair angular distributions for nonlinear molecules. With the theory, the molecular orientation distribution can be retrieved from the measured atom-pair angular distribution, which we demonstrate experimentally using ultrafast electron diffractive imaging of impulsively aligned trifluoro-iodomethane molecules. Here, the retrieved molecular orientation distribution is in good agreement with direct numerical simulations of the time-dependent Schrodinger equation using the experimental conditions. Unlike the existing retrieval methods, the retrieval method does not require solving Schrodinger equation, works for any alignment method, and is in principle applicable to asymmetric top molecules.

74 ATOMIC AND MOLECULAR PHYSICS↗

NASA NTRS1992

Combustion instability analysis for liquid propellant rocket engines

The multi-dimensional numerical model has been developed to analyze the nonlinear combustion instabilities in liquid-fueled engines. The present pressure-based approach can handle the implicit pressure-velocity coupling in a non-iterative way. The additional scalar conservation equations for the chemical species, the energy, and the turbulent transport quantities can be handled by the same predictor-corrector sequences. This method is time-accurate and it can be applicable to the all-speed, transient, multi-phase, and reacting flows. Special emphasis is given to the acoustic/vaporization interaction which may act as the crucial rate-controlling mechanism in the liquid-fueled rocket engines. The subcritical vaporization is modeled to account for the effects of variable thermophysical properties, non-unitary Lewis number in the gas-film, the Stefan flow effect, and the effect of transient liquid heating. The test cases include the one-dimenisonal fast transient non-reacting and reacting flows, and the multi-dimensional combustion instabilities encountered in the liquid-fueled rocket thrust chamber. The present numerical model successfully demonstrated the capability to simulate the fast transient spray-combusting flows in terms of the limiting-cycle amplitude phenomena, correspondence between combustion and acoustics, and the steep-fronted wave and flame propagation. The investigated parameters include the spray initial conditions, air-fuel mixture ratios, and the engine geometry. Stable and unstable operating conditions are found for the liquid-fueled combustors. Under certain conditions, the limiting cycle behavior of the combusting flowfields is obtained. The numerical results indicate that the spray vaporization processes play an important role in releasing thermal energy and driving the combustion instability.

Kim, Y. M.↗

DOE OSTI2026

Convergence Criteria for Multiphysics Simulations

The behavior of engineered systems is often influenced by multiple physical phenomena, such as mechanical deformation, heat transfer, and chemical species transport and reactions. There are often strong interactions between these phenomena, and there is increasing interest in applying coupled-physics models to improve understanding of physical behavior under complex environmental conditions. Multiple simulation frameworks that facilitate coupled-physics simulations are in widespread use, and these employ a variety of techniques to account for interactions between those physics. Many frameworks solve the physics models independently and transfer results between them. Alternatively, a single monolithic system of equations for every physics model can be formed and solved. Each of these approaches has its benefits and drawbacks, and the optimal approach varies depending on the nature of the problem. The open-source MOOSE framework was developed targeting solution of large-scale multiphysics problems. Although it provides options for all these coupling approaches, its standard approach for multiphysics solutions is to form and solve a single monolithic system of equations containing the unknowns for all physics models. MOOSE provides a streamlined approach for users to define the solution variables, the terms in the partial differential equations pertaining to each variable, and interactions between solution variables. One aspect of the monolithic solution approach that can be problematic, however, is defining appropriate convergence criteria for the nonlinear system. A standard approach is to determine convergence is to simply take a norm of the residual vector corresponding to the full vector of unknowns. However, if the residual vector contains variables for multiple physics models, the magnitudes of those variables can differ significantly, and the variables can converge at significantly different rates from each other. It is important to ensure that the variables for each of the physics are converged, and also ensure that the convergence criteria are not excessively stringent in cases when there is little change in the solution. This talk presents representative multiphysics problems to highlight these issues, and shows strategies for convergence criteria in MOOSE that are robust for multiphysics models under a variety of conditions.

97 - MATHEMATICS AND COMPUTING↗

NASA NTRS2014

An Experimental and Computational Study on Soot Formation in a Coflow Jet Flame Under Microgravity and Normal Gravity

Upon the completion of the Structure and Liftoff in Combustion Experiment (SLICE) in March 2012, a comprehensive and unique set of microgravity coflow diffusion flame data was obtained. This data covers a range of conditions from weak flames near extinction to strong, highly sooting flames, and enabled the study of gravitational effects on phenomena such as liftoff, blowout and soot formation. The microgravity experiment was carried out in the Microgravity Science Glovebox (MSG) on board the International Space Station (ISS), while the normal gravity experiment was performed at Yale utilizing a copy of the flight hardware. Computational simulations of microgravity and normal gravity flames were also carried out to facilitate understanding of the experimental observations. This paper focuses on the different sooting behaviors of CH4 coflow jet flames in microgravity and normal gravity. The unique set of data serves as an excellent test case for developing more accurate computational models.Experimentally, the flame shape and size, lift-off height, and soot temperature were determined from line-of-sight flame emission images taken with a color digital camera. Soot volume fraction was determined by performing an absolute light calibration using the incandescence from a flame-heated thermocouple. Computationally, the MC-Smooth vorticity-velocity formulation was employed to describe the chemically reacting flow, and the soot evolution was modeled by the sectional aerosol equations. The governing equations and boundary conditions were discretized on an axisymmetric computational domain by finite differences, and the resulting system of fully coupled, highly nonlinear equations was solved by a damped, modified Newtons method. The microgravity sooting flames were found to have lower soot temperatures and higher volume fraction than their normal gravity counterparts. The soot distribution tends to shift from the centerline of the flame to the wings from normal gravity to microgravity.

microgravity↗

DOE OSTI2021

Thermomagnetic properties of Bi 2 Te 3 single crystal in the temperature range from 55 K to 380 K

Magnetothermoelectric transport provides an understanding of coupled electron-hole-phonon current in topological materials and has applications in energy conversion and cooling. In this paper, we study the Nernst coefficient, the magneto-Seebeck coefficient, and the magnetoresistance of single-crystalline Bi 2 Te 3 under external magnetic field in the range of –3T to 3T and in the temperature range of 55 K to 380 K. Moreau's relation is employed to justify both the overall trend of the Nernst coefficient and the temperature at which this coefficient changes sign. We observe a nonlinear relationship between the Nernst coefficient and the applied magnetic field in the temperature range of 55 K to 255 K. An increase in both the Nernst coefficient and the magneto-Seebeck coefficient is observed as the temperature is reduced which can be attributed to the increased mobility of the carriers at lower temperatures. First-principles density functional theory calculations were carried out to physically model the experimental data including electronic and transport properties. Simulation findings agreed with the experiments and provide a theoretical insight to justify the measurements.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Bhatnagar-Gross-Krook subgrid model for neutrino quantum kinetics

Here, we present a new subgrid model for neutrino quantum kinetics, which is primarily designed to incorporate effects of collective neutrino oscillations into neutrino-radiation-hydrodynamic simulations for core-collapse supernovae and mergers of compact objects. We approximate the neutrino oscillation term in a quantum kinetic equation using the Bhatnagar-Gross-Krook (BGK) relaxation-time prescription, and the transport equation is directly applicable for classical neutrino transport schemes. The BGK model is motivated by recent theoretical indications that nonlinear phases of collective neutrino oscillations settle into quasisteady structures. We explicitly provide basic equations of the BGK subgrid model for both multiangle and moment-based neutrino transport to facilitate the implementation of the subgrid model in the existing neutrino transport schemes. We also show the capability of our BGK subgrid model by comparing it to fully quantum kinetic simulations for fast neutrino-flavor conversion. We find that the overall properties can be well reproduced in the subgrid model; the error of angular-averaged survival probability of neutrinos is within ~20%. By identifying the source of error, we also discuss perspectives to improve the accuracy of the subgrid model.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Photoinduced frustration modulation in 𝜅-type quantum spin liquid candidates

Geometric frustration is a key parameter controlling electronic and magnetic properties of quantum spin liquid systems, yet remains challenging to tune. Here, we coherently drive molecular vibrations with midinfrared pulses in two organic quantum spin liquid candidates, the insulating 𝜅−(BEDT−TTF) 2⁢ Cu 2 ⁢(CN) 3 and the metallic 𝜅−(BEDT−TTF) 4 ⁢Hg 2.89 ⁢Br 8 , and probe their electronic response through ultrafast reflectivity measurements. We observe a nonlinear coupling between local molecular vibrations and nonlocal phonons, which is expected to directly modulate the geometric frustration of their triangular lattice. Furthermore, our findings establish a promising route to dynamically control frustration in nonbipartite quantum materials.

Frustrated magnetism↗

NASA NTRS1993

Multiple steady states in atmospheric chemistry

The equations describing the distributions and concentrations of trace species are nonlinear and may thus possess more than one solution. This paper develops methods for searching for multiple physical solutions to chemical continuity equations and applies these to subsets of equations describing tropospheric chemistry. The calculations are carried out with a box model and use two basic strategies. The first strategy is a 'search' method. This involves fixing model parameters at specified values, choosing a wide range of initial guesses at a solution, and using a Newton-Raphson technique to determine if different initial points converge to different solutions. The second strategy involves a set of techniques known as homotopy methods. These do not require an initial guess, are globally convergent, and are guaranteed, in principle, to find all solutions of the continuity equations. The first method is efficient but essentially 'hit or miss' in the sense that it cannot guarantee that all solutions which may exist will be found. The second method is computationally burdensome but can, in principle, determine all the solutions of a photochemical system. Multiple solutions have been found for models that contain a basic complement of photochemical reactions involving O(x), HO(x), NO(x), and CH4. In the present calculations, transitions occur between stable branches of a multiple solution set as a control parameter is varied. These transitions are manifestations of hysteresis phenomena in the photochemical system and may be triggered by increasing the NO flux or decreasing the CH4 flux from current mean tropospheric levels.

Stewart, Richard W.↗

NASA NTRS2010

The Continual Intercomparison of Radiation Codes (CIRC) Assessing Anew the Quality of GCM Radiation Algorithms

The simulation of changes in the Earth's climate due to solar and thermal radiative processes with global climate models (GCMs) is highly complex, depending on the parameterization of a multitude of nonlinearly coupled physical processes. In contrast, the germ of global climate change, the radiative forcing from enhanced abundances of greenhouse gases, is relatively well understood. The impressive agreement between detailed radiation calculations and highly resolved spectral radiation measurements in the thermal infrared under cloudless conditions (see, for example, Fig. 1) instills confidence in our knowledge of the sources of gaseous absorption. That the agreement spans a broad range of temperature and humidity regimes using instruments mounted on surface, aircraft, and satellite platforms not only attests to our capability to accurately calculate radiative fluxes under present conditions, but also provides confidence in the spectroscopic basis for computation of fluxes under conditions that might characterize future global climate (e.g., radiative forcing). Alas, the computational costs of highly resolved spectral radiation calculations cannot be afforded presently in GCMs. Such calculations have instead been used as the foundation for approximations implemented in fast but generally less accurate algorithms performing the needed radiative transfer (RT) calculations in GCMs. Credible climate simulations by GCMs cannot be ensured without accurate solar and thermal radiative flux calculations under all types of sky conditions: pristine cloudless, aerosol-laden, and cloudy. The need for accuracy in RT calculations is not only important for greenhouse gas forcing scenarios, but is also profoundly needed for the robust simulation of many other atmospheric phenomena, such as convective processes.

Oreopoulos, Lazaros↗

DOE OSTI2024

Single crystal growth, chemical defects, magnetic and transport properties of antiferromagnetic topological insulators ( Ge 1 – δ – x Mn x ) 2 Bi 2 Te 5 ( x ≤ 0.47 , 0.11 ≤ δ ≤ 0.20 )

Magnetic topological insulators provide a platform for emergent phenomena arising from the interplay between magnetism and band topology. Here we report the single crystal growth, crystal structure, magnetic and transport properties, as well as the neutron scattering studies of topological insulator series (Ge 1–δ–x Mn x ) 2 ⁢Bi 2 ⁢Te 5 (x ≤ 0.47, 0.11 ≤ δ ≤ 0.20). Upon doping up to x = 0.47, the lattice parameter c decreases by 0.8%, while the lattice parameter a remains nearly unchanged. Significant Ge vacancies and Ge/Bi site mixing are revealed via elemental analysis as well as refinements of the neutron and x-ray diffraction data, resulting in holes dominating the charge transport. At x = 0.47, below 10.8 K, a bilayer A-type antiferromagnetic ordered state emerges, featuring an ordered moment of 3.0(3) μ B /Mn at 5 K, with the c axis as the easy axis. Magnetization data unveils a much stronger effective interlayer antiferromagnetic exchange interaction and a much smaller uniaxial anisotropy compared to MnBi 2 ⁢Te 4 . We attribute the former to the shorter nearest-neighbor Mn-Mn interlayer superexchange path and the latter to the smaller ligand-field splitting in (Ge 1–δ–x ⁢Mn x ) 2 ⁢Bi 2 ⁢Te 5 . Finally, our study demonstrates that this series of materials holds promise for the investigation of the layer Hall effect and quantum metric nonlinear Hall effect.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Multiparticle cumulant mapping for Coulomb explosion imaging: Calculations and algorithm

We present a versatile cumulant mapping algorithm for analyzing correlated particle emission, offering insights into complex electronic and nuclear dynamics. Recently, we have demonstrated the use of cumulant mapping to extract information-rich correlations between the momenta of multiple fragments produced in Coulomb explosion imaging experiments [C. Cheng et al., Phys. Rev. Lett. 130, 093001 (2023)]. We define cumulant mapping in terms of histograms, enabling fast computation of linear (additive) observables. However, applying the same algorithm to nonlinear (nonadditive) observables poses challenges, as the computation time of conventional estimators scales nonlinearly with data size. To overcome this, we develop estimators and an accompanying algorithm to enable computationally efficient estimation of the cumulant of interest. Comparisons of computation times and signal-to-noise ratios reveal the superior performance of our approach. This method is demonstrated on the (D+, D+, C+, O+) dissociation channel of CD 2 ⁢O 4+ produced in a strong-field ionization experiment. Additionally, Poisson statistics are used to simulate the two methods and provide insights into the efficiency of our algorithm. The proposed methodology unlocks efficient computation of cumulant mapping for a broader range of complex systems and observables, such as the laser pulse dependence of ionization dynamics.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2025

Toward engineering lattice structures with the material point method (MPM)

This study examines the potential of two variants of the material point method—the generalized interpolation material point (GIMP) and dual domain material point (DDMP) methods—in developing a robust computational framework for engineering lattice structures under different loading conditions. The study begins with assessing the ability of the two methods in predicting elastic buckling phenomena using column geometries with and without initial geometric imperfections. The results indicate that both methods effectively capture buckling phenomena when initial geometric imperfections are introduced. After this verification step, we create several models of tetrahedral lattice structures with varying strut diameter and orientation and subject them to quasi-static loading. We then validate the numerical results using laboratory test results. The results show that, while both methods accurately predict load–displacement curves in the pre-buckling regime, their predictive capabilities diminish in the post-buckling regime. Through visual comparison between the numerical and experimental deformed shapes, it appears that the discrepancies between model and experimental results are attributed to initial geometric imperfections in the lattices that occurred during 3D printing. We then establish a second set of lattice models where different types of initial geometric imperfections are considered. The results from these models show that imperfections have a negligible influence in the pre-buckling regime but affect the behavior considerably in the post-buckling regime. As a final step in this work, we subject the lattice models to impact loading and employ hypothetical soft and stiff materials. These results show that the lattice stiffness, which depends on material stiffness, strut diameter, and orientation, significantly influences the ability of a lattice structure to resist impact. In particular, we find that a stiffer lattice (i.e., one made with a stiff material and thicker struts) is capable of absorbing more energy than a softer one during impact. Although material nonlinearities, inelasticity, and detailed contact formulations are not considered in this study, the findings obtained herein lay the groundwork for engineering lattice structures under extreme loading conditions through a simulation-driven framework based on particle-based methods.

97 MATHEMATICS AND COMPUTING↗