A multiscale 3D finite element analysis of fluid/solute transport in mechanically loaded bone

The transport of fluid, nutrients, and signaling molecules in the bone lacunar-canalicular system （LCS） is critical for osteocyte survival and function. We have applied the fluorescence recovery after photobleaching （FRAP） approach to quantify load-induced fluid and solute transport in the LCS in situ, but the measurements were limited to cortical regions 30-50 μm underneath the periosteum due to the constrains of laser penetration. With this work, we aimed to expand our understanding of load-induced fluid and solute transport in both trabecular and cortical bone using a multiscaled image-based finite element analysis （FEA） approach. An intact murine tibia was first re-constructed from microCT images into a three-dimensional （3D） linear elastic FEA model, and the matrix deformations at various locations were calculated under axial loading. A segment of the above 3D model was then imported to the biphasic poroelasticity analysis platform （FEBio） to predict load-induced fluid pressure fields, and interstitial solute/fluid flows through LCS in both cortical and trabecular regions. Further, secondary flow effects such as the shear stress and/or drag force acting on osteocytes, the presumed mechano-sensors in bone, were derived using the previously developed ultrastructural model of Brinkman flow in the canaliculi. The material properties assumed in the FEA models were validated against previously obtained strain and FRAP transport data measured on the cortical cortex. Our results demonstrated the feasibility of this computational approach in estimating the fluid flux in the LCS and the cellular stimulation forces （shear and drag forces） for osteocytes in any cortical and trabecular bone locations, allowing further studies of how the activation of osteocytes correlates with in vivo functional bone formation. The study provides a promising platform to reveal potential cellular mechanisms underlying the anabolic power of exercises and physical activities in treating patients with skeletal deficiencies.

Solute transport and geochemical modeling of the coastal quaternary aquifer, Delta Dahab Basin, South Sinai, Egypt

The wadi dahab delta is in a dry, arid coastal zone within Egypt’s south Sinai Peninsula’s eastern portion. The primary water source is the Quaternary coastal alluvial aquifer. The groundwater salinity varies from 890to 8213 mg/L, with a mean value of 3417 mg/L. The dissolved major ions have been used to calculate the seawater mixing index(SWMI) using a linear equation that discriminates the groundwater mostly affected by water–rock interaction(SWMI 1>) and other samples mixed with Seawater(SWMI < 1). The isotopic composition of groundwater for specifically chosen groundwater samples ranges from-0.645‰ to +5.212‰ for δ^(18)O and from-9.582‰ to + 22.778‰ for δ^(2)H, where the seawater represented by a Red Sea water sample(δ^(18)O + 1.64‰-δ^(2)H + 9.80‰) and reject brine water are considerably enriched the isotopic groundwater values. The geochemical NETPATH model constrained by the dissolved significant ions, isotopes, and the rock aquifer forming minerals as phases indicate the mixing percent with the seawater ranges from 9% to 97% of seawater from 91% to 3% of original recharge water. According to the SEAWAT 3-D flow models, seawater has penetrated the Northeastern Dahab delta aquifer, with the intrusion zone extending1500 m inland. The salt dissolution, upwelling of saline water, recharge from the upstream mountain block, and seawater encroachment are the primary aspects contributing to the deterioration of groundwater quality. These findings may have significance for effective groundwater withdrawal management in arid locations worldwide with similar hydrogeological systems.

Exact solutions for magnetohydrodynamic nanofluids flow and heat transfer over a permeable axisymmetric radially stretching/shrinking sheet

We report on the magnetohydrodynamic impact on the axisymmetric flow of Al_(2)O_(3)/Cu nanoparticles suspended in H_(2)O past a stretched/shrinked sheet.With the use of partial differential equations and the corresponding thermophysical characteristics of nanoparticles,the physical flow process is illustrated.The resultant nonlinear system of partial differential equations is converted into a system of ordinary differential equations using the suitable similarity transformations.The transformed differential equations are solved analytically.Impacts of the magnetic parameter,solid volume fraction and stretching/shrinking parameter on momentum and temperature distribution have been analyzed and interpreted graphically.The skin friction and Nusselt number were also evaluated.In addition,existence of dual solution was deduced for the shrinking sheet and unique solution for the stretching one.Further,Al_(2)O_(3)/H_(2)O nanofluid flow has better thermal conductivity on comparing with Cu/H_(2)O nanofluid.Furthermore,it was found that the first solutions of the stream are stable and physically realizable,whereas those of the second ones are unstable.

Modeling impact of culture facilities on hydrodynamics and solute transport in marine aquaculture waters of North Yellow Sea

An increasing number of marine aquaculture facilities have been placed in shallow bays and open sea,which might significantly affect hydrodynamic and solute transport processes in marine aquaculture waters.In this study,a coupled hydrodynamic and solute transport model was developed with high-resolution schemes in marine aquaculture waters based on depth-averaged shallow water equations.A new expression of drag force was incorporated into the momentum equations to express the resistance of suspended culture cages.The coupled model was used to simulate the effect of suspended structures on tidal currents and the movement of a contaminant cloud in the marine aquaculture of the North Yellow Sea,China.The simulation results showed a low-velocity area appearing inside the aquaculture cage area,with a maximum reduction rate of velocity close to 45%under high-density culture.The results also showed that tidal currents were sensitive to the density of suspended cages,the length of cages,and the drag coefficients of cages.The transport processes of pollutants inside aquaculture facilities were inhibited away from the vicinity of the culture cage area because of the diminished tidal currents.Therefore,the suspended cages significantly affected the transport processes of pollutants in the coastal aquaculture waters.Furthermore,the reduced horizontal velocity significantly decreased the food supply for the aquaculture areas from the surrounding sea.

GLOBAL WEAK SOLUTIONS OF COMPRESSIBLE NAVIER-STOKES-LANDAU-LIFSHITZ-MAXWELL EQUATIONS FOR QUANTUM FLUIDS IN DIMENSION THREE

In this paper,we consider the weak solutions of compressible Navier-StokesLandau-Lifshitz-Maxwell(CNSLLM)system for quantum fluids with a linear density dependent viscosity in a 3D torus.By introducing the cold pressure Pc,we prove the global existence of weak solutions with the pressure P+Pc,where P=Aργwithγ≥1.Our main result extends the one in[13]on the quantum Navier-Stokes equations to the CNSLLM system.

TIME PERIODIC SOLUTIONS TO THE EVOLUTIONARY OSEEN MODEL FOR A GENERALIZED NEWTONIAN INCOMPRESSIBLE FLUID

In this paper we study a nonstationary Oseen model for a generalized Newtonian incompressible fluid with a time periodic condition and a multivalued,nonmonotone friction law.First,a variational formulation of the model is obtained;that is a nonlinear boundary hemivariational inequality of parabolic type for the velocity field.Then,an abstract first-order evolutionary hemivariational inequality in the framework of an evolution triple of spaces is investigated.Under mild assumptions,the nonemptiness and weak compactness of the set of periodic solutions to the abstract inequality are proven.Furthermore,a uniqueness theorem for the abstract inequality is established by using a monotonicity argument.Finally,we employ the theoretical results to examine the nonstationary Oseen model.

Investigation of the solute transportation coupled with heat transfer and fluid flow during twin-roll strip casting process

Mathematical model of solute [C] distribution in twin-roll strip casting process has been setup successfully with Calcosoft for the first time. Simulation result shows that in the center of the molten steel pool between the two rolls there is a vortex flow, which is a solute enriched area. But the highest solute concentration position is at the solidification front of the columnar grain zone near the cooling roll surface. Another solute enriched position is in the back flow above the nip point. Combined with the formation mechanism of microstruoture in final as cast strip, analysis shows that solute enriched area is in the transitional area between columnar and equiaxed grain zone.

GLOBAL SOLUTIONS TO 1D COMPRESSIBLE NAVIER-STOKES/ALLEN-CAHN SYSTEM WITH DENSITY-DEPENDENT VISCOSITY AND FREE-BOUNDARY

This paper is concerned with the Navier-Stokes/Allen-Cahn system,which is used to model the dynamics of immiscible two-phase flows.We consider a 1D free boundary problem and assume that the viscosity coefficient depends on the density in the form ofη(ρ)=ρ^(α).The existence of unique global H^(2m)-solutions(m∈N)to the free boundary problem is proven for when 0<α<1/4.Furthermore,we obtain the global C^(∞)-solutions if the initial data is smooth.

Effect of bubble morphology and behavior on power consumption in non-Newtonian fluids’aeration process

Due to a prolonged operation time and low mass transfer efficiency, the primary challenge in the aeration process of non-Newtonian fluids is the high energy consumption, which is closely related to the form and rate of impeller, ventilation, rheological properties and bubble morphology in the reactor. In this perspective, through optimal computational fluid dynamics models and experiments, the relationship between power consumption, volumetric mass transfer rate(kLa) and initial bubble size(d0) was constructed to establish an efficient operation mode for the aeration process of non-Newtonian fluids. It was found that reducing the d0could significantly increase the oxygen mass transfer rate, resulting in an obvious decrease in the ventilation volume and impeller speed. When d0was regulated within 2-5 mm,an optimal kLa could be achieved, and 21% of power consumption could be saved, compared to the case of bubbles with a diameter of 10 mm.

GLOBAL WEAK SOLUTIONS TO A THREE-DIMENSIONAL QUANTUM KINETIC-FLUID MODEL

In this paper,we study a quantum kinetic-fluid model in a three-dimensional torus.This model is a coupling of the Vlasov-Fokker-Planck equation and the compressible quantum Navier-Stokes equations with degenerate viscosity.We establish a global weak solution to this model for arbitrarily large initial data when the pressure takes the form p(ρ)=ργ+pc(ρ),whereγ>1 is the adiabatic coefficient and pc(ρ)satisfies■for k≥4 and some constant c>0.

Endoscopic intramural cystogastrostomy for treatment of peripancreatic fluid collection: A viewpoint from a surgeon

Percutaneous or endoscopic drainage is the initial choice for the treatment of peripancreatic fluid collection in symptomatic patients.Endoscopic transgastric fenestration(ETGF)was first reported for the management of pancreatic pseu-docysts of 20 patients in 2008.From a surgeon’s viewpoint,ETGF is a similar procedure to cystogastrostomy in that they both produce a wide outlet orifice for the drainage of fluid and necrotic debris.ETGF can be performed at least 4 wk after the initial onset of acute pancreatitis and it has a high priority over the surgical approach.However,the surgical approach usually has a better success rate because surgical cystogastrostomy has a wider outlet(>6 cm vs 2 cm)than ETGF.However,percutaneous or endoscopic drainage,ETGF,and surgical approach offer various treatment options for peripancreatic fluid collection patients based on their conditions.

An inverse analysis of fluid flow through granular media using differentiable lattice Boltzmann method

This study presents a method for the inverse analysis of fluid flow problems.The focus is put on accurately determining boundary conditions and characterizing the physical properties of granular media,such as permeability,and fluid components,like viscosity.The primary aim is to deduce either constant pressure head or pressure profiles,given the known velocity field at a steady-state flow through a conduit containing obstacles,including walls,spheres,and grains.The lattice Boltzmann method(LBM)combined with automatic differentiation(AD)(AD-LBM)is employed,with the help of the GPU-capable Taichi programming language.A lightweight tape is used to generate gradients for the entire LBM simulation,enabling end-to-end backpropagation.Our AD-LBM approach accurately estimates the boundary conditions for complex flow paths in porous media,leading to observed steady-state velocity fields and deriving macro-scale permeability and fluid viscosity.The method demonstrates significant advantages in terms of prediction accuracy and computational efficiency,making it a powerful tool for solving inverse fluid flow problems in various applications.

Solute carrier-related signature for assessing prognosis and immunity in patients with clear-cell renal cell carcinoma

Background:Clear-cell renal cell carcinoma(ccRCC)is the most common malignant kidney cancer.However,the tumor microenvironment and crosstalk involved in metabolic reprogramming in ccRCC are not well-understood.Methods:We used The Cancer Genome Atlas to obtain ccRCC transcriptome data and clinical information.The EMTAB-1980 cohort was used for external validation.The GENECARDS database contains the first 100 solute carrier(SLC)-related genes.The predictive value of SLC-related genes for ccRCC prognosis and treatment was assessed using univariate Cox regression analysis.An SLC-related predictive signature was developed through Lasso regression analysis and used to determine the risk profiles of patients with ccRCC.Patients in each cohort were separated into high-and low-risk groups based on their risk scores.The clinical importance of the signature was assessed through survival,immune microenvironment,drug sensitivity,and nomogram analyses using R software.Results:SLC25A23,SLC25A42,SLC5A1,SLC3A1,SLC25A37,SLC5A6,SLCO5A1,and SCP2 comprised the signatures of the eight SLCrelated genes.Patients with ccRCC were separated into high-and low-risk groups based on the risk value in the training and validation cohorts;the high-risk group had a significantly worse prognosis(p<0.001).The risk score was an independent predictive indicator of ccRCC in the two cohorts according to univariate and multivariate Cox regression(p<0.05).Analysis of the immune microenvironment showed that immune cell infiltration and immune checkpoint gene expression differed between the two groups(p<0.05).Drug sensitivity analysis showed that compared to the low-risk group,the high-risk group was more sensitive to sunitinib,nilotinib,JNK-inhibitor-VIII,dasatinib,bosutinib,and bortezomib(p<0.001).Survival analysis and receiver operating characteristic curves were validated using the E-MTAB-1980 cohort.Conclusions:SLC-related genes have predictive relevance in ccRCC and play roles in the immunological milieu.Our results provide insight into metabolic reprogramming in ccRCC and identify promising treatment targets for ccRCC.

Solute drag-controlled grain growth in magnesium investigated by quasi in-situ orientation mapping and level-set simulations

Critical properties of metallic materials,such as the yield stress,corrosion resistance and ductility depend on the microstructure and its grain size and size distribution.Solute atoms that favorably segregate to grain boundaries produce a pinning atmosphere that exerts a drag pressure on the boundary motion,which strongly affects the grain growth behavior during annealing.In the current work,the characteristics of grain growth in an annealed Mg-1 wt.%Mn-1 wt.%Nd magnesium alloy were investigated by advanced experimental and modeling techniques.Systematic quasi in-situ orientation mappings with a scanning electron microscope were performed to track the evolution of local and global microstructural characteristics as a function of annealing time.Solute segregation at targeted grain boundaries was measured using three-dimensional atom probe tomography.Level-set computer simulations were carried with different setups of driving forces to explore their contribution to the microstructure development with and without solute drag.The results showed that the favorable growth advantage for some grains leading to a transient stage of abnormal grain growth is controlled by several drivers with varying importance at different stages of annealing.For longer annealing times,residual dislocation density gradients between large and smaller grains are no longer important,which leads to microstructure stability due to predominant solute drag.Local fluctuations in residual dislocation energy and solute concentration near grain boundaries cause different boundary segments to migrate at different rates,which affects the average growth rate of large grains and their evolved shape.

The Conversion of Non-Dispersed Polymers into Low-Potassium Anti-Collapse Drilling Fluids

Different drillingfluid systems are designed according to mineral composition,lithology and wellbore stability of different strata.In the present study,the conversion of a non-dispersed polymer drillingfluid into a low potas-sium anti-collapsing drillingfluid is investigated.Since the two drillingfluids belong to completely different types,the key to this conversion is represented by new inhibitors,dispersants and water-loss agents by which a non-dispersed drillingfluid can be turned into a dispersed drillingfluid while ensuring wellbore stability and reason-able rheology(carrying sand—inhibiting cuttings dispersion).In particular,the(QYZ-1)inhibitors and(FSJSS-2)dispersants are used.The former can inhibit the hydration expansion capacity of clay,reduce the dynamic shear force and weaken the viscosity;the latter can improve the sealing effect and reduce thefiltrate loss.The results have shown that after adding a reasonable proportion of these substances(QYZ-1:FSJSS-2)to the non-dispersed polymer drillingfluid,while the apparent viscosity,plastic viscosity,structural viscosity andfluidity index under-went almost negligible changes,the dynamic plastic ratio increased,and thefiltration loss decreased significantly,thereby indicating good compatibility.According to the tests(conducted in the Leijia area),the density was 1.293 g/cm3,and after standing for 24 h,the SF(static settlement factor)was 0.51.Moreover,thefiltration loss was reduced to 4.0 mL,the rolling recovery rate reached 96.92%,with excellent plugging and anti-collapse performances.

Functional thermal fluids and their applications in battery thermal management:A comprehensive review

With the increasing requirements for fast charging and discharging,higher requirements have been put forward for the thermal management of power batteries.Therefore,there is an urgent need to develop efficient heat transfer fluids.As a new type of heat transfer fluids,functional thermal fluids mainly includ-ing nanofluids(NFs)and phase change fluids(PCFs),have the advantages of high heat carrying density,high heat transfer rate,and broad operational temperature range.However,challenges that hinder their practical applications remain.In this paper,we firstly overview the classification,thermophysical prop-erties,drawbacks,and corresponding modifications of functional thermal fluids.For NFs,the high ther-mal conductivity and high convective heat transfer performance were mainly elaborated,while the stability and viscosity issues were also analyzed.And then for PCFs,the high heat carrying density was mainly elaborated,while the problems of supercooling,stability,and viscosity were also analyzed.On this basis,the composite fluids combined NFs and PCFs technology,has been summarized.Furthermore,the thermal properties of traditional fluids,NFs,PCFs,and composite fluids are compared,which proves that functional thermal fluids are a good choice to replace traditional fluids as coolants.Then,battery thermal management system(BTMS)based on functional thermal fluids is summarized in detail,and the thermal management effects and pump consumption are compared with that of water-based BTMS.Finally,the current technical challenges that parameters optimization of functional thermal fluids and structures optimization of BTMS systematically are presented.In the future,it is necessary to pay more attention to using machine learning to predict thermophysical properties of functional thermal fluids and their applications for BTMS under actual vehicle conditions.

Averaged Dynamics of Fluids near the Oscillating Interface in a Hele-Shaw Cell

The steady flow in a Hele-Shaw cell filled with fluids with a high viscosity contrast in the presence of fluid oscillations is experimentally studied.The control of oscillatory dynamics of multiphase systems with interfaces is a challenging technological problem.We consider miscible(water and glycerol)and immiscible(water and high-viscosity silicone oil PMS-1000)fluids under subsonic oscillations perpendicular to the interface.Observations show that the interface shape depends on the amplitude and frequency of oscillations.The interface is undisturbed only in the absence of oscillations.Under small amplitudes,the interface between water and glycerol widens due to mixing.When the critical amplitude is reached,the interface becomes unstable to the fingering instability:Aqueous fingers penetrate the high-viscosity glycerol and induce intensive mixing of miscible fluids and associated decay of the instability.After the disappearance of the fingers,the interface takes a U-shape in the central part of the cell.A similar effect is observed for immiscible fluids:The oscillating interface tends to bend to the side of a high-viscosity fluid.Again,when the critical amplitude is reached,the fingering instability arises at the convex interface.This paper focuses on the causes of bending of the initially undisturbed interface between miscible or immiscible fluids.For this purpose,we measure the steady flow velocity near the interface and in the bulk of a high-viscosity fluid using Particle Image Velocimetry(PIV).

Localized wave solutions and interactions of the (2+1)-dimensional Hirota-Satsuma-Ito equation

This paper studies the(2+1)-dimensional Hirota-Satsuma-Ito equation.Based on an associated Hirota bilinear form,lump-type solution,two types of interaction solutions,and breather wave solution of the(2+1)-dimensional Hirota-Satsuma-Ito equation are obtained,which are all related to the seed solution of the equation.It is interesting that the rogue wave is aroused by the interaction between one-lump soliton and a pair of resonance stripe solitons,and the fusion and fission phenomena are also found in the interaction between lump solitons and one-stripe soliton.Furthermore,the breather wave solution is also obtained by reducing the two-soliton solutions.The trajectory and period of the one-order breather wave are analyzed.The corresponding dynamical characteristics are demonstrated by the graphs.

Magmatic-hydrothermal Evolution and Mineralization Mechanisms of the Wangjiazhuang Cu(-Mo)Deposit in the Zouping Volcanic Basin,Shandong Province,China:Constraints from Fluid Inclusions

The Wangjiazhuang Cu(-Mo)deposit,located within the Zouping volcanic basin in western Shandong Province,China,is unique in this area for having an economic value.In order to expound the metallogenetic characteristics of this porphyry-like hydrothermal deposit,a detailed fluid inclusion study has been conducted,employing the techniques of representative sampling,fluid inclusion petrography,microthermometry,Raman spectroscopy,LA-ICP-MS analysis of single fluid inclusions,as well as cathode fluorescence spectrometer analysis of host mineral quartz.The deposit contains mainly two types of orebodies,i.e.veinlet-dissemination-stockwork orebodies in the K-Si alteration zone and pegmatiticquartz sulfide veins above them.In addition,minor breccia ore occurs locally.Four types of fluid inclusions in the deposit and altered quartz monzonite are identified:L-type one-or two-phase aqueous inclusions,V-type vapor-rich inclusions with V/L ratios greater than 50%-90%,D-type multiphase fluid inclusions containing daughter minerals or solids and S-type silicate-bearing fluid inclusions containing mainly muscovite and biotite.Ore petrography and fluid inclusion study has revealed a three-stage mineralization process,driven by magmatic-hydrothermal fluid activity,as follows.Initially,a hydrothermal fluid,separated from the parent magma,infiltrated into the quartz monzonite,resulting in its extensive K-Si alteration,as indicated by silicate-bearing fluid inclusions trapped in altered quartz monzonite.This is followed by the early mineralization,the formation of quartz veinlets and dissemination-stockwork ores.During the main mineralization stage,due to the participation and mixing of meteoric groundwater with magmatic-sourced hydrothermal fluid,the cooling and phase separation caused deposition of metals from the hydrothermal fluids.As a result,the pegmatitic-quartz sulfide-vein ores formed.In the late mineralization stage,decreasing fluid salinity led to the formation of L-type aqueous inclusions and chalcopyrite-sulfosalt ore.Coexistence of V-type and D-type inclusions and their similar homogenization temperatures with different homogenization modes suggest that phase separation or boiling of the ore-forming fluids took place during the early and the main mineralization stages.The formation P-T conditions of S-type inclusions and the early and the main mineralization stages were estimated as ca.156-182 MPa and 450-650℃,350-450℃,18-35 MPa and 280-380℃,8-15 MPa,respectively,based on the microthermometric data of the fluid inclusions formed at the individual stages.

PDE Standardization Analysis and Solution of TypicalMechanics Problems

A numerical approach is an effective means of solving boundary value problems(BVPs).This study focuses on physical problems with general partial differential equations(PDEs).It investigates the solution approach through the standard forms of the PDE module in COMSOL.Two typical mechanics problems are exemplified:The deflection of a thin plate,which can be addressed with the dedicated finite element module,and the stress of a pure bending beamthat cannot be tackled.The procedure for the two problems regarding the three standard forms required by the PDE module is detailed.The results were in good agreement with the literature,indicating that the PDE module provides a promising means to solve complex PDEs,especially for those a dedicated finite element module has yet to be developed.