Nonlinear Simulation Research Articles

Splitting of triply quantized vortices in Bose-Einstein condensates of finite temperature

Abstract Utilizing the dissipative Gross-Pitaevskii equation, we investigated the splitting dynamics of triply quantized vortices at finite temperature. Through linear perturbation analysis, \textcolor{red}{we determined the excitation modes of these vortices across various dissipation parameters. We identified three unstable modes with $p=2, 3, $ and $4$ fold rotational symmetries, } revealing a significant dynamic transition of the most unstable mode. That is, as the dissipation parameter increases, the most unstable mode transitions from $p=2$ mode to $p=3$ mode. Throughout the entire range of dissipation parameters, the $p=4$ unstable mode is never the dominant mode. Subsequently, we performed nonlinear numerical simulations of the vortex splitting process. Under random perturbations, we confirmed the dynamical transition, and under specific perturbations, we confirmed the instability of the $p=4$ mode. Our findings on the finite temperature dependence of the splitting dynamics of triply quantized vortices are expected to be verifiable in experiments.

Anand Model Constants of Sn-Ag-Cu Solders: What Do They Actually Mean?

Abstract Finite element simulations are broadly utilized to calculate the thermally induced mechanical deformations of interconnecting solder materials in electronics. The Anand unified viscoplasticity model is a prevalent choice for simulating the mechanical responses of solder joints, comprising nine parameters whose individual effects are not fully addressed in the literature. For this reason, this study examines the influence of each Anand parameter on the mechanical response of leadless Tin-Silver-Copper solder alloys with varying silver content through comprehensive nonlinear finite element simulations. Anand model coefficients for different SAC alloys, including SAC105, SAC205, SAC305, and SAC405, were sourced from existing literature and used to establish a systematic study matrix for each parameter. These coefficients were then integrated into thermomechanical simulations to induce inelastic deformations in the solder interconnections. The response of the solder interconnects is analyzed using equivalent inelastic strains and inelastic strain energy density. Results indicated that specific Anand parameters could skew the solder joint stress-strain response towards brittle behavior, while other parameters could lead to a more ductile response. Through statistical factorial analysis, the significance of each parameter is found to vary considerably, ranging from negligible to highly significant. The findings of this study are crucial for understanding the behavior of different SAC solder configurations and their expected thermal fatigue performance based on Anand creep constants. Furthermore, this paper provides foundational insights into the interpretation of Anand coefficients and their influence on the mechanical response of solder materials, a topic not yet explored in the existing literature.

The role of ion-scale micro-turbulence in pedestal width of the DIII-D wide-pedestal QH mode

Abstract The low-edge rotation, intrinsically ELM-free, and improved confinement wide-pedestal quiescent H-mode (QH-mode), discovered in DIII-D tokamak, has pedestal widths exceeding the EPED-KBM model scaling typically by at least 25%. Ion-scale (k_y ρ_s<1) microturbulence and its role in setting the pedestal structure is investigated using the radially local δf gyrokinetic code CGYRO. The electromagnetic trapped electron mode (TEM) is unstable at the pedestal top, while plasma beta (β_e) is ~60% below the kinetic-ballooning mode (KBM) onset threshold and the electron temperature gradient (ETG) mode is found to be unstable in the peak gradient region. Nonlinear simulation reveals that the ion-scale turbulence could produce electron energy flux consistent with the flux inferred from power balance at the pedestal top, with a reasonable variation of the local E×B shearing rate; and the local neoclassical transport from NEO is dominant over the simulated turbulent transport in the ion energy flux channel. The simulated ion-scale turbulence produces much lower electron energy flux than inferred from experiment in the pedestal peak gradient region. A correction to the EPED-KBM pedestal width scaling is obtained based on the two-dimensional scan of pedestal top plasma beta (β_e) and normalized electron density and temperature scale lengths, a⁄L_(n_e ) , a⁄L_(T_e ) using CGYRO linear simulations. Mode transitions among TEM, micro-tearing mode (MTM), ion-temperature gradient (ITG) mode and KBM, are observed in the 2D scan at the pedestal top. A fixed normalized growth rate for these drift-type modes is taken to determine the pedestal width scaling, which shows good consistency with the QH experimental database on pedestal heights and widths. The onset of KBM instabilities and the local ExB shear suppression criterion set the lower and upper limit for the pedestal width of standard QH-mode, wide-pedestal QH-mode and type-I ELMy H mode. A potentially higher and wider pedestal is expected from the new scaling of pedestal width. This work presents an improved understanding of the ion-scale micro-turbulence of wide-pedestal QH-mode and sheds light on a promising scenario for future reactors, including ITER and beyond.

Prediction of performance and turbulence in ITER burning plasmas via nonlinear gyrokinetic profile prediction

Abstract Burning plasma performance, transport, and the effect of hydrogen isotope (H, D, D-T fuel mix) on confinement has been predicted for ITER baseline scenario (IBS) conditions using nonlinear gyrokinetic profile predictions. Accelerated by surrogate modeling (Rodriguez-Fernandez et al 2022 Nucl. Fusion 62 076036), high fidelity, nonlinear gyrokinetic simulations performed with the CGYRO code (Candy et al 2016 J. Comput. Phys. 324 73), were used to predict profiles of Ti , Te , and ne while including the effects of alpha heating, auxiliary power (NBI + ECH), collisional energy exchange, and radiation losses inside of r / a = 0.9. Predicted profiles and resulting energy confinement are found to produce fusion power and gain that are approximately consistent with mission goals ( P fusion = 500 MW at Q = 10) for the baseline scenario and exhibit energy confinement that is within 1σ of the H-mode energy confinement scaling. The power of the surrogate modeling technique is demonstrated through the prediction of alternative ITER scenarios with reduced computational cost. These scenarios include conditions with maximized fusion gain and an investigation of potential resonant magnetic perturbation (RMP) effects on performance with a minimal number of gyrokinetic profile iterations required (3–6). These predictions highlight the stiff ITG nature of the core turbulence predicted in the ITER baseline and demonstrate that Q > 17 conditions may be accessible by reducing auxiliary input power while operating in IBS conditions. Prediction of full kinetic profiles allowed for the projection of hydrogen isotope effects around ITER baseline conditions. The gyrokinetic fuel ion species was varied from H, D, and 50/50 D-T and kinetic profiles were predicted. Results indicate that a weak or negligible isotope effect will be observed to arise from core turbulence in IBS conditions. The resulting energy confinement, turbulence, and density peaking, and the implications for ITER operations will be discussed.

A magnetohydrodynamic mechanism for the formation of solar polar vortices

Polar vortices are ubiquitous features of planetary atmospheric flows, from the Earth-like rocky planets to Jupiter- and Saturn-like gas giant planets. Very little is known about their existence or dynamics on the Sun. What should be expected near the Sun’s pole for the upcoming solar multi-viewpoint and polar missions? Here, we report the magnetohydrodynamic (MHD) nonlinear simulations for the formation and evolution of solar polar vortices using a near-surface MHD shallow-water model. Our findings indicate that the rush to the poles, the migration of magnetic fields toward the pole following the Sun’s magnetic cycle, can positively contribute to the formation of polar vortices. The mechanism proposed here for the formation of polar vortices involves the role of magnetic fields and may be relevant to any star with a magnetic cycle. The Sun’s polar vortices resulting from this mechanism are predominantly MHD, consisting of a tight pair of cyclonic and anticyclonic swirls. This mechanism is likely to operate during all solar cycle phases except the peak, when the polar field reverses. Polar vortices can impact dynamical evolution of global flows and polar fields, which seed the next activity cycle, hence better knowledge of physics of polar regions may lead to improved solar cycle and space weather forecasts.

Numerical study of turbulent eddy self-interaction in tokamaks with low magnetic shear. Part I: Linear simulations

Abstract This study examines the impact of low and zero magnetic shear, along with rational values of safety factor, on linear mode behavior in tokamak plasmas. We focus on how magnetic field line topology and parallel domain length influence linear mode structures and growth rates as magnetic shear decreases. Our results
show that as magnetic shear is reduced to low values, linear modes can transition between different branches, significantly affecting their characteristics. At zero magnetic shear, accurate modeling requires precisely capturing magnetic topology, as small deviations near rational surfaces can greatly impact growth rates. These findings are extended to nonlinear simulations in a companion paper (Volčokas, 2024), revealing that long turbulent eddies and magnetic field line topology play a crucial role in turbulence self-interaction and plasma transport.

Experimental study and finite element analysis on the shear performance of steel-bamboo composite cellular beam

To study the influence of opening ratio on the shear performance of steel-bamboo composite cellular beams (SBCC), taking the opening size as the basic parameter, experimental research was carried out on five test beams. The failure modes and deformation characteristics of the specimens were observed, and the influence and law of the change of hole size on the shear performance of steel-bamboo composite cellular beams were analyzed. Through nonlinear finite element simulation calculation, multi-parameter analysis was carried out on the composite cellular beams, and the contribution of flanges to the shear performance of cellular composite beams was studied. The results show that the integrity of steel-bamboo composite cellular beams is good. The failure is mainly manifested as the fracture of bamboo web at the beam hole area, the buckling of section steel and the debonding and separation of steel-bamboo interface. The opening ratio has a great influence on the shear performance of steel-bamboo composite cellular beams. The thickness of bamboo flange cannot be ignored for the shear contribution of cellular beams. The shear capacity of cellular beams is shared by flanges and webs. The larger the opening ratio, the greater the shear contribution of flanges to composite cellular beams.

An affine formulation of eigenstrain-based homogenization method and its application to polycrystal plasticity

Abstract Reduced order models (ROMs) are typically incorporated into concurrent multiscale approaches to allow for efficient nonlinear multiscale simulations and to alleviate high cost of direct nonlinear computational homogenization schemes. ROMs based on the ideas of transformation field analysis are among the most popular in the literature since they only require linear elastic simulations for model construction and typically have low number of degrees of freedom. However, these models have been shown to deliver overly stiff response in simulating wide range of materials. The present study focuses on mitigating this problem in the context of eigenstrain-based homogenization method (EHM) using instantaneous moduli information for polycrystal elastoviscoplasticity. For this purpose, a new EHM model is developed with the intention of using affine moduli for recomputation of the instantaneous localization tensors. The accuracy of the method is compared to the original EHM and direct crystal plasticity finite element simulations for several synthetic polycrystal microstructures, loading conditions and varying phase contrast. We show that the affine model delivers consistently softer response compared to the original EHM model. In particular, the affine model delivers notably more accurate response in the presence of high phase contrast. The affine EHM is able to capture local load redistribution through recomputation of the localization tensors.

Static actuator-sharing algorithm for concurrent control of multiple plasma properties

Abstract Simultaneous regulation of multiple properties in next-generation tokamaks like ITER and Fusion Pilot Plant (FPP) may require the integration of different plasma control algorithms. Such integration requires the conversion of individual controller commands into physical actuator requests while accounting for the coupling between different plasma properties. This work proposes a tokamak and scenario-agnostic actuator-sharing algorithm (ASA) to perform the above-mentioned command-request conversion and, hence, integrate multiple plasma controllers. The proposed algorithm implicitly solves a quadratic programming (QP) problem formulated to account for the saturation limits and the relation between the controller commands and physical actuator requests. Since the constraints arising in the QP program are linear, the proposed ASA is highly computationally efficient and can be implemented in the tokamak plasma control system (PCS) in real time. Furthermore, the proposed algorithm is designed to handle real-time changes in the control objectives and actuators’ availability. Nonlinear simulations carried out using the Control Oriented Transport SIMulator (COTSIM) illustrate the effectiveness of the proposed algorithm in achieving multiple control objectives simultaneously.

Numerical study of turbulent eddy self-interaction in tokamaks with low magnetic shear. Part II: Nonlinear simulations.

Abstract Using local nonlinear gyrokinetic simulations and building up on linear results from a companion paper (Vol\v{c}okas, 2024), we demonstrate that turbulent eddies can extend along magnetic field lines for hundreds of poloidal turns in tokamaks with weak or zero magnetic shear $\hat{s}$. We observe that this parallel eddy length scales inversely with magnetic shear and at $\hat{s}=0$ is limited by the thermal speed of electrons $v_{th,e}$. We examine the consequences of these ``ultra long" eddies on turbulent transport, in particular, how field line topology mediates strong parallel self-interaction. Our investigation reveals that, through this process, field line topology can strongly affect transport. It can cause transitions between different turbulent instabilities and in some cases triple the logarithmic gradient needed to drive a given amount of heat flux. We also extensively study a novel ``eddy squeezing" effect introduced in (Vol\v{c}okas, 2023), which reduces the perpendicular size of eddies and their ability to transport energy, thus representing a novel approach to improve confinement. Finally, we investigate the triggering mechanism of Internal Transport Barriers (ITBs) using low magnetic shear simulations, shedding light on why ITBs are often easier to trigger where the safety factor has a low-order rational value.

Sliding plane formalism for aeroacoustic and adjoint-based sensitivity calculations

This paper demonstrates a methodology for time-domain and time-accurate nonlinear, direct and adjoint simulations of unsteady flows and aeroacoustics for multi-component systems in relative motion. Here, the principal effort is directed towards mitigating the problem of distortion and contamination of the adjoint field at the moving interface, through a computationally lightweight, high-order sliding plane approach for which the adjoint equivalent is simple to obtain. This effort requires an attentive treatment of the interface conditions that surpasses the requirements of the more common forward (primary) problem. Sensitivity of a given quantity of interest from a time-varying flow with respect to a large number of parameters is then obtained through the adjoint operator, which is evaluated using nonlinear-adjoint looping. This technique is implemented using checkpointing and the PETSc TSAdjoint library and, after validation, applications including a rotor–stator interaction problem are presented.

Dosing strategy of tacrolimus when co-administered with Wuzhi tablet in renal transplant recipients.

Tacrolimus (TAC) is a first-line immunosuppressant to prevent allograft rejection. Wuzhi tablet is widely used as a TAC-sparing agent in China that could significantly elevate TAC exposure. However, insufficient data support the dose recommendation of TAC when co-administered with Wuzhi. A total of 305 adult renal transplant patients with 2,541 TAC trough concentrations (C0) were enrolled for population pharmacokinetic (PPK) modeling. CYP3A5 polymorphism was genotyped, and corresponding clinical factors were recorded. Nonlinear mixed-effects modeling and Monte Carlo simulation were used for dose recommendation. PK parameters were calculated based on one-compartment model with first-order absorption and elimination. The estimated total clearance (CL/F) and volume of distribution (Vd/F) of TAC were 23.84 L/h and 1,075.96 L, respectively. Wuzhi, CYP3A5 genotype, hematocrit (HCT), and weight were found to have a significant influence on CL/F. CL/F was significantly lower in the individuals who were CYP3A5 non-expressers and received TAC together with Wuzhi. CYP3A5 genotype (expressers or non-expressers), body weight (40-80 kg), and hematocrit (20 - 40%) were selected as the specific clinical scenarios, and the starting dose of TAC ranged from 1.5 to 4.5 mg when co-administered with Wuzhi. We establish a TAC PPK model comprising Wuzhi as a covariate in renal transplant recipients and recommend an initial dose of TAC when co-administered with Wuzhi, which could provide reference for the individualized regimens of TAC.

Rational mutual interactions in starch-lipid-chlorogenic acid ternary systems enable formation of ordered starch structure for anti-digestibility during hot-extrusion 3D printing: Based on the nonlinear rheology and molecular simulation

Rational mutual interactions in starch-lipid-chlorogenic acid ternary systems enable formation of ordered starch structure for anti-digestibility during hot-extrusion 3D printing: Based on the nonlinear rheology and molecular simulation

Experimental and numerical investigation on hyperelastic sealing disc contact behavior in pipeline, a comparison between fluid-driven and pull-through approaches

Pipelines are widely recognized as the safest and most efficient means of fluid transportation. Pipeline pigging, a secure and reliable technique, is employed for both cleaning and sealing pipelines to optimize efficiency. Since pigging operation requires differential pressure (ΔP) to facilitate pig movement, this pressure must correspond to the pipeline's operational conditions. Consequently, studying the required ΔP for sealing discs, which serve as primary sealing elements of pigs, is crucial. The primary focus of this study is to examine the difference between pulling the pig using a cable and propelling it with fluid. Furthermore, the main novelty of this research is to experimentally investigate a single sealing disc and conduct fluid-driven tests on it. To study more precisely, three sealing discs with various thicknesses and the same hardness have been launched into a 6-inch spool test, containing five pipes with different wall thicknesses, resulting in multiple oversize ratios (%OSz). In the experimental study, both fluid-driven and pull-through tests were conducted. The range of discrepancy between two methods varies from 15% to 26% for different oversize ratios which is considerable. Additionally, A 2-D axisymmetric nonlinear numerical simulation was conducted in a finite element software ABAQUS in order to study the behavior of sealing discs. Using a pressure-dependent friction coefficient with the proper hyperelastic model was key to achieving simulations that have good agreements with experimental results, with discrepancy of less than 10 percent in all extracted pressures.

Verification of electromagnetic simulation capabilities in global gyrokinetic particle-in-cell code GTS

Recently, the numerical scheme presented by Mishchenko et al. [Phys. Plasmas 21, 052113 (2014); 21, 092110 (2014)] enabled explicit gyrokinetic simulations of low-frequency electromagnetic instabilities in tokamaks at experimentally relevant values of plasma β. This scheme resolved the long-standing cancellation problem that previously hindered gyrokinetic particle-in-cell code simulations of magnetohydrodynamic phenomena with inherently small parallel electric fields. Moreover, the scheme did not employ approximations that eliminate critical tearing-type instabilities. Here, we report on the implementation of this numerical scheme in the global gyrokinetic particle-in-cell code GTS. This implementation allows for a more complete and accurate picture of interaction between small scale turbulence and MHD modes in tokamaks. Additionally, we present a comprehensive set of verification simulations of numerous electromagnetic instabilities relevant to present-day tokamaks. These simulations encompass the kinetic ballooning mode, the internal kink mode, the tearing mode, the micro-tearing mode, and the toroidal Alfven eigenmode destabilized by energetic ions, which are all instrumental in understanding tokamak physics. We will also showcase the preliminary nonlinear simulations of kinetic ballooning instabilities and (2,1) island formation due to tearing mode instability. These simulations validate the accuracy of the scheme implementation and pave the way for studying how these instabilities affect plasma confinement and performance.

Improved General Relativity Search Algorithm (IGRSA) for Designing Power System Stabilizers

This paper proposes a novel optimization algorithm for designing power system stabilizers (PSSs). The prospect of attaining higher stability motivated the authors to creat a new optimization algorithm for this study. A novel algorithm named the Improved General. Relativity Search Algorithm (IGRSA) was also developed by cloud theory to improve the performance of GRSA. The supremacy of the proposed approach is tested by comparing it with an introduced objective function in a medium multi-machine power system. The nonlinear simulation results and eigenvalues analysis has demonstrated that the proposed approach in this study is highly effective in enhancing the dynamic stability of the power system.

Onset of global instability in a premixed annular V-flame

We investigate self-excited axisymmetric oscillations of a lean premixed methane–air V-flame in a laminar annular jet. The flame is anchored near the rim of the centrebody, forming an inverted cone, while the strongest vorticity is concentrated along the outer shear layer of the annular jet. Consequently, the reaction and vorticity dynamics are largely separated, except where they coalesce near the flame tip. The global eigenmodes corresponding to the linearised reacting flow equations around the steady base state are computed in an axisymmetric setting. We identify an arc branch of eigenmodes exhibiting strong oscillations at the flame tip. The associated eigenvalues are robust with respect to domain truncation and numerical discretisation, and they become destabilised as the Reynolds number increases. The frequency of the leading eigenmode is found to correspond to the Lagrangian disturbance advection time from the nozzle outlet to the flame tip. The essential role of this convective mechanism is also supported by resolvent analysis, which finds that the same flame-tip disturbance structure and frequency are optimally amplified when the flame is subjected to external white noise forcing. Strong non-modal effects in the form of pseudo-resonance are not found. Nonlinear time-resolved simulation further reveals notable hysteresis phenomena in the subcritical regime prior to instability. Hence, even when the flame is linearly stable, perturbations of sufficient amplitude can trigger limit-cycle oscillations and higher-dimensional dynamics sustained by nonlinear feedback. A Monte Carlo simulation of passive tracers in the unsteady flame suggests a nonlinear non-local instability mechanism. Notably, linear analysis of the subcritical time-averaged limit-cycle state yields eigenvalues that do not match the nonlinear periodic oscillation frequencies. This mismatch is attributed to the fundamentally nonlinear dynamics of the subcritical V-flame instability, where the dichromatic, non-local interaction between the heat release rate along the flame surface and the vortex dynamics in the jet shear layer cannot be approximated as a simple distortion of the mean flow.

Developing Programs for Converting MIDAS GEN to ANSYS Models Based on Python

The reasonableness and accuracy of engineering design are often assessed through the use of a variety of structural design analysis software, which are then compared and verified. However, it is challenging for a single analysis software to meet the diverse and complex design requirements. In order to meet the specific engineering requirements, it is necessary to convert the MIDAS result model into an ANSYS structural model and conduct a nonlinear analysis and simulation in ANSYS. Nevertheless, the existing interface is unable to facilitate direct conversion of the model. Accordingly, this paper presents a Python-based ANSYS APDL program that enables the complete conversion of MIDAS GEN structural models to ANSYS finite element models. The program is capable of converting a range of data, including material, section, element, connection, load, node mass, constraint, time history function, and so forth. The program is capable of converting specific connection units, including elastic and general connection units. Additionally, the beam-column section direction, beam end freedom release, rigid element, and special anti-rocking structure of the structure can be considered. Ultimately, the theatre model is transformed. Following a comparison of the analysis results, it was found that the mass and mode of the model before and after the transformation were essentially identical. The maximum error of the first six orders of the structure is 2.95%, with the structural displacement under gravity load remaining essentially unchanged. The research and analysis demonstrate the accuracy and reliability of the MIDAS GEN conversion ANSYS program. The conversion program significantly reduces the time required for direct modeling in ANSYS, enhancing work efficiency. The study has considerable practical significance for the seismic sway design and analysis of buildings based on vibration isolation design.

A Nonlinear Approach in the Quantification of Numerical Uncertainty by High-Order Methods for Compressible Turbulence with Shocks

This is a comprehensive overview on our research work to link interdisciplinary modeling and simulation techniques to improve the predictability and reliability simulations (PARs) of compressible turbulence with shock waves for general audiences who are not familiar with our nonlinear approach. This focused nonlinear approach is to integrate our “nonlinear dynamical approach” with our “newly developed high order entropy-conserving, momentum-conserving and kinetic energy-preserving methods” in the quantification of numerical uncertainty in highly nonlinear flow simulations. The central issue is that the solution space of discrete genuinely nonlinear systems is much larger than that of the corresponding genuinely nonlinear continuous systems, thus obtaining numerical solutions that might not be solutions of the continuous systems. Traditional uncertainty quantification (UQ) approaches in numerical simulations commonly employ linearized analysis that might not provide the true behavior of genuinely nonlinear physical fluid flows. Due to the rapid development of high-performance computing, the last two decades have been an era when computation is ahead of analysis and when very large-scale practical computations are increasingly used in poorly understood multiscale data-limited complex nonlinear physical problems and non-traditional fields. This is compounded by the fact that the numerical schemes used in production computational fluid dynamics (CFD) computer codes often do not take into consideration the genuinely nonlinear behavior of numerical methods for more realistic modeling and simulations. Often, the numerical methods used might have been developed for weakly nonlinear flow or different flow types other than the flow being investigated. In addition, some of these methods are not discretely physics-preserving (structure-preserving); this includes but is not limited to entropy-conserving, momentum-conserving and kinetic energy-preserving methods. Employing theories of nonlinear dynamics to guide the construction of more appropriate, stable and accurate numerical methods could help, e.g., (a) delineate solutions of the discretized counterparts but not solutions of the governing equations; (b) prevent numerical chaos or numerical “turbulence” leading to FALSE predication of transition to turbulence; (c) provide more reliable numerical simulations of nonlinear fluid dynamical systems, especially by direct numerical simulations (DNS), large eddy simulations (LES) and implicit large eddy simulations (ILES) simulations; and (d) prevent incorrect computed shock speeds for problems containing stiff nonlinear source terms, if present. For computation intensive turbulent flows, the desirable methods should also be efficient and exhibit scalable parallelism for current high-performance computing. Selected numerical examples to illustrate the genuinely nonlinear behavior of numerical methods and our integrated approach to improve PARs are included.

Robustness of Dry‐Assembled Precast Wall and Coupled Wall‐Frame High‐Rise Structures With Voided Prestressed Slabs

ABSTRACTPrecast concrete structures are recognized among the structural typologies that may mostly suffer from progressive collapse as a consequence of local unpredicted damage, also due to past accidents occurred in buildings employing old precast technologies neglecting the basic robustness criteria and characterized by poorly detailed and/or executed joints and connections. Although vertical ties may relatively easily be provided in both frame and wall panel structures by employing ductile connection devices, dry‐assembled precast systems avoiding concrete pouring in the whole superstructure typically struggle to provide horizontal ties as strong as required by the current standards concerning progressive collapse, in which they are not fully framed yet. Nevertheless, these systems may exploit alternative sources of robustness to stop collapse progression arising from their structural arrangement, the shape of the slab elements employed, and the slab mechanical connections. This article presents the results of a numerical investigation focused on the structural behavior of a dry‐assembled precast floor system consisting of prestressed box‐section slab elements employed in either wall panel or coupled wall‐frame structural systems following the loss of one or multiple primary load‐bearing vertical elements. Nonlinear static simulations of different loss scenarios occurring in prototypes of 100‐m‐tall high‐rise buildings were carried out based upon spread plasticity modeling of the structural elements and experimentally calibrated nonlinear behavior of mechanical connections.