Moving Boundary Research Articles

Stimulated Brillouin Scattering Gain of Waveguides Calculated with Acoustic Perturbation Method

Abstract Stimulated Brillouin Scattering (SBS) is a phenomenon of energy transfer from an optical pump beam to longer wavelengths of light through its interaction with the medium via acoustic phonons. Previous studies have reported a calculation of the SBS gain based on the optical force approach, where the gain is not only affected by the classical paradigm of intrinsic material’s nonlinear effects in the form of electrostriction but also due to radiation pressure. In this work, the acoustic perturbation approach is applied in analyzing the optical and mechanical response of the waveguide, wherein calculating the scattering process, the changes of the waveguide structure due to photoelasticity (PE), termed Moving Boundary (MB), are considered. This acoustic perturbation method avoids uncertainties related to optical force contributions. To validate its applicability, the present method is used to calculate the SBS gain of various waveguide structures, namely an unsuspended Silicon nano-waveguide, ridge Lithium Niobate (LiNBO3) nano-waveguide on Sapphire (Al2O3) substrate and buried Arsenic trisulfide (As2S3) in Silicon dioxide (SiO2). The forward and backward SBS gain obtained using the present method of acoustic perturbation are similar to the reported values obtained from other methods of calculation as well as experiments.

Dynamic Contact Analysis of Flexible Telescopic Boom Systems with Moving Boundary

Dynamic Contact Analysis of Flexible Telescopic Boom Systems with Moving Boundary

Oxygen Switchable Photo-Hydrovoltaic Effect along the Silicon-Water Interface.

Moving boundaries of electrical double layers have shown promising capability in driving directional electron flows in solids, leading to a range of hydrovoltaic effects. The recent discovery of a photohydrovoltaic phenomenon utilizes a moving illumination zone to generate moving boundaries with different properties at the solid-water interface, referred to as the kinetic photovoltaic effect. Here, oxygen was found to act as a chemical switch to turn on and off the kinetic photovoltaic effect. Introducing oxygen would rapidly diminish the kinetic photovoltage in p-Si. On the contrary, degassing oxygen leads to a gradual recovery, whose rate can be facilely speeded up by more than one order through electrostatic gating. Mechanistic investigations of the oxygen switch behavior uncovered a dependence of surface band bending intensity of silicon on oxygen adsorption, which highlights the role of gas molecules, often overlooked, in applications based on semiconductor-liquid interfaces, such as photoelectrochemistry.

An Inverse Method for Solving Problems about Oscillations of Mechanical Systems with Moving Boundaries

An Inverse Method for Solving Problems about Oscillations of Mechanical Systems with Moving Boundaries

Integro-Interpolation Method of Constructing a Difference Scheme in a Problem with a Moving Boundary

Integro-Interpolation Method of Constructing a Difference Scheme in a Problem with a Moving Boundary

Day ahead demand response model with algorithm-based consumption classification and tariff planning

Demand Response (DR) programs effectively support the supply-demand balance in power systems. This study proposes a new Time-of-Use (ToU) tariff model, a type of price-based demand response (DR). The proposed model aims to shape the daily consumption profile of the distribution region by encouraging consumers through tariffs. Generally, in ToU demand response, various presumptions are being adopted to define consumption periods and consumer classes. It is unlikely that these presumptions are suitable for every day of the year and every consumer in the system. Achieving these definitions and classifications based on algorithms, free from human error, will increase the success and comprehensiveness of the DR program. To maximize success in this context, a deep learning model has been leveraged for day-ahead consumption forecasting. The forecasted daily consumption profile is partitioned into ToUs by using the Moving Boundary method. In addition, consumers within the region have been segmented based on the similarities in consumption behavior. Tariffs have been priced for each cluster through an optimization model that considers their participation rates in total consumption and contributions to consumption in ToU regions. The proposed model has been tested over a dataset belonging to a distribution system where different types of consumers coexist. Two sample days have been chosen as one is a weekday and the other is a weekend. The results prove that shifting consumption as encouraged by the proposed TOU method reduces consumers' costs. While, the profitability of the system operator is also preserved.

LANDAU LEGENDRE WAVELET GALERKIN METHOD APPLIED TO STUDY TWO-PHASE MOVING BOUNDARY PROBLEM OF HEAT TRANSFER IN FINITE REGION

In this paper, we developed a mathematical model of solidification where specific heat and thermal conductivity are temperature dependent. This model is a two-phase moving boundary problem (MBP) of heat transfer in finite region and represents as MBP of system of parabolic nonlinear second-order partial differential equations (PDEs). We developed a Landau Legendre wavelet Galerkin method for finding the solution of the problem. The MBP of a system of PDEs is transformed into a variable boundary value problem of nonlinear ordinary defferential equations (ODEs) by the use of dimensionless variables and the Landau transform. The problem is converted into a system of algebraic equations with the application of Legendre wavelet Galerkin method. In particular case, we compared present solution with Laplace transform solution and found approximately the same. The whole investigation has been done in dimensionless form. When the specific heat and thermal conductivity exponentially vary in temperatures, it is the effect of dimensionless parameters: thermal diffusivity (α<sub>12</sub>), ratio of thermal conductivity (<i>k</i><sub>12</sub>), dimensionless temperature (θ<sub><i> f</i></sub>), Fourier number (F<sub>0</sub>), Stefan number (Ste), and ratio of densities (ρ<sub>1</sub> / ρ<sub>2</sub>) are discussed in detail.

The dynamic behaviour of ORC-based power units fed by exhaust gases of internal combustion engines in mobile applications

Exhaust gases Waste Heat Recovery (WHR) in Internal Combustion Engines (ICEs) performed through ORC-based Power Units (ORC-bPU) ensures a reduction of fuel consumption and consequently CO2 emissions. Nevertheless, high fluctuations of hot source push the ORC-bPU to work frequently in severe off-design conditions. Hence, a deep knowledge of the unit dynamic behaviour is fundamental to develop control and operating strategies allowing to achieve a successful recovery. Dynamic response assessment of ORC-bPUs has been treated in literature but making mainly reference to theoretical approaches or to variations of input variables without considering the matching to the engine and the necessary constraints which are present in a real mobile application. Moreover, Moving Boundary (MB) models are generally used as theoretical tools to support the dynamic analysis focusing the attention on Heat Recovery Vapor Generator (HRVG) which rules the plant dynamic response. In MB models, boundary conditions at the HRVG exit are mainly fixed by the expander and by the remaining circuit. To date, the attention was focused on the presence as expander of dynamic machines (turbine), deserving an evident lower attention to volumetric expanders nonetheless their higher suitability for low-medium scale ORC applications (few kWs). To fill these knowledge gaps, a novel approach is pursued in the present paper. The dynamic response of an ORC-bPU was experimentally observed for a step change of torque and revolution speed of an ICE directly connected to the recovery unit. The variation of engine operating conditions produced an unsteady variation of the temperature and flow rate of the exhaust gas feeding the unit. The dynamics of the engine and the connecting pipes is considered as well as all the typical constraints of a real situation. Results show HRVG and exhaust pipe ruling the dynamic behaviour of the unit, leading to long time constant and raising time of gas temperature of 150 s and 300 s respectively. These results show that the unit almost always works in transient conditions being engine torque and speed variations very fast. Nevertheless, thanks to the self-regulation capacity and the robustness of the scroll expander, the unit reacts in a satisfactory way. A further novelty of the present study concerns the integration of the widely used Moving Boundary (MB) approach with Permeability theory which considers the presence of a volumetric expander and all the ither downstream components, finding a match not available in literature between the HRVG and a volumetric expander and pump. The results of the new Permeability-Moving Boundary model (P-MB) was compared with those of a more comprehensive mathematical model. The comparisons show a comparable accuracy and a significantly lower computational time achieved with P-MB. Hence, P-MB could be implemented on board to set up control strategies enhancing the recovery tuned by a model-based approach.

Numerical investigation and optimization on dynamic power generation performance of enhanced geothermal system

The Enhanced geothermal system (EGS) is a crucial technology for the development and utilization of hot dry rock (HDR). However, current research typically overlooks the strong coupling characteristics of fracture seepage, wellbore transport, and thermoelectric conversion, leading to a lack of precision in describing the dynamic operational performance of EGS throughout its entire lifecycle. A dynamic numerical model of evaporator and condenser is established by the moving boundary (MB) method. Subsequently, the dynamic power generation process is integrated with a coupled thermal-hydraulic-mechanical (THM) model of the fractured reservoir, proposing a methodology for calculating the dynamic performance of EGS throughout its lifecycle.Traditional EGS models exhibit significant deviations between calculated results and actual performance due to the neglect of instantaneous changes in injection temperature and the failure to consider the nonlinear characteristics of heat exchangers. The production pressure difference is a primary factor influencing the net output power of EGS. Therefore, selecting an appropriate injection pressure is crucial to ensuring continuous power generation throughout the EGS lifecycle. By optimizing the operating parameters of the actual EGS project, the optimal injection flow rate, fracture aperture, initial fracture permeability and well spacing are 50 kg/s, 0.03 m, 1.5 × 10−10 m2 and 350 m respectively.

Experimental Validation of Switched Moving Boundary Modeling of Phase Change Thermal Energy Storage Systems

Abstract Thermal Energy Storage (TES) devices use Phase Change Materials (PCMs) to store and release thermal energy. Control-oriented models are needed to predict the behavior of TES devices and experimental validation is necessary to demonstrate the predictive capabilities of these models. This paper presents an experimental validation of a switched Moving Boundary (MB) approach for modeling TES devices, where the dynamics of the device are captured with fewer states than traditional models. A graph-based modeling approach is used to model heat flow, while the moving boundary captures the time-varying liquid and solid regions of the TES. The model uses a Finite State Machine (FSM) to switch between four modes of operation based on the State-of-Charge (SOC) of the TES. Results show that the switched MB approach has similar accuracy and lower computational cost compared to traditional modeling approaches when predicting the SOC of an experimental TES device.

The Effect of a Moving Boundary on the Shear Strength of Granular Materials in a Direct Shear Test

The boundary state significantly influences the soil shear strength. Therefore, it is necessary to overcome the limitations of existing indoor test instruments and determine the differences in the shear properties of granular materials to ensure the economic feasibility and mechanical integrity of engineering structures. In this study, the core formula for the direct shear test was derived from the static balancing analysis of the shear box, the external force on the specimen, and the internal force on the shear surface. Three loading methods were then developed by the staggered state of the upper and lower boxes: the upper box moving shear loading method (UM), the lower box moving shear loading method (LM), and the bidirectional moving shear loading method (BM). Finally, by manipulating the motion boundary, the discrete element method (DEM) was employed to simulate the shear test of granular materials. Among the three loading methods, the order of the peak shear stresses was as follows: UM > BM > LM. Moreover, the order of the sample post-peak stress uniformities was as follows: LM > BM > UM. A shear strength conversion formula was then proposed. The findings of this study promote the advancement of the shear mechanics theory of granular materials in direct shear testing and can serve as a scientific basis for the design and manufacture of shear equipment.

Passive Gate-Tunable Kinetic Photovoltage along Semiconductor-Water Interfaces.

Moving boundaries of electrical double layer at solid-liquid interface enables unprecedented persistent energy conversion and provokes a kinetic photovoltaic effect by moving an illumination region along the semiconductor-water interface. Here, we report a transistor-inspired gate modulation of kinetic photovoltage by applying a bias at the semiconductor-water interface. The kinetic photovoltage of both p-type and n-type silicon samples can be facilely switched on/off, stemming from the electrical-field-modulated surface band bending. In contrast to the function of solid-state transistors relying on external sources, passive gate modulation of the kinetic photovoltage is achieved simply by introducing a counter electrode with materials of desired electrochemical potential. This architecture provides the ability to modulate the kinetic photovoltage over three orders of magnitude and opens up a new way for self-powered optoelectronic logic devices.

The study on the impact of AAA wall motion on the hemodynamics based on 4D CT image data.

Purpose: To analyze the effect of the physiological deformation of the vessel wall on the hemodynamics in the abdominal aortic aneurysm (AAA), this paper compared the hemodynamics in AAA based on the moving boundary (MB) simulation and the rigid wall (RW) simulation. Method: Patient-specific models were reconstructed to generate mesh based on four-dimensional computed tomography angiography (4D CT) data. The dynamic mesh technique was used to achieve deformation of the vessel wall, surface mesh and volume mesh of the fluid domain were successively remeshed at each time step. Besides, another rigid wall simulation was performed. Hemodynamics obtained from these two simulations were compared. Results: Flow field and wall shear stress (WSS) distribution are similar. When using the moving boundary method (MBM), mean time-averaged wall shear stress (TAWSS) is lower, mean oscillatory shear index (OSI) and mean relative residence time (RRT) are higher. When using the 10th and 20th percentile values for TAWSS and 80th and 90th percentile values for RRT, the ratios of areas with low TAWSS, high OSI and high RRT to the entire vessel wall are higher than those assuming the vessel as rigid. In addition, one overlapping region of low TAWSS, high OSI and high RRT by using the MBM is consistent with the location of thrombus obtained from the follow-up imaging data. Conclusion: The hemodynamics results by using the MBM reflect a higher blood retention effect. This paper presents a potential tool to assess the risk of intraluminal thrombus (ILT) formation based on the MBM.

Equivalent Formulations of the Oxygen Depletion Problem, Other Implicit Moving Boundary Value Problems, and Implications for Numerical Approximation

Equivalent Formulations of the Oxygen Depletion Problem, Other Implicit Moving Boundary Value Problems, and Implications for Numerical Approximation

The Method of Integral Transformations for Solving Boundary-Value Problems for the Heat Conduction Equation in Limited Areas Containing a Moving Boundary

The Method of Integral Transformations for Solving Boundary-Value Problems for the Heat Conduction Equation in Limited Areas Containing a Moving Boundary

Unstructured, Pressure-Based Algorithm Using Implicit Runge–Kutta for Moving Boundary Problems

An unstructured, pressure-based algorithm, IRK-SIMPLER, is developed for simulating two-dimensional incompressible flows over moving bodies. The Navier–Stokes equations in the arbitrary Lagrangian Eulerian form are discretized using collocated, median-dual control volumes with a flux correction method. The semidiscrete form is integrated in time using a diagonally implicit Runge–Kutta method. The unstructured discretization of the spatial computational domain is composed of a circular body-fitted grid, a rectangular patch grid, and a background grid. The motion of the bodies is decomposed into a series of rotations and translations. Rotation is achieved by rotating the body-fitted grid along with the body and flipping the edges of the appropriate cells. Translation is achieved by morphing the patch grid as the body moves. For large translational displacements, the rectangular patch grid is moved with the body and rematched with the background grid to obtain a new conforming grid. The algorithm is validated for multiple incompressible flow cases.

Modified Variational Iteration Method for Solving Sedimentary Ocean Basin Moving Boundary Value Problem

The motion of the shoreline in sedimentary ocean basin were investigated and solved using the modified variational iteration method since it provides an efficient approximate solution of the problem under consideration, which appears to be accurate when the consecutive error between iterative solutions of the problem is evaluated. A diffusion equation with specified beginning and boundary conditions, with a moving interface type reflecting sediment transportation, is used to simulate the physical phenomena of sand particle migration in sea water. This article also discusses the shifting boundary value problem that occurs from shoreline change in a sedimentary ocean basin.

Fractional Variational iteration Method for Solving Fractional Order Shoreline Moving Boundary Value Problem

In this article, the mathematical model of shoreline problem will be presented, which includes Caputo sense space-time fractional order derivative. The solution to this problem would be an approximate solution using fractional variational iteration method (FVIM). Graphically the results will be compared with an exact solution to the integer order derivatives. Three cases are discussed, time fractional, space fractional and standard diffusion problem.

Travelling-Wave and Asymptotic Analysis of a Multiphase Moving Boundary Model for Engineered Tissue Growth

We derive a multiphase, moving boundary model to represent the development of tissue in vitro in a porous tissue engineering scaffold. We consider a cell, extra-cellular liquid and a rigid scaffold phase, and adopt Darcy’s law to relate the velocity of the cell and liquid phases to their respective pressures. Cell–cell and cell–scaffold interactions which can drive cellular motion are accounted for by utilising relevant constitutive assumptions for the pressure in the cell phase. We reduce the model to a nonlinear reaction–diffusion equation for the cell phase, coupled to a moving boundary condition for the tissue edge, the diffusivity being dependent on the cell and scaffold volume fractions, cell and liquid viscosities and parameters that relate to cellular motion. Numerical simulations reveal that the reduced model admits three regimes for the evolution of the tissue edge at large time: linear, logarithmic and stationary. Employing travelling-wave and asymptotic analysis, we characterise these regimes in terms of parameters related to cellular production and motion. The results of our investigation allow us to suggest optimal values for the governing parameters, so as to stimulate tissue growth in an engineering scaffold.

Aeroelastic Topology Optimization of Wing Structure Based on Moving Boundary Meshfree Method

The increasing structural flexibility of large aircraft leads to significant aeroelastic effects. More efficient topology optimization techniques are required for the design to further take advantage of aeroelasticity and obtain lightweight structures. This paper proposes a moving boundary meshfree topology optimization that combines the Galerkin method of weighted residuals and non-uniform rational B-splines (NURBS). The solution domain is described by the control points of NURBS and its property is calculated adaptively with an integration subtraction technique. The minimal compliance is searched for using the globally convergent method of moving asymptotes (GCMMA) by designing the locations of control points as subject to volume and flux constraints. The method is first applied to a typical two-dimensional design example with symmetric boundary conditions. The results show that the shape constraints can be conveniently applied, and smoother boundaries are obtained with fewer parameters. Then, a three-dimensional wing structure with asymmetric boundary conditions is optimized. A three-dimensional flight load that combines the high-order-panel and meshfree methods is employed to calculate the elastic loads and update asymmetric external loads during the optimization process. The designed wing satisfies engineering requirements and the presented method can solve the practical topology optimization problems of three-dimensional structures.