To save content items to your account,
please confirm that you agree to abide by our usage policies.
If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your account.
Find out more about saving content to .
To save content items to your Kindle, first ensure no-reply@cambridge.org
is added to your Approved Personal Document E-mail List under your Personal Document Settings
on the Manage Your Content and Devices page of your Amazon account. Then enter the ‘name’ part
of your Kindle email address below.
Find out more about saving to your Kindle.
Note you can select to save to either the @free.kindle.com or @kindle.com variations.
‘@free.kindle.com’ emails are free but can only be saved to your device when it is connected to wi-fi.
‘@kindle.com’ emails can be delivered even when you are not connected to wi-fi, but note that service fees apply.
State-of-the-art controllers for active back exosuits rely on body kinematics and state machines. These controllers do not continuously target the lumbosacral compression forces or adapt to unknown external loads. The use of additional contact or load detection could make such controllers more adaptive; however, it can be impractical for daily use. Here, we developed a novel neuro-mechanical model-based controller (NMBC) that uses a personalized electromyography (EMG)-driven musculoskeletal (MSK) model to estimate lumbosacral joint loading. NMBC provided adaptive, subject- and load-specific assistive forces proportional to estimates of the active part of biological joint moments through a soft back support exosuit. Without a priori information, the maximum assistive forces of the cable were modulated across weights. Simultaneously, we applied a non-adaptive, kinematic-dependent, trunk inclination-based controller (TIBC). Both NMBC and TIBC reduced the mean and peak biomechanical metrics, although not all reductions were significant. TIBC did not modulate assistance across weights. NMBC showed larger reductions of mean than peak values, significant reductions during the erect stance and the cumulative compressive loads by 21% over multiple cycles in a cohort of 10 participants. Overall, NMBC targeted mean lumbosacral compressive forces during lifting without a priori information of the load being carried. This may facilitate the adoption of non-hindering wearable robotics in real-life scenarios. As NMBC is informed by an EMG-driven MSK model, it is possible to tune the timing of NMBC-generated torque commands to the exosuit (delaying or anticipating commands with respect to biological torques) to target further reduction of peak or mean compressive forces and muscle fatigue.
Low-density polymer foams pre-ionized by a well-controlled nanosecond pulse are excellent plasma targets to trigger direct laser acceleration (DLA) of electrons by sub-picosecond relativistic laser pulses. In this work, the influence of the nanosecond pulse on the DLA process is investigated. The density profile of plasma generated after irradiating foam with a nanosecond pulse was simulated with a two-dimensional hydrodynamic code, which takes into account the high aspect ratio of interaction and the microstructure of polymer foams. The obtained plasma density profile was used as input to the three-dimensional particle-in-cell code to simulate energy, angular distributions and charge carried by the directional fraction of DLA electrons. The modelling shows good agreement with the experiment and in general a weak dependence of the electron spectra on the plasma profiles, which contain a density up-ramp and a region of near-critical electron density. This explains the high DLA stability in pre-ionized foams, which is important for applications.
Ocean turbulence at meso- and submesocales affects the propagation of surface waves through refraction and scattering, inducing spatial modulations in significant wave height (SWH). We develop a theoretical framework that relates these modulations to the current that induces them. We exploit the asymptotic smallness of the ratio of typical current speed to wave group speed to derive a linear map – the U2H map – between surface current velocity and SWH anomaly. The U2H map is a convolution, non-local in space, expressible as a product in Fourier space by a factor independent of the magnitude of the wavenumber vector. Analytic expressions of the U2H map show how the SWH responds differently to the vortical and divergent parts of the current, and how the anisotropy of the wave spectrum is key to large current-induced SWH anomalies. We implement the U2H map numerically and test its predictions against WAVEWATCH III numerical simulations for both idealised and realistic current configurations.
Machine learning has become a dominant problem-solving technique in the modern world, with applications ranging from search engines and social media to self-driving cars and artificial intelligence. This lucid textbook presents the theoretical foundations of machine learning algorithms, and then illustrates each concept with its detailed implementation in Python to allow beginners to effectively implement the principles in real-world applications. All major techniques, such as regression, classification, clustering, deep learning, and association mining, have been illustrated using step-by-step coding instructions to help inculcate a 'learning by doing' approach. The book has no prerequisites, and covers the subject from the ground up, including a detailed introductory chapter on the Python language. As such, it is going to be a valuable resource not only for students of computer science, but also for anyone looking for a foundation in the subject, as well as professionals looking for a ready reckoner.
In the present research, the effect of streamwise finlets on the coherent structures of a turbulent boundary layer and their relation with pressure fluctuations and trailing-edge noise is investigated experimentally over a NACA0018 airfoil. A synthetic measurement is performed using time-resolved particle image velocimetry, wall-pressure transducers and a far-field microphone. The finlets induce strong momentum transport within the boundary layer, leading to the formation of a detached shear layer and backward flow separation. A strong velocity deficit is produced close to the wall. The instantaneous flow organisation reveals the formation of hairpin-like vortices on top of the finlets and spanwise rollers in the near-wall separation bubble. The newly generated vortices disrupt the turbulent coherent structures of the untreated case remarkably. An overall lift-up process of the unsteady turbulent structures is produced, bringing the most energetic turbulent structures away from the wall and reducing the near-wall shear stress. The spatial and temporal relation between instantaneous unsteady flow features and wall-pressure fluctuations is analysed quantitatively. A notable reduction of the correlation and coherence intensity in the mid- and high-frequency bands is achieved due to the modification of the turbulent structures. The former frequency ranges agree with that of the pressure fluctuations and far-field noise suppression, revealing the noise-reduction mechanisms.
Regression is a fundamental prediction task common in data-centric engineering applications that involves learning mappings between continuous variables. In many engineering applications (e.g., structural health monitoring), feature-label pairs used to learn such mappings are of limited availability, which hinders the effectiveness of traditional supervised machine learning approaches. This paper proposes a methodology for overcoming the issue of data scarcity by combining active learning (AL) for regression with hierarchical Bayesian modeling. AL is an approach for preferentially acquiring feature-label pairs in a resource-efficient manner. In particular, the current work adopts a risk-informed approach that leverages contextual information associated with regression-based engineering decision-making tasks (e.g., inspection and maintenance). Hierarchical Bayesian modeling allow multiple related regression tasks to be learned over a population, capturing local and global effects. The information sharing facilitated by this modeling approach means that information acquired for one engineering system can improve predictive performance across the population. The proposed methodology is demonstrated using an experimental case study. Specifically, multiple regressions are performed over a population of machining tools, where the quantity of interest is the surface roughness of the workpieces. An inspection and maintenance decision process is defined using these regression tasks, which is in turn used to construct the active-learning algorithm. The novel methodology proposed is benchmarked against an uninformed approach to label acquisition and independent modeling of the regression tasks. It is shown that the proposed approach has superior performance in terms of expected cost—maintaining predictive performance while reducing the number of inspections required.
Sea ice is a mushy layer, a porous material whose properties depend on the relative proportions of solid and liquid. The growth of sea ice is governed by heat transfer through the ice together with appropriate boundary conditions at the interfaces with the atmosphere and ocean. The salinity of sea ice has a major effect on its thermal properties so might naïvely be expected to have a major effect on its growth rate. However, previous studies observed a low sensitivity throughout the winter growth season. The goal of this study is to identify the controlling physical mechanisms that explain this observation. We develop a simplified quasi-static framework by applying a similarity transformation to the underlying heat equation and neglecting the explicit time dependence. We find three key processes controlling the sensitivity of growth rate to salinity. First, the trade-off between thermal conductivity and (latent) heat capacity leads to low sensitivity to salinity even at moderately high salinity and brine volume fraction. Second, the feedback on the temperature profile reduces the sensitivity relative to models that assume a linear profile, such as zero-layer Semtner models. Third, thicker ice has the opposite sensitivity of growth rate to salinity compared with thinner ice, sensitivities that counteract each other as the ice grows. Beyond its use in diagnosing these sensitivities, we show that the quasi-static approach offers a valuable sea-ice model of intermediate complexity between zero-layer Semtner models and full partial-differential-equation-based models such as Maykut–Untersteiner/Bitz–Lipscomb and mushy-layer models.
The lattice Boltzmann method has become a popular tool for simulating complex flows, including incompressible turbulent flows; however, as an artificial compressibility method, it can generate spurious pressure oscillations whose impact on the statistics of incompressible turbulence has not been systematically examined. In this work, we propose a theoretical approach to analyse the origin of compressibility-induced oscillations (CIOs) and explore ways to suppress or remove them. We begin by decomposing the velocity field and pressure field each into the solenoidal component and the compressive component, and then study the evolution of these two components analytically and numerically. The analysis yields an evolution equation of the mean-square pressure fluctuation which reveals several coupling effects of the two components. The evolution equation suggests that increasing the bulk-to-shear viscosity ratio can suppress CIOs, which is confirmed by numerical simulations. Furthermore, based on the derived evolution equation and data from the simulation, a model is developed to predict the long-term behaviours of the mean-square pressure fluctuations. In the case of decaying turbulence in a periodic domain, we show that the Helmholtz–Hodge decomposition can be used to obtain the solenoidal components reflecting the true evolution of incompressible turbulent flow, from the mesoscopic artificial compressibility approach. The study provides general theoretical guidelines to understand, suppress and even remove CIOs in other related pseudo-compressibility methods.
We investigate flame–acoustic interactions in a turbulent combustor during the state of intermittency before the onset of thermoacoustic instability using complex networks. Experiments are performed in a turbulent bluff-body stabilised dump combustor where the inlet airflow rate is varied quasi-statically and continuously. We construct a natural visibility graph from the local heat release rate fluctuations ($\dot {q}'$) at each location. Comparing the average degree during epochs of high- and low-amplitude acoustic pressure oscillations ($p'$) during the state of intermittency, we detect frequency modulation in $\dot {q}'$. Through this approach, we discover unique spatial patterns of cross-variable coupling between the frequency of $\dot {q}'$ and the amplitude of $p'$. The frequency of $\dot {q}'$ increases in regions of flame anchoring owing to high-frequency excitation of the flow and flame during epochs of high-amplitude $p'$ dynamics. However, the frequency of $\dot {q}'$ decreases in regions associated with flame-front distortions by large coherent vortices. In experiments with continuously varying airflow rates, the spatial pattern of frequency modulation varies with an increase in the average amplitude of $p'$ owing to an increase in the epochs of periodic $p'$ dynamics and the size of vortices forming in the flow. Dynamic shifts in the location of flame anchoring induce low-frequency fluctuations in $\dot {q}'$ during very-high-amplitude intermittent $p'$ dynamics. Our approach using conditional natural visibility graphs thus reveals the spatial pattern of amplitude-frequency coupling between the co-evolving flame and the acoustic field dynamics in turbulent reacting flows.
In this experimental and numerical study, we revisit the question of the onset of the elastic regime in viscoelastic pinch-off. This is relevant to all modern filament thinning techniques, which aim to measure the extensional properties of low-viscosity polymer solutions. Examples are the slow retraction method (SRM) for capillary breakup extensional rheometry (CaBER), or the dripping method, in which a drop detaches from a nozzle. As part of these techniques, a stable liquid bridge is brought slowly to its stability threshold, where capillary-driven thinning starts. This thinning slows down dramatically at a critical radius $h_1$, marking the onset of the elasto-capillary regime, characterised by a filament of nearly uniform radius. While a theoretical scaling exists for this transition in the case of the classical step-strain CaBER protocol, where polymer chains stretch without relaxing during the fast plate separation, we show that this theory is not necessarily valid for a slow protocol such as the SRM. In that case, polymer chains start stretching (beyond their equilibrium coiled configuration) only when the bridge thinning rate becomes comparable to the inverse of their relaxation time. We derive a universal scaling for $h_1$, valid for both low- and high-viscosity polymer solutions. This scaling is validated by CaBER experiments with a slow plate separation protocol using different polymer solutions, plate diameters and sample volumes, as well as by numerical simulations using the FENE-P model.
Novel methods of data collection and analysis can enhance traditional risk management practices that rely on expert engineering judgment and established safety records, specifically when key conditions are met: Analysis is linked to the decisions it is intended to support, standards and competencies remain up to date, and assurance and verification activities are performed. This article elaborates on these conditions. The reason engineers are required to perform calculations is to support decision-making. Since humans are famously weak natural statisticians, rather than ask stakeholders to implicitly assimilate data, and arrive at a decision, we can instead rely on subject matter experts to explicitly define risk management decision problems. The results of engineering calculation can then also communicate which interventions (if any) are considered to be risk-optimal. It is also proposed that the next generation of engineering standards should learn from the success of open source software development in community building. Interacting with open datasets and code can promote engagement, identification (and resolution) of errors, training and ultimately competence. Finally, the profession’s tradition of independent verification should also be applied to the complex models that will increasingly contribute to the safety of the built environment. Model assurance will be required to keep pace with model development to identify suitable use cases as adequately safe. These are considered to be increasingly important components in ensuring that methods of data-centric engineering can be safely and appropriately adopted in industry.
Safe and efficient flight operations depend on effective air traffic management and the decision-making skills of air traffic control officers (ATCOs). However, managing air traffic in terminal control areas (TMAs), especially in approach control units, is challenging due to the complexity of the airspace. This is particularly evident in metroplex airspaces like the Istanbul TMA, which features multiple civil and military airports, diverse approach systems, and heavy traffic volumes, all contributing to an exceptionally complex operational environment. This study examines how experienced ATCOs perceive airspace complexity, focusing on approach control units within TMAs. Using Istanbul TMA as a case study, the research combines qualitative and quantitative methods to analyse the factors contributing to complexity. In the first phase, the Content Validity Method (CVM) is used to identify and confirm the key factors influencing airspace complexity. In the second phase, the Best-Worst Method (BWM) is applied to measure the importance of these factors. The study involves two groups of ATCOs: 40 in the first group and 20 in the second. The results reveal that ‘conflicts’ are the most critical factor affecting airspace complexity, highlighting the importance of conflict resolution in air traffic control. Other significant factors include rules and procedures, airspace design and traffic density. This study provides clear insights into the challenges of managing TMA, especially in metroplex airspace. Identifying and analysing the key factors offers valuable guidance for improving air traffic management and supporting ATCOs in making better decisions.
Accurate characterization of high-power laser parameters, especially the near-field and far-field distributions, is crucial for inertial confinement fusion experiments. In this paper, we propose a method for computationally reconstructing the complex amplitude of high-power laser beams using modified coherent modulation imaging. This method has the advantage of being able to simultaneously calculate both the near-field (intensity and wavefront/phase) and far-field (focal-spot) distributions using the reconstructed complex amplitude. More importantly, the focal-spot distributions at different focal planes can also be calculated. To verify the feasibility, the complex amplitude optical field of the high-power pulsed laser was measured after static aberrations calibration. Experimental results also indicate that the near-field wavefront resolution of this method is higher than that of the Hartmann measurement. In addition, the far-field focal spot exhibits a higher dynamic range (176 dB) than that of traditional direct imaging (62 dB).
Research in lower limb wearable robotic control has largely focused on reducing the metabolic cost of walking or compensating for a portion of the biological joint torque, for example, by applying support proportional to estimated biological joint torques. However, due to different musculotendon unit (MTU) contractile speed properties, less attention has been given to the development of wearable robotic controllers that can steer MTU dynamics directly. Therefore, closed-loop control of MTU dynamics needs to be robust across fiber phenotypes, that is ranging from slow type I to fast type IIx in humans. The ability to perform closed-loop control the in-vivo dynamics of MTUs could lead to a new class of wearable robots that can provide precise support to targeted MTUs for preventing onset of injury or providing precision rehabilitation to selected damaged tissues. In this paper, we introduce a novel closed-loop control framework that utilizes nonlinear model predictive control to keep the peak Achilles tendon force within predetermined boundaries during diverse range of cyclic force production simulations in the human ankle plantarflexors. This control framework employs a computationally efficient model comprising a modified Hill-type MTU contraction dynamics component and a model of the ankle joint with parallel actuation. Results indicate that the closed-form muscle-actuation model’s computational time is in the order of microseconds and is robust to different muscle contraction velocity properties. Furthermore, the controller achieves tendon force control within a time frame below $ 18\mathrm{ms} $, aligning with the physiological electromechanical delay of the MTU and facilitating its potential for future real-world applications.
We present a higher-order spectral element method (SEM) for analyzing eccentric anisotropic multilayered waveguides. The formulation uses conformal transformation optics to map the original eccentric waveguide to an equivalent concentric one. This transformation is extended to handle anisotropic non-reciprocal media characterized by non-symmetric and non-Hermitian tensors. In the transformed concentric domain, higher-order two-dimensional basis functions associated with the zeros of the completed Lobatto polynomial are used to expand the fields. To simulate radially unbounded problems, we employ a complex-stretched perfectly matched layer to mimic open space. The method is validated with radially bounded two-layer and three-layer waveguides and radially unbounded three-layer waveguides. Comparisons with the finite element method (FEM) demonstrate that our SEM approach requires significantly fewer degrees of freedom than FEM.
The betatron radiation source features a micrometer-scale source size, a femtosecond-scale pulse duration, milliradian-level divergence angles and a broad spectrum exceeding tens of keV. It is conducive to the high-contrast imaging of minute structures and for investigating interdisciplinary ultrafast processes. In this study, we present a betatron X-ray source derived from a high-charge, high-energy electron beam through a laser wakefield accelerator driven by the 1 PW/0.1 Hz laser system at the Shanghai Superintense Ultrafast Laser Facility (SULF). The critical energy of the betatron X-ray source is 22 ± 5 keV. The maximum X-ray flux reaches up to 4 × 109 photons for each shot in the spectral range of 5–30 keV. Correspondingly, the experiment demonstrates a peak brightness of 1.0 × 1023 photons·s−1·mm−2·mrad−2·0.1%BW−1, comparable to those demonstrated by third-generation synchrotron light sources. In addition, the imaging capability of the betatron X-ray source is validated. This study lays the foundation for future imaging applications.
Traditional wavefront control in high-energy, high-intensity laser systems usually lacks real-time capability, failing to address dynamic aberrations. This limits experimental accuracy due to shot-to-shot fluctuations and necessitates long cool-down phases to mitigate thermal effects, particularly as higher repetition rates become essential, for example, in inertial fusion research. This paper details the development and implementation of a real-time capable adaptive optics system at the Apollon laser facility. Inspired by astronomical adaptive optics, the system uses a fiber-coupled 905 nm laser diode as a pilot beam that allows for spectral separation, bypassing the constraints of pulsed lasers. A graphics processing unit-based controller, built on the open-source Compute And Control for Adaptive Optics framework, manages a loop comprising a bimorph deformable mirror and a high-speed Shack–Hartmann sensor. Initial tests showed excellent stability and effective aberration correction. However, integration into the Apollon laser revealed critical challenges unique to the laser environment that must be resolved to ensure safe operation with amplified shots.
In this work, the stability and transition to turbulence over blunt flat plates with different leading-edge radii are investigated computationally. The benchmark experimental work for comparative studies is conducted by Borovoy et al. (AIAA J., vol. 60, 2022, pp. 497–507). The freestream Mach number is 5, the unit Reynolds number is $6\times 10^7$ m$^{-1}$, and the maximum nose-tip radius 3 mm exceeds the experimental reversal value. High-resolution numerical simulation and stability analysis are performed. Three-dimensional broadband perturbation is added on the far field boundary to initiate the transition. The highlight of this work is that the complete physical process is considered, including the three-dimensional receptivity, linear and nonlinear instabilities, and transition. The experimental reversal phenomenon is reproduced favourably in the numerical simulation for the first time. Linear stability analysis shows that unstable first and second modes are absent in the blunt-plate flows owing to the presence of the entropy layer, although these modes are evident in the sharp-leading-edge case. Therefore, the transition on the blunt plate is due to non-modal instabilities. Numerical results for all the blunt-plate cases reveal the formation of streamwise streaky structures downstream of the nose (stage I) and then the presence of intermittent turbulent spots in the transitional region (stage II). In stage I, a preferential spanwise wavelength approximately 0.9 mm is selected for all the nose-tip radii, and low-frequency components are dominant. In stage II, high-frequency secondary instabilities appear to grow, which participate in the eventual breakdown. By contrast, leading-edge streaks are not remarkable in the sharp-leading-edge case, where transition is induced by oblique first and Mack second modes. The transition reversal beyond the critical nose-tip radius arises from an increasing magnitude of the streaky response in the early stage, while the transition mechanism stays similar qualitatively.