Numerical Analysis Research Articles

Numerical Analysis of Spoilers and Chamfered Corners for Mitigating Wind Loads on Low-Rise Flat Roof Buildings

Numerical Analysis of Spoilers and Chamfered Corners for Mitigating Wind Loads on Low-Rise Flat Roof Buildings

An experimental and numerical analysis of slit plate as replaceable fuse in moment resisting frame

This study contributes to the ongoing progress in developing cost-effective structural fuse solutions for steel structures, emphasizing enhanced strength, ductility, and energy dissipation capacity. These fuses are strategically placed to coincide with plastic hinge locations under lateral load. During minor seismic events, they provide substantial rotational stiffness, while during medium to large earthquakes, they leverage their post-yield strength for effective energy dissipation. Three varieties of beam-column sub-assemblages viz. single oval slit damper fuse (SOSDF), dual oval slit damper fuse (DOSDF) and single elliptical slit damper fuse (SESDF), featuring advanced replaceable fuses were evaluated using both experimental and numerical analyses under quasi-static loading conditions. The seismic performance of these specimens was comprehensively assessed, focusing on characteristics such as moment-carrying capacity concerning rotation, energy dissipation, ductility, and failure modes. The findings indicate that all specimens exhibited stable hysteresis, with connections utilizing SESDF and DOSDF showing slightly higher moment carrying capacities (19 % and 10 %, respectively) compared to those employing SOSDF. According to AISC 2016 standards, all three connections were classified as full-strength connections. Moreover, the ductility of specimens incorporating SESDF was notably higher as 3.5, marking a 28 % and 26 % increase over the ductility of SOSDF and DOSDF connections, respectively. In major seismic events, damage was confined solely to the fuse components, with no damage observed in the beam and column. Detailed finite element analysis using ABAQUS CAE further corroborated these findings, demonstrating good alignment between experimental and simulated results.

Experimental and numerical study on dynamic response of a photovoltaic support structural platform with a U-shaped tuned liquid column damper

A tuned liquid column damper (TLCD) is widely used in offshore structures as a passive energy dissipation device to reduce the harm of wind, wave, and current loads to the safety of offshore platforms. This investigation explores the dynamic response and interaction mechanism of a photovoltaic support structural platform (SSP) equipped with a TLCD by experimental and numerical analysis. The vibration mitigation performance of the TLCD under varying liquid depth ratios, external excitation frequencies, and amplitudes is systematically evaluated. The experiments focus on assessing the structural dynamic response and wall pressure to determine the vibration reduction efficiency of the TLCD for structures with fundamental frequencies more than 2 Hz. The experimental results indicate that the TLCD can reduce the structural response maximally near the natural frequency (2.52 Hz) of the structure, with vibration reduction effectiveness ranging from 42.0% to 85.1%. Additionally, one develops a two-way fluid-structure interaction (FSI) model to investigate further the impact of hydrodynamic pressure caused by sloshing inside the TLCD on the SSP motion response, which providing a more detailed explanation for the different phenomena observed in the experimental results.

Enhancing (quasi-)long-range order in a two-dimensional driven crystal.

It has been recently shown that 2D systems can exhibit crystalline phases with long-range translational order showcasing a striking violation of the Hohenberg-Mermin-Wagner (HMW) theorem, which is valid at equilibrium. This is made possible by athermal driving mechanisms that inject energy into the system without exciting long wavelength modes of the density field, thereby inducing hyperuniformity. However, as thermal fluctuations are superimposed on the non-equilibrium driving, long-range translational order is inevitably lost. Here, we discuss the possibility of exploiting non-equilibrium effects to suppress arbitrarily large density fluctuations even when a global thermal bath is coupled to the system. We introduce a model of a harmonic crystal driven both by a global thermal bath and by a momentum conserving noise, where the typical observables related to density fluctuations and long-range translational order can be analytically derived and put in relation. This model allows us to rationalize the violation of the HMW theorem observed in previous studies through the prediction of large-wavelength phonons, which thermalize at a vanishing effective temperature when the global bath is switched off. The conceptual framework introduced through this theory is then applied to numerical simulations of a hard-disk solid in contact with a thermal bath and driven out-of-equilibrium by active collisions. Our numerical analysis demonstrates how varying driving and dissipative parameters can lead to an arbitrary enhancement of the quasi-long-range order in the system regardless of the applied global noise amplitude. Finally, we outline a possible experimental procedure to apply our results to a realistic granular system.

Vibration-based damage detection of beams using Hybrid Genetic Algorithm with combined l0 and l1 regularization

Vibration-based damage detection is an effective way of identifying beam damages before they become severe and potentially catastrophic. This paper introduces a novel and robust method for detecting multiple cracks in beams using a Hybrid Genetic Algorithm (HGA). The objective function of the proposed method, based on dynamic properties, can be employed with partial modal information and is enhanced by including the initial residuals between numerical and experimental models along with combined l0 and l1 regularization. This enhancement mitigates the impact of limited modal data, experimental noise, and modeling inaccuracies. Additionally, a method of selecting an appropriate regularization parameter by analyzing residual and solution norm curves is presented. The performance of the proposed method is evaluated and verified using both numerical simulations and experimental studies on a steel cantilever beam under four damage scenarios. In the experimental study, natural frequencies and mode shapes were determined under ambient vibration conditions using Enhanced Frequency Domain Decomposition (EFDD) and the Stochastic Subspace Identification (SSI) method. Finite element models of the beam were developed in ABAQUS software for numerical analysis, and a comparative study of numerical and experimental data was conducted. The proposed approach successfully identifies, locates and quantifies single and multiple damages in the beam, demonstrating its reliability and effectiveness, particularly in scenarios with limited modal data and significant experimental noise.

Two-dimensional flat-band solitons in superhoneycomb lattices

Abstract Flat-band periodic materials are characterized by a linear spectrum containing at least one band where the propagation constant remains nearly constant irrespective of the Bloch momentum across the Brillouin zone. These materials provide a unique platform for investigating phenomena related to light localization. Meantime, the interaction between flat-band physics and nonlinearity in continuous systems remains largely unexplored, particularly in continuous systems where the band flatness deviates slightly from zero, in contrast to simplified discrete systems with exactly flat bands. Here, we use a continuous superhoneycomb lattice featuring a flat band in its spectrum to theoretically and numerically introduce a range of stable flat-band solitons. These solutions encompass fundamental, dipole, multi-peak, and even vortex solitons. Numerical analysis demonstrates that these solitons are stable in a broad range of powers. They do not bifurcate from the flat band and can be analyzed using Wannier function expansion leading to their designation as Wannier solitons. These solitons showcase novel possibilities for light localization and transmission within nonlinear flat-band systems.

Porosity estimation based on the shear modulus inversion of seismic shear wave

Reliable estimation of subsurface porosity is necessary for hydrocarbon reservoir characterization and fluid identification. Porosity estimations from seismic data can provide the lateral distribution of subsurface porosity, but the results may be highly nonunique because subsurface elastic properties (such as velocity and density) can be affected by porosity and pore fluids. Because shear modulus is insensitive to pore fluid content, it can be effectively used to estimate porosity. We develop a novel porosity estimation method using the reflected SH wave (SH-SH wave), whose propagation characteristics depend mainly on shear modulus and S-wave velocity. We derive a new analytical expression for the SH-SH-wave reflection coefficient, which acts as a function of shear modulus and S-wave velocity in its natural logarithm form. This new approximation has high accuracy, and both coefficients corresponding to shear modulus and S-wave velocity are “model-parameter independent”; thus, there is no need for prior estimation of any model parameter during inversion. Numerical analysis indicates that shear modulus inverted from the SH-SH wave has lower uncertainty, the problem is better conditioned, and the method requires data with fewer incidence angles than algorithms that invert for the shear modulus from the PP wave. Furthermore, the highly correlated rock-physics relationships between porosity and shear modulus facilitate accurate porosity estimation. A field data application indicates that high-resolution porosity of fine structures can be reliably recovered using our method.

Master Curves for Poroelastic Spherical Indentation With Step Displacement Loading

Abstract Theoretical and numerical analyses are conducted to rigorously construct master curves that can be used for interpretation of displacement-controlled poroelastic spherical indentation test. A fully coupled poroelastic solution is first derived within the framework of Biot’s theory using the McNamee–Gibson displacement function method. The fully saturated porous medium is assumed to consist of slightly compressible solid and fluid phases and the surface is assumed to be impermeable over the contact area and permeable everywhere else. In contrast to the cases in our previous studies with an either fully permeable or impermeable surface, the mixed drainage condition yields two coupled sets of dual integral equations instead of one in the Laplace transform domain. The theoretical solutions show that for this class of poroelastic spherical indentation problems, relaxation of the normalized indentation force is affected by material properties through weak dependence on a single-derived material constant only. Finite element analysis is then performed in order to examine the differences between the theoretical solution, obtained by imposing the normal displacement over the contact area, and the numerical results where frictionless contact between a rigid sphere and the poroelastic medium is explicitly modeled. A four-parameter elementary function, an approximation of the theoretical solution with its validity supported by the numerical analysis, is proposed as the master curve that can be conveniently used to aid the interpretation of the poroelastic spherical indentation test. Application of the master curve for the ramp-hold loading scenario is also discussed.

Particle manipulation based on the controllable light field of graphene nano-butterfly structures

In this Letter, we have proposed a nano-optical convey belt consisting of three sets of graphene nano-butterfly structures (GNBS) with different lengths arranged in a periodic pattern to trap, transport, and sort particles depending on controlling the Fermi level of graphene. The graphene nano-butterfly structure, based on a couple of trapezoids, can form a potential well to capture particles by localized surface plasmon resonance under the excitation of an incident laser. There are three adjustable geometric parameters containing width (W1, W2) and length (L) that make the GNBS have a high degree of design flexibility. Graphene can support surface plasmons in the mid-infrared region, and at the same time, it possesses a more flexible tuning mechanism compared to metals. When the incident light power is 1 mW/μm2 and the resonance wavelength is 10.6 μm, the depth of the potential well is greater than 10KbT for 50 nm particles. Moreover, the sorting of particles with different diameters has also been validated under appropriate power and switching time. Based on the capabilities of capture and sorting verified by numerical analysis, our design offers a new scheme for on-chip optofluidic applications.

Investigations into dynamic variation characteristics of near-bottom boundary flows over a long term at strait areas

Dynamic processes at the bottom boundary layers are crucial for various marine applications such as marine energy exploration and environmental protection. This paper explored hydrodynamic characteristics of near-bottom boundary currents based on the continuous measured velocity data at a site near the bottom boundary layer in the Strait of Georgia. Previous studies on the bottom boundary layer encompasses various aspects, including physical, chemical, and biological dimensions. However, the mechanism underlying the hierarchical structure of flow velocity remain poorly understood. Thus, the power spectrum and wavelet analysis were used to analyze the significant period of the averaged velocity and tidal characteristics, and to calculate the thickness of near-bottom boundary layer. The results indicate that the current velocities at the boundary layer tend to flow in south-eastern and north-western directions. There is a significant period of 15 days in August and September. At a depth of 301 m, the vertical velocities are stratified, indicating different velocity layers. The boundary layer thickness varies from 5.23 to 14.74 m, as indicated by the vertical structural characteristics. The structural characteristics of the seabed boundary layer are conducive to deepening our understanding of sediment's incipient motion, transport, and sedimentation process. The findings provide a theoretical foundation for future studies on sediment movement and the numerical simulation analysis of the bottom boundary layer. It is of great significance for the study of wave flow interaction and marine dynamic structure, and can provide scientific support for the development and utilization of the strait area.

Power Forecasting for Photovoltaic Microgrid Based on MultiScale CNN-LSTM Network Models

Photovoltaic (PV) microgrids comprise a multitude of small PV power stations distributed across a specific geographical area in a decentralized manner. Computational services for forecasting the output power of power stations are crucial for optimizing resource deployment. This paper proposes a deep-learning-based architecture for short-term prediction of PV power. Firstly, in order to make full use of the spatial information between different power stations, a spatio–temporal feature fusion method is proposed. This method is capable of exploiting both the power information of neighboring power stations with strong correlations and meteorological information with the PV feature data of the target power station. By using a multiscale convolutional neural network–long short-term memory (CNN-LSTM) network model, it is capable of generating a PV feature dataset containing spatio–temporal attributes that expand the data source and enhance the feature constraints. It is capable of predicting the output power sequences of power stations in PV microgrids with high model generalization and responsiveness. To validate the effectiveness of the proposed framework, an extensive numerical analysis is also conducted based on a real-world PV dataset.

Performance Prediction of GFRP-Reinforced Concrete Deep Beams Containing a Web Opening in the Shear Span

This study aimed to investigate the nonlinear structural behavior of concrete deep beams internally reinforced with glass fiber-reinforced polymer (GFRP) reinforcing bars and containing a web opening of various sizes and locations within the shear span. Three-dimensional (3D) numerical simulation models were developed for large-scale GFRP-reinforced concrete deep beams (300 mm × 1200 mm × 5000 mm) with a shear span-to-depth ratio (a/h) of 1.04. Predictions of the numerical models were validated against published experimental data. A parametric study was conducted to examine the effect of varying the opening size and location on the shear response. Results of the numerical analysis indicated that the strength of the deep beam models with an opening in the middle of the shear span decreased with an increase in either the opening width or height. The rate of the strength reduction caused by increasing the opening height was, however, more significant than that produced by increasing the opening width. Placing a web opening in the compression zone close to the load plate was very detrimental to the beam strength. Conversely, a negligible strength reduction was recorded when the web opening was placed in the tension side above the flexural reinforcement and away from the natural load path. Data of the parametric study were utilized to introduce simplified analytical formulas capable of predicting the shear capacity of GFRP-reinforced concrete deep beams with a web opening in the shear span.

On the Challenges of Integrating a Rotating Detonation Combustor with an Industrial Gas Turbine and Important Design Considerations for Row-1 Blades

Abstract In an effort to increase the efficiency and performance of gas turbine power cycles, pressure gain combustion (PGC) has gained significant interest. Since Rotating Detonation Combustors (RDC) can provide a quasi-steady mode of operation, research has been triggered to integrate RDC with power-generating gas turbines. However, the presence of subsonic and supersonic flow fields which are generated due to the shock waves that stem from the detonation wave front and the highly non-uniform temperature and velocity profiles may cause a depreciation in the turbine performance. The current study seeks to investigate the challenges of integrating the RDC with nozzle guide vanes (NGV) of an industrial, can-annular gas turbine and attempts to understand the major contributors that impact efficiency and identify the key areas of optimization that need to be considered for maximizing performance. The RDC was integrated with the NGVs through a non-optimized straight duct-type geometry with a diffuser cone. 3-Dimensional Numerical analyses were performed to investigate sources of total pressure loss and to understand the unsteady effects of RDC which contribute towards the deterioration of performance. The entropy generation at different regions of interest was calculated to identify the major irreversibilities in the system. Also, total pressure and temperature distribution along the radial direction at the exit of the transitional duct is presented.

Numerical analysis of magnetohydrodynamic (MHD) convection using CNT-based nanofluids in a wavy-shaped enclosure with a double lid

Numerical analysis of magnetohydrodynamic (MHD) convection using CNT-based nanofluids in a wavy-shaped enclosure with a double lid

A contact-based constitutive model for the numerical analysis of masonry structures using the distinct element method

This study presents a robust contact constitutive model in the distinct element method (DEM) framework for simulating the mechanical behavior of masonry structures. The model is developed within the block-based modeling strategy, where the masonry unit is modeled as deformable blocks with potential crack surfaces in the middle of the bricks, while the mortar joints are defined as zero-thickness interfaces. The modeling strategy implements multi-surface plasticity with damage mechanics, including a tension cut-off, Coulomb failure criterion, and an elliptical compressive cap for the damage in tension, shear, and compression, respectively. Two new features are introduced in this contact model: a piecewise linear softening function for strength degradation in tension and shear and a hardening/softening function to phenomenologically define the compressive damage of masonry composite into the unit-mortar interface. The constitutive model is implemented in commercial DEM software using the small displacement configuration and validated against material and experimental tests on masonry walls subjected to constant pre-compression and monotonically increasing in-plane load. The experimental and numerical results regarding the force-displacement relationship and damage pattern produced by the proposed constitutive model are compared and critically discussed, demonstrating the capability of DEM coupled with the suitable constitutive law in simulating the behavior of masonry structures.

Uncovering Pattern-Transformable Soft Granular Crystals Induced by Microscopic Instability

Abstract Upon compression, some soft granular crystals undergo pattern transformation. Recent studies have unveiled that the underlying mechanism of this transformation is closely tied to microscopic instability, resulting in symmetry breaking. This intriguing phenomenon gives rise to unconventional mechanical properties in the granular crystals, paving the way for potential metamaterial application. However, no consistent approach has been reported for studying other unexplored transformable granular crystals. In this study, we present a systematic approach to identify a new set of pattern-transformable diatomic granular crystals induced by microscopic instability. After identifying the kinematic constraints for diatomic soft granular crystals, we have generated a list of feasible particle arrangements for instability-induced pattern transformation under compression. Instead of computationally intensive finite element models (FEMs) with continuum elements, we adopt a simplified mass-spring model derived from granular contact networks to efficiently evaluate these feasible particle arrangements for pattern transformation. Our numerical analysis encompasses quasi-static analysis and microscopic/macroscopic instability analyses within the framework of linear perturbation. Subsequently, the pattern transformation of the identified particle arrangements is confirmed through quasi-static analyses employing detailed finite element (FE) simulations with continuum elements. Additional numerical simulations with continuum elements reveal that the pattern transformations of particle arrangements are significantly influenced by the initial void volume and some transformed granular crystals may exhibit strong low-frequency directional phononic band-gaps, which were not observed in the initial granular crystals.

Enhancing cerebral arteriovenous malformation analysis: Development and application of patient-specific lumped parameter models based on 3D imaging data

Objectives:Cerebral arteriovenous malformations (AVMs) present complex neurovascular challenges, characterized by direct arteriovenous connections that disrupt normal brain blood flow dynamics. Traditional lumped parameter models (LPMs) offer a simplified angioarchitectural representation of AVMs, yet often fail to capture the intricate structure within the AVM nidus. This research aims at refining our understanding of AVM hemodynamics through the development of patient-specific LPMs utilizing three-dimensional (3D) medical imaging data for enhanced structural fidelity. Methods:This study commenced with the meticulous delineation of AVM vascular architecture using threshold segmentation and skeletonization techniques. The AVM nidus’s core structure was outlined, facilitating the extraction of vessel connections and the formation of a detailed fistulous vascular tree model. Sampling points, spatially distributed and derived from the pixel intensity in imaging data, guided the construction of a complex plexiform tree within the nidus by generating smaller Y-shaped vascular formations. This model was then integrated with an electrical analog model to enable precise numerical simulations of cerebral hemodynamics with AVMs. Results:The study successfully generated two distinct patient-specific AVM networks, mirroring the unique structural and morphological characteristics of the AVMs as captured in medical imaging. The models effectively represented the intricate fistulous and plexiform vessel structures within the nidus. Numerical analysis of these models revealed that AVMs induce a blood shunt effect, thereby diminishing blood perfusion to adjacent brain tissues. Conclusion:This investigation enhances the theoretical framework for AVM research by constructing patient-specific LPMs that accurately reflect the true vascular structures of AVMs. These models offer profound insights into the hemodynamic behaviors of AVMs, including their impact on cerebral circulation and the blood steal phenomenon. Further incorporation of clinical data into these models holds the promise of deepening the theoretical comprehension of AVMs and fostering advancements in the diagnosis and treatment of AVMs.

Unveiling Mispricing Risks: Nonlarge Homogeneous Portfolio Factor Copula Models for Enhanced Valuation of Subordinated Loan Securitization

AbstractThis paper presents an innovative factor copula model for collateralized loan obligation (CLO) tranche valuation, incorporating non‐Gaussian distributions and dynamic correlations without relying on the large homogeneous portfolio (LHP) assumption. Through numerical analyses and comparisons with LHP models, I find that non‐LHP models produce higher tranche spreads, especially for lower‐rated tranches. Sensitivity analysis reveals varying sensitivities to changes in the number of collaterals, risk‐free rate, average collateral ratings, recovery rates, and time to maturity. The non‐LHP one‐factor copula models, including stochastic correlation and random factor loading models, outperform LHP models in root mean squared errors when calibrated to market data. The results underscore the importance of considering model limitations in CLO tranche pricing and highlight potential mispricing of spread risk in higher‐rated tranches using LHP models. The proposed models contribute to a more comprehensive understanding of CLO tranche pricing by accounting for various factors and assumptions influencing fair premiums.

Local buckling behaviour of Q690 high strength steel cold-formed round-ended oval hollow section stub columns under axial compression: Experimental investigation, numerical modelling and design

This study presents a comprehensive investigation of Q690 high strength steel (HSS) cold-formed round-ended oval hollow section (ROHS) stub columns subjected to concentric compression, utilising a combination of experimental and numerical analyses. Material properties were determined through tensile coupon tests. Residual stresses in both membrane and bending types were measured and recorded and their impact on structural behavior was evaluated. Initial local geometric imperfection measurements were conducted for all specimens. In addition to experimental studies, the numerical investigation was conducted to develop a finite element (FE) model of Q690 HSS cold-formed ROHS stub column. The accuracy of the FE model for HSS cold-formed ROHS stub column was validated through a comparison of the numerical predictions with the experimental test results. An extensive parametric study was conducted to produce structural performance data for HSS cold-formed ROHS with an expanded range of cross-sectional slenderness, aiming to enhance the existing data repository. The cross-sectional resistance of stub columns obtained from experimental and numerical investigations were compared with the design strengths predicted by the advanced design methods of Direct Strength Method (DSM) and Continuous Strength Method (CSM). The comparison indicate that the existing design methods provide over-conservative and unreliable cross-sectional resistance predictions. Modifications on DSM and CSM are proposed, based upon which notable improvements in cross-sectional resistance in predictions are obtained in a reliable manner.

Postbuckling behaviour of anisotropic shear deformable laminated cylindrical shells with general imperfection distribution subjected to torsion

Composite cylindrical shells belong to typical basic structure in engineering, whose buckling and postbuckling behaviours determine the stability and safety of overall structures under complex service loadings. In this paper, a nonlinear mechanical model for anisotropic laminated composite cylindrical shells with general distributed initial imperfection is derived to analyse the buckling and postbuckling behaviours under torsional loads, in which the general distributed imperfections are represented by a generic two-dimensional imperfection function reconstructed by using measured data. Based on a higher order shear deformation shell theory, the governing equations are established with consideration of von Kármán–Donnell-type of kinematic nonlinearity and the coupling effects of extension/twist, extension/flexural and flexural/twist. By using a two-step perturbation technique, a new set of explicit solutions is obtained to determine the buckling loads and postbuckling equilibrium paths of the torsional composite cylindrical shells. Furthermore, the internal quantitative relationship between the shear stress with torsional characteristics along with an associate compressive stress of anisotropic laminated shell is for the first time obtained. Numerical analyses are performed to validate the accuracy and effectiveness of the present explicit solutions, and obtain the postbuckling response of different types of imperfections, anisotropic laminated cylindrical shells with different values of shell parameters.