Body Geometry Research Articles

3D CFD analysis of pressure, boundary layer and shear stresses on a gudgeon (Gobio gobio)

The fish’s mechanosensory lateral line system detects non-acoustic hydrodynamic stimuli required for feeding, schooling, predator avoidance and underwater object detection. Biological investigations have established that flow stimuli are detected through the boundary layer as pressure gradients by canal neuromasts and as shear stresses acting on the superficial neuromasts. Previous works have also shown that the spatial distribution of neuromasts is strongly correlated with the pressure coefficient. Despite these fundamental insights, substantial knowledge gaps persist in understanding how fish body geometry influences the boundary layer, the pressure distribution and shear stresses. To address these gaps, we provide a set of numerical models based on the open-source CFD toolkit OpenFOAM which are experimentally validated using velocity measurements obtained in a laboratory fish swim tunnel. Specifically, we investigate the mid dorsal-ventral planar flow fields around a 3D fish-shaped body of gudgeon (Gobio gobio), a common freshwater bottom-dwelling fish. The contributions of this work are two-fold: First, we provide a comparison of the boundary layer thicknesses and velocity profiles at flow velocities ranging from 0.25 to 1.25 m/s. Second, we qualitatively compare the spatial distributions of the pressure coefficient, dynamic pressure and shear stresses to biological observations of the neuromast locations of adult gudgeon.

Three-Dimensional Subsurface Model of Luk-Ulo Melange Complex, Karangsambung, Indonesia: Insights from Gravity Modeling

The Luk-Ulo Melange Complex (LMC) is characterized by a chaotic assemblage of mixed rocks with a block-in-matrix fabric. The exposed blocks consist of various scattered rock types, trending in an ENE-WSW direction. In the case of Mt. Parang, the origin of the diabase remains uncertain, with ongoing debate as to whether it is associated with in situ volcanic activity or represents an exotic block within the melange deposit. Subsurface data obtained through geophysical investigation can aid in modeling the geometry of intrusive bodies using inverse modeling techniques. In this study, we conducted a gravity survey and performed 3D inverse modeling to investigate the subsurface beneath Karangsambung. A total of 818 gravity data points and 28 rock density measurements were integrated with existing geological data to construct an a priori 3D geological model. To ensure the results align with geological concepts, the 3D inversion utilized a stochastic approach, allowing for the incorporation of multiple geological constraints over fifty million iterative procedures. Ultimately, the inversion successfully reduced the misfit between observed and calculated data from 2.71 to 0.55 mGal. Based on the inverted 3D model, the diabase rock in Mt. Parang is identified as having an intrusive origin. The intrusion model exhibited minimal changes in density, volume, and shape during the inversion process. Additionally, the model suggests the presence of a solidified magma reservoir at a depth of approximately 3 km, potentially related to Dakah volcanism. The inverted model also reveals the block-in-matrix structure of the Luk-Ulo Melange Complex in the northern area.

Impact of Vehicle Engine Supply System on In-Service Changes in Body Geometry and Wheel Alignment

Impact of Vehicle Engine Supply System on In-Service Changes in Body Geometry and Wheel Alignment

Dyke emplacement under mixed loading conditions: Insights from the Dharwar Craton, India

Dykes are intrusive igneous bodies that play crucial role in the supply and ascent of magma to the Earth’s crust. Magma can intrude along pre-existing anisotropies such as fractures or foliations present within the host rock or it may create its own path by fracturing the host rock. In the latter scenario, when fractures are formed by the pressure exerted by the invading magma, a dyke’s outcrop shape and geometry are diagnostic of the conditions under which it evolved. Here, we report mafic dykes emplaced within the younger granites of Dharwar Craton, peninsular India. Outcrop attributes of these dykes are characteristic of emplacement under conditions of mixed mode loading. We discuss different discrete modes of fracture formation and their possible combinations to understand the generation and eventual emplacement of dykes under mixed mode loading. This leads to the development of a comprehensive sequence of progressive dyke evolution under mixed mode I-III loading and thereby distinguishing incremental orders of dyke horn formation. We further apply this knowledge along with collected field evidence on dyke body geometries to propose an evolutionary model of dyke formation and emplacement within the Chitradurga granite under varying regional stress fields of the Chitradurga Schist Belt, Western Dharwar Craton, India. We infer, that the dykes initiated as extensional fractures within an earlier NE-SW directed compressive stress field and were subsequently sheared sinistrally by the effect of the adjacent Chitradurga Shear Zone on account of a later E-W to ESE-WNW directed compression.

Stiffness in compression therapy: Analytical estimation of pressure changes beneath textile compression devices

Compression pressure changes in dynamic conditions under textile compression devices have a critical impact on the success of compression therapy. Models exist to predict the level of compression pressure, but not the actual change in pressure. This paper aims to derive a formula to accurately determine the pressure change under textile compression devices, to investigate the factors influencing the pressure change and to verify their effects through theoretical analysis. Firstly, a formula based on Laplace’s law is presented which mathematically describes the dependencies of pressure changes. Secondly, a simulation is carried out to demonstrate the effect of these dependencies on pressure changes using theoretical textile curve functions. Finally, the effect of these dependencies is demonstrated by testing short- and long-stretch bandages in a tensile testing machine and using the recorded material curves to simulate a theoretical application of these bandages to the lower limbs. The results show that the change in pressure is not solely determined by the intrinsic properties of the material, but is influenced by several variables, including the mechanical performance of the textile materials during stretching, the target pressures for application of the textile material, and the body geometries to which the material is applied. Pressure change cannot be a constant for textile compression devices such as bandages. The research increases the understanding of the factors that influence pressure changes in compression device materials. The findings may have implications for the design and selection of compression textiles in clinical applications.

Evaluating thermal sensation in outdoor environments - Different methods of coupling CFD and radiation modelling with a human body thermoregulation model

Thermal comfort is vital in outdoor spaces for social interaction and recreation. With increases in heat wave events and global warming, it becomes more essential to be able to assess heat stresses. Outdoor spaces have asymmetrical and transient environmental factors such as solar radiation, long-wave radiation and natural wind, which makes the thermal perception outdoors largely different from indoors. Traditional thermal comfort models have some limitations and less accuracy when applied to outdoor spaces. The current study focuses on developing a method to assess the human thermal interaction with the complex outdoor environment and the consequent thermal sensation. Computational Fluid Dynamics (CFD) and radiation modelling are coupled to the Joint System Thermoregulation (JOS-3) model to assess the skin temperature, based on which thermal sensation is predicted using the updated thermal comfort model developed by UC Berkeley. The proposed methodology is applied to assess the thermal sensation of human subjects standing in front of an outdoor radiant cooling hub. The important finding is that the one-way coupling, i.e., no skin temperature feedback to the CFD simulation, works reasonably well in the outdoor environment, and that convective heat transfer coefficients available from experimental studies can be coupled with CFD simulated environmental parameters to obtain more reliable convective heat exchange to avoid meshing the geometry of a human body in CFD and further reduce the computing cost.

Improving the energy harvesting performance of clamped flexible plates using vortex shedding behind bluff bodies: A computational study

This study assesses flow-induced vibrations of clamped flexible plates with the objective of improving their energy harvesting performance. Toward this, a rectangular bluff body is placed between the two clamped flexible plates to harness the vortices generated behind the bluff body. The strain energy of the plate is used as a measure of energy harvesting performance. Computational studies are performed for different parameters of interest, such as dimensions and material properties of the plate and dimensions and locations of the bluff body. The effects of these stated parameters on flow-induced vibration response and vortical structures are investigated, and the optimal values for the location and geometry of the bluff body, as well as the aspect ratio and Young's modulus of the plate for energy harvesting performance, are determined. The results show that vortex shedding from the bluff body strongly influences the dynamic behavior and energy output of the flexible structures inside the wake region of the bluff body at different locations. Additionally, the aspect ratio and its effect on vorticity and energy harvesting are discussed in detail, along with the average displacement and average lift force experienced by the plates. The outcomes of this work offer significant insights into optimizing the design of clamped flexible plates for optimal energy generation by cleverly exploiting the vortex shedding behind fixed bluff bodies.

Seismic amplification of Late Quaternary paleovalley systems: 2D seismic response analysis of the Pescara paleovalley (Central Italy)

Robust site characterization and ground response analysis require a thorough understanding of subsurface features, including geophysical properties and geometries of sediment bodies. Late Quaternary paleovalley systems, often overlooked in seismic hazard assessments, represent a potential threat due to their unconsolidated infill (with shear wave velocities <200 m/s) and sharp contrast with the adjacent substrate. Through an integrated approach that combined geophysical and stratigraphic data, we characterized the subsurface of the Pescara paleovalley system. Geostatistical interpolation of microtremor measurements enabled mapping resonance frequencies, highlighting abrupt changes and delineating the paleovalley boundaries. High-resolution core descriptions were then correlated with resonance frequencies, enabling the reconstruction of a 3D geophysical depth model of the buried paleovalley morphology. Furthermore, analyzing velocity profiles from down-hole tests led to the identification of five main seismic/stratigraphic layers within the valley fill. The geometry and facies architecture were reconstructed through a cross-section transversal to the paleovalley axis and then implemented into a 2D finite element model. Seismic response was computed, revealing significant amplification factors at frequencies closely matching the direct observations. Amplification factors peaked at frequencies between 0.9 and 1.1 Hz in the paleovalley center and up to 5.5 Hz towards the flanks, reaching a factor of 4.6. These findings suggest a notable increase in amplification amplitude compared to simpler geological contexts and emphasize the potential impact on common building types. Response spectra show strong amplifications in the paleovalley system, potentially leading to an underestimation of spectral accelerations compared to NTC18 guidelines. The comparisons of 1D and 2D modeling approaches revealed minimal differences, indicating that the generally flat geometry of the valley may not exhibit clear 2D effects. However, local subsurface stratigraphy strongly influences lateral changes in seismic response, emphasizing the importance of detailed subsurface knowledge for realistic seismic response estimates.

A fresh look at the Gulf of St. Lawrence, Eastern Canada: Geophysical imaging of salt bodies and volcanic rocks in the upper Paleozoic Magdalen Basin

In the Gulf of St. Lawrence, the late Middle Devonian to early Permian Magdalen Basin hosts numerous salt bodies, but their number and geometry remained poorly constrained due to the poor quality of industry-seismic data. Salt originates from the Middle Mississippian (Visean) Windsor Group, which includes marine carbonates and evaporites. In this study, we use multibeam bathymetry as well as unpublished aerogravity and aeromagnetic data to better characterize the geometry of salt bodies. These datasets show a high degree of coherency and illustrate deformation on the seafloor, gravity lows and magnetic highs, all associated with salt bodies and overlying volcanic rocks. In map view, the geometry of individual salt bodies varies from nearly circular to ovoid and strongly elongated. Some salt bodies coalesce in complex-shaped salt walls elongated NNE, parallel to some of the faults that controlled the early stages of the Magdalen Basin development. Forward modeling indicates that low density salt bodies overlain by a 100–300 m thick volcanic pile may account for both gravity and magnetic signals. The ubiquity of mafic volcanic rocks indicates that the late Visean to Serpukhovian volcanic episode is more voluminous than previously thought. This episode is possibly associated with a series of dykes continuous for up to 225 km and with a period of underplating now recorded by a lower crustal layer with high seismic velocity.

Space–time isogeometric analysis of tire aerodynamics with complex tread pattern, road contact, and tire deformation

AbstractThe space–time (ST) computational method “ST-SI-TC-IGA” and recently-introduced complex-geometry isogeometric analysis (IGA) mesh generation methods have enabled high-fidelity computational analysis of tire aerodynamics with near-actual tire geometry, road contact, tire deformation, and aerodynamic influence of the car body. The tire geometries used in the computations so far included the longitudinal and transverse grooves. Here, we bring the tire geometry much closer to an actual tire geometry by using a complex, asymmetric tread pattern. The complexity of the tread pattern required an updated version of the NURBS Surface-to-Volume Guided Mesh Generation (NSVGMG) method, which was introduced recently and is robust even in mesh generation for complex shapes with distorted boundaries. The core component of the ST-SI-TC-IGA is the ST Variational Multiscale (ST-VMS) method, and the other key components are the ST Slip Interface (ST-SI) and ST Topology Change (ST-TC) methods and the ST Isogeometric Analysis (ST-IGA). They all play a key role. The ST-TC, uniquely offered by the ST framework, enables moving-mesh computation even with the topology change created by the contact between the tire and the road. It deals with the contact while maintaining high-resolution flow representation near the tire.The computational analysis we present is the first of its kind and shows the effectiveness of the ST-SI-TC-IGA and NSVGMG in tire aerodynamic analysis with complex tread pattern, road contact, and tire deformation.

IMPROVING APPROACHES TO DETERMINING THE LEVEL OF FILLING OF THE PIPE BODY WITH PUMPED PRODUCT

As of 2015, there were more than 200 permanent and variable gravity sections of crude oil and petroleum product pipelines in the pipeline system of PJSC Transneft, the natural monopoly for pipeline transportation of crude oil and petroleum products in the Russian Federation. The operation of gravity sections in oil pipelines in Russia is a very important safety and environmental issue, and Russian oil transportation companies are constantly working on improving approaches to monitor and control pipeline operation modes. When these sections occur, special attention is given to determining the geometric capacity of the pipeline section, as pumping is done in batches and commercial pumping product remains in the gravity section during periodic shutdowns. Due to the limited technical equipment of these sections and the likely illegal unrecorded product withdrawals in these pipeline sections, the most demanding approach is to determine the fill factor of the pipe body. The fill factor is influenced by many factors, including the presence of gas phase and water accumulations, correct geometry of the pipe body and other factors. Given the given initial data, it can be seen that the solution to the problem of determining the filling level of a stopped oil pipeline is reduced to the construction of a multifactor model. For accurate determination, in addition to the model, there is a need for relevant software, an algorithm for solving the problem ofcalculating the filling factor. Correspondingly, operating organisations need a system for accounting and control of the product filling level at these sections at the level of improvement of the state system for ensuring unity of measurements and metrological support of accounting for the long-distance pipeline transport of oil and oil products - a strategic resource, the revenues from the sale of which account for about 23 % of the budgetary system of the Russian Federation.The calculation created and proposed by the authors is divided into stages in which the initial data is sequentially collected, the voltages on the capacitors are calculated, the data is processed using the interpolation method, an image is constructed using the Newton-Gauss method, and the filling factor is calculated with various multiphase pumping products. The calculation is made based on voltage readings from sensors located orthogonally to the axis of the pipe body. Using the approximation method and the Newton-Gauss method, you can get rid of the «noise» and level the phase boundary. When tested and implemented at pipeline pumping facilities, it will be possible to clarify the existing dependences of the influence of gas phase friction on the fluid flow rate, and, therefore, evaluate the effect on the filling factor. This calculation can become a tool in the fight against illegal tappings.

Integrated GPR analysis and numerical modelling for discovering archaeological mysteries at Kom Ombo temple

ABSTRACT Archaeological research is made accurate and detailed by non-invasive geophysical techniques. Hence, the thorough utilisation of GPR data sets alongside numerical modelling is imperative for reducing uncertainty surrounding archaeological sites. This study aims to assess the practicality and efficacy of GPR numerical modelling in analysing concealed archaeological elements within the renowned Kom Ombo temple in southern Egypt. Using 80 MHz and 270 MHz antennas to detect any buried structures on the temple grounds. From the processed 2D GPR sections, two distinct anomalies (deep and shallow) were identified, possibly indicating hidden archaeological structures. Given the likelihood of these features possessing complex geometries, four synthetic models were devised, maintaining consistent electrical properties while varying in buried body geometry, to evaluate the GPR response. Convolution perfectly matched layer (CPML) absorbing media was designed to absorb waves within the modelling grid, alongside a transverse magnetic (TM) mode formulation, were utilised. The findings suggest that the deep anomaly could potentially represent a tunnel or significant geological formation, while the shallow anomaly may correspond to an extension of the Turkish fort. Overall, numerical modelling emerges as an enhanced mapping tool for uncovering hidden archaeological features in areas worldwide where excavation via traditional methods is not feasible.

Combined spin orientation and phase function of asteroids

Context. Large sky surveys provide numerous non-targeted observations of small bodies of the Solar System. The upcoming LSST of the Vera C. Rubin Observatory will be the largest source of small body photometry in the next decade. With non-coordinated epochs of observation, colors – and therefore taxonomy and composition – can only be computed by comparing absolute magnitudes obtained in each filter by solving the phase function (evolution of brightness of the small body against the solar phase angle). Current models in use in the community (HG, HG12*, and HG1G2), however, fail to reproduce the long-term photometry of many targets due to the change in the aspect angle between apparitions. Aims. We aim to derive a generic yet simple phase function model accounting for the variable geometry of the small bodies over multiple apparitions. Methods. As a spinoff of the HG1 G2 model, we propose the sHG1G2 phase function model in which we introduce a term describing the brightness changes due to spin orientation and polar oblateness. We applied this new model to 13 245 908 observations of 122 675 Solar System objects (SSOs). These observations were acquired in the g and r filters with the Zwicky Transient Facility between November 1, 2019 and December 1, 2023. We retrieved them and implemented the new sHG1G2 model in FINK, a broker of alerts designed for the LSST. Results. The sHG1G2 model leads to smaller residuals than other phase function models, providing a better description of the photometry of asteroids. We determined the absolute magnitude, H, and phase function coefficients (G1, G2) in each filter, the spin orientation (α0, δ0), and the polar-to-equatorial oblateness, R, for 95 593 SSOs, which constitutes about a tenfold increase in the number of characterized objects compared to the current census. Conclusions. The application of the sHG1 G2 model to ZTF alert data using the FINK broker shows that the model is appropriate for extracting physical properties of asteroids from multi-band and sparse photometry, such as the forthcoming LSST survey.

Inversion of marine magnetotelluric and controlled source electromagnetic data with site gap and low signal-to-noise ratio

SUMMARY Marine magnetotelluric (MT) and controlled source electromagnetic (CSEM) methods have been routinely applied to survey the crustal and upper mantle electrical resistivity structures beneath the sea floor. In practice, there are inevitably site gaps and contamination by noises in the collected data because of lost ocean bottom electromagnetic receivers, unusable data, and difficulties in deploying instruments near deep trenches. So far, it remains unclear to what degree those factors will lower the resolution and the credibility of marine MT and CSEM inversion models. In this paper, we investigate the individual and combined effects of site gaps and data noises on the inversion models through synthetic analyses based on a simple block resistivity model and a realistic resistivity structure derived from the Mariana Trench. The results suggest that data with a sufficiently high signal-to-noise ratio can reasonably recover the subseafloor structures in the area of data gap. The transverse electric mode and tipper data from the MT method are much more sensitive to the structure near the site gap. The joint inversion of MT and CSEM data would improve the model's resolution at the site gap area. The inversion of data with a relatively low signal-to-noise ratio, for example, 10 per cent, can recover the structures with few artefacts if there is no site gap. But if the site gap and noisy data are combined, even a joint inversion cannot correctly recover the burial depths and geometries of the anomalous bodies beneath the site gap where vertical strips are likely present. To improve the model's resolution and suppress inversion artefacts, we propose constraining part of the model with as much a prior information as possible. Specifically, for a survey in the subduction zone, we could reduce the penalties on the model's smoothness at the upper and low interfaces of the resistive subduction slab, or even fix the resistivity of the resistive slab with the help of other information, if any. The inversion models shown in this paper provide valuable references for the site design before marine MT and CSEM surveys as well as for interpreting real data inversion models that may be subject to the same biases introduced by the site gap and noise.

Anisotropic Mach Cone Aligned Mesh Adaptation for Low Boom Simulations

An automated off-body Mach cone aligned structured curvilinear grid generation procedure is presented for near-field computational fluid dynamics simulations. This procedure combines output-based indicators and mesh redistribution to perform anisotropic mesh adaptation while maintaining Mach cone alignment. Automation is achieved through a novel direction-based adaptation indicator formulation. The adaptation procedure is demonstrated on the JAXA Wing Body geometry and X-59 C608 demonstrator models from the Second and Third AIAA Sonic Boom Prediction Workshops, respectively. It is demonstrated that anisotropic mesh adaptation may result in a greater than 50% reduction in the resource usage required to achieve the same level of accuracy as uniform and user-constructed Mach cone aligned grids for near-field pressure signatures, ground-level overpressure signatures, and loudness metrics.

Body schema: Resolving the conundrum of the distorted body

Humans show perceptual biases that suggest distorted internal representations of their own body. New research reveals that these perceptual biases can reflect integration of prior assumptions about body posture rather than a misshaped representation of the body's geometry.

Biases in hand perception are driven by somatosensory computations, not a distorted hand model

To sense and interact with objects in the environment, we effortlessly configure our fingertips at desired locations. It is therefore reasonable to assume that the underlying control mechanisms rely on accurate knowledge about the structure and spatial dimensions of our hand and fingers. This intuition, however, is challenged by years of research showing drastic biases in the perception of finger geometry.1,2,3,4,5 This perceptual bias has been taken as evidence that the brain's internal representation of the body's geometry is distorted,6 leading to an apparent paradox regarding the skillfulness of our actions.7 Here, we propose an alternative explanation of the biases in hand perception-they are the result of the Bayesian integration of noisy, but unbiased, somatosensory signals about finger geometry and posture. To address this hypothesis, we combined Bayesian reverse engineering with behavioral experimentation on joint and fingertip localization of the index finger. We modeled the Bayesian integration either in sensory or in space-based coordinates, showing that the latter model variant led to biases in finger perception despite accurate representation of finger length. Behavioral measures of joint and fingertip localization responses showed similar biases, which were well fitted by the space-based, but not the sensory-based, model variant. The space-based model variant also outperformed a distorted hand model with built-in geometric biases. In total, our results suggest that perceptual distortions of finger geometry do not reflect a distorted hand model but originate from near-optimal Bayesian inference on somatosensory signals.

Neural Garment Dynamics via Manifold‐Aware Transformers

AbstractData driven and learning based solutions for modeling dynamic garments have significantly advanced, especially in the context of digital humans. However, existing approaches often focus on modeling garments with respect to a fixed parametric human body model and are limited to garment geometries that were seen during training. In this work, we take a different approach and model the dynamics of a garment by exploiting its local interactions with the underlying human body. Specifically, as the body moves, we detect local garment‐body collisions, which drive the deformation of the garment. At the core of our approach is a mesh‐agnostic garment representation and a manifold‐aware transformer network design, which together enable our method to generalize to unseen garment and body geometries. We evaluate our approach on a wide variety of garment types and motion sequences and provide competitive qualitative and quantitative results with respect to the state of the art.

Numerical Investigation of Butterfly Valve Performance in Variable Valve Sizes, Positions and Flow Regimes

Reliability and efficiency of valves are necessary for precise control and sufficient heat-flow to heat application plants for the integrated energy systems of nuclear power plants (NPPs). Strategic Management Analysis Requirement and Technology (SMART) valves’ ability to control flow and assess environmental parameters stands out for these requirements. Their ability to sustain the downstream flow rate, prevent reverse flow, and maintain pressure in the heat transport loop is much more efficient with the integration of sensors and intelligent algorithms. For assessing valve performance and monitoring, mechanical design and operating conditions are two important parameters. In this study, the butterfly valves of three different sizes are simulated with water and steam using STAR-CCM+ in various flow regimes and positions to analyze performance parameters to strategize an automated control system for efficiently balancing the heat–transport network. Also, flow behavior is studied using velocity and pressure fields for valve–body geometry optimization. It can be observed, through performance parameters, that the valves are suitable for operation between 30° and 90° positions with significantly low loss coefficients and high flow coefficients, and the performance parameters follow a certain pattern in both water and steam flow in each scenario.

Capacitive biophysical stimulation improves the healing of vertebral fragility fractures: a prospective multicentre randomized controlled trial

BackgroundCapacitively coupling electric fields (CCEF) is a method of non-invasive biophysical stimulation that enhances fracture repair and spinal fusion. This multicentre randomized controlled trial aimed to further examine the roles of CCEF in (1) the resolution of vertebral bone marrow oedema (VBME) using a follow-up MRI study and (2) pain relief, analgesic drug consumption and quality of life improvement in stimulated patients who were referred with acute vertebral fragility fractures (VFFs) compared to non-stimulated patients.MethodsBetween September 2016 and December 2019, patients who were referred to the spine centres that participated in this multicentre randomized clinical study with acute VFFs of type OF1 or OF2 were included in the present study. All the VFFs were conservatively managed according to Good Clinical Practice. Moreover, the patients were randomized into two groups: the CCEF group received, as an adjunct to the clinical study protocol, biophysical stimulation with a CCEF device (Osteospine, IGEA) for 8 h per day for 60 days, whereas the control group was treated according to the clinical study protocol. At baseline (T0), the 30-day follow-up (T1), the 60-day follow-up (T2), and the 6-month follow-up (T3), each patient underwent clinical evaluation using the Visual Analogue Scale (VAS) for Pain and the Oswestry Disability Index (ODI). Analgesic therapy with paracetamol 1000 mg tablets for 7 days—or longer, depending on the pain intensity—was performed; patients were required to report their paracetamol consumption on a specific sheet between study day 8 to 180 days of follow-up. MRI studies of the thoracolumbar spine were performed at 0 (T0), 30 (T1) and 60 days of follow-up (T2) using a 1.5-T MRI system in all of the centres that took part in the study. For each VBME area examined via MRI, the vertebral body geometry (i.e. anterior wall height/posterior wall height and vertebral kyphosis) were assessed.ResultsA total of 66 patients (male: 9, 13.63%; mean age: 73.15 years old) with 69 VFFs were included in the present study and randomized as follows: 33 patients were included in the control group and the remaining 33 patients were randomized into the CCEF group. In the CCEF group, good compliance with CCEF therapy was observed (adherence = 94%), and no adverse effects were recorded. In the stimulated patients, faster VBME resolution and significantly less vertebral body collapse during follow-up were observed compared to the control patients. Moreover, in the active group, faster pain reduction and improvement in the ODI mean score were observed. Stimulated patients also reported a significantly lower paracetamol consumption rate from the third follow-up after treatment until the 6-month follow-up. In terms of sex-related differences, in the CCEF group, VBME showed a faster resolution in male patients compared with females.ConclusionBiophysical stimulation with CCEF, as an adjunct to traditional conservative treatment, is a useful tool to hasten the VBME resolution process and prevent vertebral body deformation. These MRI findings also correlate with faster back pain resolution and quality of life improvement. From the third follow-up after treatment untilthe 6-month follow-up, stimulated patients reported a significantly lower paracetamol consumption than control patients, even though back pain and quality of life showed no significant differences between the two groups.Level of evidenceII.Trial Registration Register: ClinicalTrials.gov, number: NCT05803681.