These are the test problems reported in the FEBio paper that introduced the fluid-solutes solver.
Owner: ateshian
Jay J. Shim, Steve A. Maas, Jeffrey A. Weiss, & Gerard A. Ateshian (2023). Finite Element Implementation Of Computational Fluid Dynamics With Reactive Neutral And Charged Solute Transport In FEBio. Journal of Biomechanical Engineering, (), 1-26. DOI: 10.1115/1.4062594
fluid-solutes
Files
Filename
Description
Tags
Download
Fickian Diffusion/FickianDiffusion.fs2
This model analyzes axisymmetric diffusion of a neutral solute into a disk, from the top and lateral boundaries. It models the solvent as inviscid. It uses natural boundary conditions on boundaries that are impermeable to the solvent and solute, and prescribed concentration on the top and lateral boundaries.
fluid-solutes
Solute Convection in Planar Poiseuille Flow/ConvectionSolutePoiseuilleFMPe1e2.fs2
This problem demonstrates solute convection due to fully-developed Poiseuille flow in a two-dimensional domain. The Peclet number for this flow is 100, implying significant solute diffusion along with convection. For this relatively low Peclet number it is possible to use the fully-coupled solution method.
fluid-solutes
Solute Convection in Planar Poiseuille Flow/ConvectionSolutePoiseuilleFMsequentialPe1e4.fs2
This problem demonstrates solute convection due to fully-developed Poiseuille flow in a two-dimensional domain. The Peclet number for this flow is 1e4, implying insignificant solute diffusion during convection. For this elevated Peclet number the analysis uses the sequential solution method with the concentration tolerance set to 0, to solve for solute concentration in a single step at each time point.
fluid-solutes
Solute Convection in Planar Poiseuille Flow/ConvectionSolutePoiseuilleFMsequentialPe1e8.fs2
This problem demonstrates solute convection due to fully-developed Poiseuille flow in a two-dimensional domain. The Peclet number for this flow is 1e8, implying insignificant solute diffusion during convection. For this elevated Peclet number the analysis uses the sequential solution method with the concentration tolerance set to 0, to solve for solute concentration in a single step at each time point.
fluid-solutes
Flow Past Block in Narrow Channel/FlowPastBlockPe1e3.fs2
This problem demonstrates solute convection for flow past a square block in a narrow channel, in two dimensions. The Peclet number for this flow is 1000. For this relatively elevated Peclet number the solver uses use the fully-coupled method, but with the convergence tolerance on the solute concentration set to ctol = 0. Because the fluid flow experiences a backflow at the outlet boundary, this problem uses solute backflow stabilization.
fluid-solutes
Flow Past Block in Narrow Channel/FlowPastBlockPe1e5.fs2
This problem demonstrates solute convection for flow past a square block in a narrow channel, in two dimensions. The Peclet number for this flow is 1e5. For this elevated Peclet number the solver uses use the fully-coupled method, but with the convergence tolerance on the solute concentration set to ctol = 0. Because the fluid flow experiences a backflow at the outlet boundary, this problem uses solute backflow stabilization.
fluid-solutes
Flow Past Block in Narrow Channel/FlowPastBlockPe1e5extended.fs2
This problem demonstrates solute convection for flow past a square block in a narrow channel, in two dimensions, using an extended domain so that the flow does not reach the outlet boundary over the requested solution time. The Peclet number for this flow is 1e5. For this elevated Peclet number the solver uses use the sequential method, with the convergence tolerance on the solute concentration set to ctol = 0 to analyze the solute response in a single step at each time point.
fluid-solutes
Flow Past Block in Narrow Channel/FlowPastBlockPe1e5sequential.fs2
This problem demonstrates solute convection for flow past a square block in a narrow channel, in two dimensions. The Peclet number for this flow is 1e5. For this elevated Peclet number the solver uses use the sequential method, with the convergence tolerance on the solute concentration set to ctol = 0 to analyze the solute response in a single step at each time point. Because the fluid flow experiences a backflow at the outlet boundary, this problem uses solute backflow stabilization.
fluid-solutes
Solute Transport through Carotid Bifurcation/BifurcatedCarotidArteryPe3e6sequential.fs2
This problem analyzes solute convection by blood flowing through a carotid bifurcation, at a Peclet number that peaks at 3e6. Due to this elevated Peclet number the model uses the sequential solver. Solute backflow stabilization is prescribed at the outlet boundaries. To prevent non-physical accumulation of solute concentration at stagnation points near the bifurcation, the solute concentration is set to zero on the no-slip lateral boundary of the artery. This model analyzes a single cycle of arterial blood flow, assuming that the outlet fluid pressure is zero.
fluid-solutes
Solute Transport through Carotid Bifurcation/BifurcatedCarotidArteryPe3e6RCRsequential.fs2
This problem analyzes solute convection by blood flowing through a carotid bifurcation, at a Peclet number that peaks at 3e6. Due to this elevated Peclet number the model uses the sequential solver. Solute backflow stabilization is prescribed at the outlet boundaries. To prevent non-physical accumulation of solute concentration at stagnation points near the bifurcation, the solute concentration is set to zero on the no-slip lateral boundary of the artery. This model analyzes six cycles of arterial blood flow, assuming that the outlet fluid pressure is regulated by a Windkessel (RCR circuit) outlet condition.
fluid-solutes
Salt Dissociation/SaltDissociationForward_k_1.fs2
This simple model analyzes a forward chemical reaction modeling salt dissociation in a fluid-solutes domain. Three solutes are included: NaCl (electrically neutral), Na+ (positive charge) and Cl- (negative charge) and the reaction dissociates NaCl into Na+ and Cl-. This model uses a forward reaction rate kF = 1.
fluid-solutes, chemical reaction
Salt Dissociation/SaltDissociationForward_k_10.fs2
This simple model analyzes a forward chemical reaction modeling salt dissociation in a fluid-solutes domain. Three solutes are included: NaCl (electrically neutral), Na+ (positive charge) and Cl- (negative charge) and the reaction dissociates NaCl into Na+ and Cl-. This model uses a forward reaction rate kF = 10.
fluid-solutes, chemical reaction
Salt Dissociation/SaltDissociationForward_k_p1.fs2
This simple model analyzes a forward chemical reaction modeling salt dissociation in a fluid-solutes domain. Three solutes are included: NaCl (electrically neutral), Na+ (positive charge) and Cl- (negative charge) and the reaction dissociates NaCl into Na+ and Cl-. This model uses a forward reaction rate kF = 0.1.
fluid-solutes, chemical reaction
Salt Dissociation/SaltDissociationReversible_kr_1.fs2
This simple model analyzes a reversible chemical reaction modeling salt dissociation in a fluid-solutes domain. Three solutes are included: NaCl (electrically neutral), Na+ (positive charge) and Cl- (negative charge) and the reaction dissociates NaCl into Na+ and Cl-. This model uses a forward reaction rate kF = 1 and a reverse reaction rate kR=1.
fluid-solutes, chemical reaction
Salt Dissociation/SaltDissociationReversible_kr_10.fs2
This simple model analyzes a reversible chemical reaction modeling salt dissociation in a fluid-solutes domain. Three solutes are included: NaCl (electrically neutral), Na+ (positive charge) and Cl- (negative charge) and the reaction dissociates NaCl into Na+ and Cl-. This model uses a forward reaction rate kF = 1 and a reverse reaction rate kR=10.
fluid-solutes, chemical reaction
Salt Dissociation/SaltDissociationReversible_kr_p1.fs2
This simple model analyzes a reversible chemical reaction modeling salt dissociation in a fluid-solutes domain. Three solutes are included: NaCl (electrically neutral), Na+ (positive charge) and Cl- (negative charge) and the reaction dissociates NaCl into Na+ and Cl-. This model uses a forward reaction rate kF = 1 and a reverse reaction rate kR=0.1.
fluid-solutes, chemical reaction
Current Flow in an Electrolyte/CurrentFlowElectrolyte.fs2
This model illustrates one-dimensional current flow in an electrolyte consisting of a mixture of a fluid solvent, and two solutes, Na+ and Cl-. It uses the effective solute flux as a prescribed boundary condition. To ground the domain electrically, the effective solute concentrations of Na+ and Cl- are fixed at the center of the one-dimensional domain.
fluid-solutes
Expiration of Different-Sized Saliva Droplets/ExpirationSalivaDroplets_v10_sequential_R1.fs2
This example is a simplified two-dimensional simulation of a person talking, breathing, or coughing where saliva droplets are exhaled into a larger domain, to analyze the flow of these droplets and their effective sedimentation time toward the ground due to gravity, given that expiration produces flow vortices. This model uses a relatively coarser mesh for saliva droplets with a radius of 1 micron. The Peclet number for this model peaks at 4e9, therefore the sequential solver is used.
fluid-solutes
Expiration of Different-Sized Saliva Droplets/ExpirationSalivaDroplets_v10_sequential_R15.fs2
This example is a simplified two-dimensional simulation of a person talking, breathing, or coughing where saliva droplets are exhaled into a larger domain, to analyze the flow of these droplets and their effective sedimentation time toward the ground due to gravity, given that expiration produces flow vortices. This model uses a relatively coarser mesh for saliva droplets with a radius of 15 microns. The Peclet number for this model peaks at 6e10, therefore the sequential solver is used.
fluid-solutes
Expiration of Different-Sized Saliva Droplets/ExpirationSalivaDroplets_v10_sequential_R30.fs2
This example is a simplified two-dimensional simulation of a person talking, breathing, or coughing where saliva droplets are exhaled into a larger domain, to analyze the flow of these droplets and their effective sedimentation time toward the ground due to gravity, given that expiration produces flow vortices. This model uses a relatively coarser mesh for saliva droplets with a radius of 30 micron. The Peclet number for this model peaks at 1e11, therefore the sequential solver is used.
fluid-solutes
Expiration of Different-Sized Saliva Droplets/ExpirationSalivaDropletsFM_v10_sequential_R1.fs2
This example is a simplified two-dimensional simulation of a person talking, breathing, or coughing where saliva droplets are exhaled into a larger domain, to analyze the flow of these droplets and their effective sedimentation time toward the ground due to gravity, given that expiration produces flow vortices. This model uses a relatively finer mesh for saliva droplets with a radius of 1 micron. The Peclet number for this model peaks at 4e9, therefore the sequential solver is used.
fluid-solutes
Expiration of Different-Sized Saliva Droplets/ExpirationSalivaDropletsFM_v10_sequential_R15.fs2
This example is a simplified two-dimensional simulation of a person talking, breathing, or coughing where saliva droplets are exhaled into a larger domain, to analyze the flow of these droplets and their effective sedimentation time toward the ground due to gravity, given that expiration produces flow vortices. This model uses a relatively finer mesh for saliva droplets with a radius of 15 microns. The Peclet number for this model peaks at 6e10, therefore the sequential solver is used.
fluid-solutes
Expiration of Different-Sized Saliva Droplets/ExpirationSalivaDropletsFM_v10_sequential_R30.fs2
This example is a simplified two-dimensional simulation of a person talking, breathing, or coughing where saliva droplets are exhaled into a larger domain, to analyze the flow of these droplets and their effective sedimentation time toward the ground due to gravity, given that expiration produces flow vortices. This model uses a relatively finer mesh for saliva droplets with a radius of 30 micron. The Peclet number for this model peaks at 1e11, therefore the sequential solver is used.
FEBio Studio Model Repository